Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 169,792 G→A 14.0% intergenic (+756/+2931) PP_RS08750 → / ← fhuA_1 tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143169,7920GA14.0% 178.6 / 20.7 73intergenic (+756/+2931)PP_RS08750/fhuA_1tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (5/5);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-02

AACTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGC  >  ACB143/169660‑169920
                                                                                                                                    |                                                                                                                                
aaCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGc                                                                                                                                 >  1:145154/1‑134 (MQ=32)
aaCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGc                                                                                                                                 <  2:145154/134‑1 (MQ=32)
aaCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCg                                                                                                                                >  2:261519/1‑135 (MQ=32)
aaCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGt                                                                                                                               <  6:173264/136‑1 (MQ=32)
aaCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGt                                                                                                                               <  4:35696/136‑1 (MQ=32)
    ttGTCTGAATATGGGGGGCCCATACTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTTAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAAAAGAACCTGAAACCGTGTCCGTACa                                                                                                                            >  5:150479/1‑135 (MQ=11)
    ttGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTAc                                                                                                                             >  3:142302/1‑134 (MQ=32)
    tagtgtGAATTTGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGATAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACaa                                                                                                                           >  6:309342/6‑136 (MQ=255)
    ctGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTgtacgtac                                                                                                                             >  6:248678/2‑134 (MQ=255)
          gAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTg                                                                                                                     >  6:337245/1‑136 (MQ=255)
           aaTATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTgg                                                                                                                    <  2:358729/136‑1 (MQ=255)
           aaTATGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTgg                                                                                                                    <  7:178063/135‑1 (MQ=11)
               tGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAg                                                                                                                 <  2:134317/135‑1 (MQ=32)
                  ggggACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACAGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTa                                                                                                             <  1:209021/136‑1 (MQ=17)
                          tcctccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTa                                                                                                     <  4:269401/136‑1 (MQ=32)
                          tcctccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTa                                                                                                     <  5:161022/136‑1 (MQ=32)
                                 aGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                              <  6:171198/136‑1 (MQ=32)
                                  ggCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                             <  4:302582/136‑1 (MQ=32)
                                  ggCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                             <  5:119250/136‑1 (MQ=32)
                                       aaTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTAc                                                                                        <  6:355561/136‑1 (MQ=32)
                                                gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGtat                                                                               >  4:125433/1‑136 (MQ=32)
                                                gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                >  7:96914/1‑135 (MQ=32)
                                                gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                >  7:79086/1‑135 (MQ=32)
                                                 actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGtat                                                                               >  1:118380/1‑135 (MQ=32)
                                                 actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                <  1:90901/134‑1 (MQ=32)
                                                    gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATg                                                                           >  5:152874/1‑136 (MQ=32)
                                                    gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGAGTACAAGCAGTGGGAGCATACTTTGTTAGGTGACTGCGTACCTTTTGTATAATg                                                                           >  6:234898/1‑136 (MQ=17)
                                                      ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATggg                                                                         <  1:69271/136‑1 (MQ=32)
                                                        gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                       <  2:191714/136‑1 (MQ=11)
                                                        gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGt                                                                        >  6:259957/1‑135 (MQ=11)
                                                        gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                       <  5:167885/136‑1 (MQ=11)
                                                        gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                       <  5:337245/136‑1 (MQ=11)
                                                        gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                       <  2:118380/136‑1 (MQ=32)
                                                        gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                       <  7:301697/136‑1 (MQ=32)
                                                           aGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGc                                                                    >  2:186318/1‑136 (MQ=32)
                                                                  cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtat                                                             <  4:256616/136‑1 (MQ=32)
                                                                    gTACCGTGAGGGAAAGGCGAAACGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                           <  2:368851/136‑1 (MQ=17)
                                                                       ccGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTct                                                        >  3:120985/1‑135 (MQ=11)
                                                                       ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAtt                                                          >  8:15260/1‑134 (MQ=32)
                                                                       ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAtt                                                          >  6:77657/1‑134 (MQ=32)
                                                                       ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCa                                                        <  4:307075/136‑1 (MQ=32)
                                                                       ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATggg                                                                         <  2:299074/118‑1 (MQ=17)
                                                                       ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATggg                                                                         >  1:299074/1‑118 (MQ=17)
                                                                        cGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTctg                                                       >  1:297749/1‑134 (MQ=11)
                                                                         gTGAGGGAAAGGCGAAACGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                      <  6:126195/136‑1 (MQ=17)
                                                                          tGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGAGACTTATATTCAGTg                                                     >  8:37622/1‑136 (MQ=17)
                                                                             gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                   >  7:339422/1‑135 (MQ=32)
                                                                             gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                   >  1:260444/1‑135 (MQ=32)
                                                                             gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                   >  8:312414/1‑135 (MQ=32)
                                                                              ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGa                                                 <  3:302971/136‑1 (MQ=32)
                                                                               gAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGAATTATATTCAGTGGCGAg                                                >  3:75190/1‑136 (MQ=17)
                                                                                 aaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCt                                              <  3:83289/136‑1 (MQ=32)
                                                                                 aaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCt                                              <  1:284537/136‑1 (MQ=32)
                                                                                     cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAAc                                          >  8:105830/1‑136 (MQ=32)
                                                                                      gAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTAGAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAAcc                                         >  8:364729/1‑136 (MQ=17)
                                                                                           gAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATa                                    <  1:300351/136‑1 (MQ=32)
                                                                                           gAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATa                                    <  7:241221/136‑1 (MQ=32)
                                                                                           gAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATa                                    <  2:276136/136‑1 (MQ=32)
                                                                                                 cGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAg                              >  2:192808/1‑136 (MQ=11)
                                                                                                 cGGAGAGGGGAGTGAAATAGAACCTGAACCCGTGTGCGTACACGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATTGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAg                              <  2:161752/136‑1 (MQ=12)
                                                                                                 cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAg                              <  2:225758/136‑1 (MQ=32)
                                                                                                 cGGAGAGGGGAGAGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAg                              <  8:354840/136‑1 (MQ=17)
                                                                                                          gAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                     <  6:362437/136‑1 (MQ=32)
                                                                                                          gAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                     <  8:308319/136‑1 (MQ=32)
                                                                                                          gAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGGGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                     <  3:167936/136‑1 (MQ=17)
                                                                                                              gAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGc                 >  1:257010/1‑136 (MQ=32)
                                                                                                                aaTAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGc                 >  2:205909/1‑134 (MQ=32)
                                                                                                                   aGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTc            >  2:139893/1‑136 (MQ=32)
                                                                                                                     aaCCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCt           >  5:22487/1‑135 (MQ=32)
                                                                                                                      aCCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTa         <  1:75833/136‑1 (MQ=32)
                                                                                                                             aCCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTATCGAAAGCTAGTCTTAATAGGGc  >  6:211743/1‑136 (MQ=17)
                                                                                                                             aCCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGc  <  1:223106/136‑1 (MQ=32)
                                                                                                                             aCCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGc  >  2:228386/1‑136 (MQ=32)
                                                                                                                             aCCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGc  >  4:7274/1‑136 (MQ=32)
                                                                                                                                    |                                                                                                                                
AACTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGC  >  ACB143/169660‑169920

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: