Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,964,035 T→G 7.0% K339T (AAA→ACA)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,0350TG 7.0% 1745.9 / 12.0 619K339T (AAA→ACA) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base T (310/249);  new base G (2/58);  total (312/307)
Fisher's exact test for biased strand distribution p-value = 6.53e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 630 total aligned reads displayed.

TGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACGCT  >  NC_000913/1963886‑1964183
                                                                                                                                                     |                                                                                                                                                    
 gCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGttt                                                                                                                                                     >  1:10120023/1‑150 (MQ=255)
 gCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGttt                                                                                                                                                     <  2:3777071/150‑1 (MQ=255)
 gCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGttt                                                                                                                                                     >  2:3331290/1‑150 (MQ=255)
  cGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTg                                                                                                                                                    <  1:1786949/150‑1 (MQ=255)
       aCGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTc                                                                                                                                               <  1:9344523/150‑1 (MQ=255)
         ggCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCgg                                                                                                                                             <  1:6872068/150‑1 (MQ=255)
         ggCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCgg                                                                                                                                             <  2:7459175/150‑1 (MQ=255)
          gCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGc                                                                                                                                            >  1:9158504/1‑150 (MQ=255)
            ccAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAg                                                                                                                                          <  2:3643596/150‑1 (MQ=255)
              aacaacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGgc                                                                                                                                        >  2:5342481/1‑150 (MQ=255)
               acaacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGgcg                                                                                                                                       >  2:962748/1‑150 (MQ=255)
                caacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGgcgc                                                                                                                                      <  1:7919038/150‑1 (MQ=255)
                 aacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGtt                                                                                                                                                <  1:1405526/139‑1 (MQ=255)
                 aacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCt                                                                                                                                     >  2:11049363/1‑150 (MQ=255)
                 aacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCt                                                                                                                                     >  2:7694311/1‑150 (MQ=255)
                   cTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCttt                                                                                                                                   <  1:6813852/150‑1 (MQ=255)
                     cACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACACCGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCttttt                                                                                                                                 <  2:6464998/150‑1 (MQ=255)
                      aCCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTg                                                                                                                                >  2:2107581/1‑150 (MQ=255)
                        cATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGtt                                                                                                                              <  2:4878201/150‑1 (MQ=255)
                         ctCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTc                                                                                                                             <  1:2020991/149‑1 (MQ=255)
                         aTCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTc                                                                                                                             <  1:8681635/150‑1 (MQ=255)
                               ctgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTcc                                                                                                                       <  2:6908326/150‑1 (MQ=255)
                               ctgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTcc                                                                                                                       >  1:4643963/1‑150 (MQ=255)
                               ctgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTcc                                                                                                                       >  1:5833405/1‑150 (MQ=255)
                               ctgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGGTCGGCAGGCGCTTTTTGTTCGCGTcc                                                                                                                       >  2:10160349/1‑150 (MQ=255)
                                tgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcg                                                                                                                      >  2:9291072/1‑150 (MQ=255)
                                tgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcg                                                                                                                      >  2:2361070/1‑150 (MQ=255)
                                 ggAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgc                                                                                                                     >  2:10541440/1‑150 (MQ=255)
                                 ggAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgc                                                                                                                     >  1:9963281/1‑150 (MQ=255)
                                 ggAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgc                                                                                                                     <  2:5148324/150‑1 (MQ=255)
                                  gAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgcg                                                                                                                    >  2:7972532/1‑150 (MQ=255)
                                  gAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgcg                                                                                                                    >  1:1247291/1‑150 (MQ=255)
                                  gAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgcg                                                                                                                    >  1:11042712/1‑150 (MQ=255)
                                           cATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATccacca                                                                                                           <  1:9828052/150‑1 (MQ=255)
                                            aTCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAg                                                                                                          >  2:4594121/1‑150 (MQ=255)
                                              cGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGc                                                                                                        <  1:3413539/150‑1 (MQ=255)
                                                ggAGCAGTCGATAACCCACTTCCATACCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGGTAGGCAGGCGCTTTTTGTTCGAGTCCGCGTAACCACCAGGCCa                                                                                                      >  1:3111981/1‑150 (MQ=255)
                                                  aTCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGc                                                                                                    >  1:706873/1‑150 (MQ=255)
                                                      gTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGa                                                                                                <  1:1617220/150‑1 (MQ=255)
                                                      gTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGa                                                                                                <  1:7785623/150‑1 (MQ=255)
                                                      gTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGa                                                                                                <  1:4007002/150‑1 (MQ=255)
                                                         gATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATccac                                                                                                             <  1:1497399/134‑1 (MQ=255)
                                                           tAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAAt                                                                                           <  1:7337093/150‑1 (MQ=255)
                                                               ccACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCgg                                                                                       >  1:10542374/1‑150 (MQ=255)
                                                               ccACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCgg                                                                                       >  1:9436491/1‑150 (MQ=255)
                                                               ccACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCgg                                                                                       >  2:972988/1‑150 (MQ=255)
                                                                 aCTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCa                                                                                     <  1:7536847/150‑1 (MQ=255)
                                                                 aCTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCa                                                                                     <  1:6155882/150‑1 (MQ=255)
                                                                 aCTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGGTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCCGGCCAGCCCGAGCAATCCGGCa                                                                                     >  1:8929666/1‑150 (MQ=255)
                                                                      cATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGaa                                                                                <  1:5265772/150‑1 (MQ=255)
                                                                       aTTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAAc                                                                               >  1:662018/1‑150 (MQ=255)
                                                                        ttCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACa                                                                              <  1:6285650/150‑1 (MQ=255)
                                                                        ttCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACa                                                                              <  1:3883498/150‑1 (MQ=255)
                                                                           ccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCa                                                                           >  2:1213228/1‑150 (MQ=255)
                                                                             agagAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCaaa                                                                         >  1:3664629/1‑150 (MQ=255)
                                                                             agagAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAAACCGGCAGTGAACAGCaaa                                                                         >  1:4546998/1‑150 (MQ=255)
                                                                               agaATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCaaaaa                                                                       >  2:8518705/1‑150 (MQ=255)
                                                                                   tCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATAcc                                                                   <  1:4295277/150‑1 (MQ=255)
                                                                                    cTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCa                                                                  >  1:4080589/1‑150 (MQ=255)
                                                                                     ttCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAg                                                                 >  1:6069041/1‑150 (MQ=255)
                                                                                           ttGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCgg                                                           <  1:21360/150‑1 (MQ=255)
                                                                                             gTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCa                                                         >  1:9678032/1‑150 (MQ=255)
                                                                                              tACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCAt                                                        >  1:6562639/1‑150 (MQ=255)
                                                                                              tACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCAt                                                        >  2:94834/1‑150 (MQ=255)
                                                                                                cATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTc                                                      <  2:9105036/150‑1 (MQ=255)
                                                                                                    gTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAgg                                                  <  2:1578902/150‑1 (MQ=255)
                                                                                                    gTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAgg                                                  >  1:5434067/1‑150 (MQ=255)
                                                                                                      tCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATATTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCTAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCa                                                <  2:4735560/150‑1 (MQ=255)
                                                                                                       ccACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGACCAGCAAAAATACCAGGTTCGGCATTCCAGGCAc                                               >  2:8312623/1‑150 (MQ=255)
                                                                                                         aCGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCa                                             <  2:679718/150‑1 (MQ=255)
                                                                                                         aCGTCGCTTCGACAACGGTATTATTATCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCATTTTGTTCGCGTCCGCGTATACCCCAGGCCAGCCCGAGCAATACGGCACTGTACAGCAAAAATACCAGTTTAGGCATTCCAGGCAccc                                             >  1:10246054/1‑149 (MQ=255)
                                                                                                            tCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGc                                          >  2:6153407/1‑150 (MQ=255)
                                                                                                             cGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGcc                                         <  2:7272580/150‑1 (MQ=255)
                                                                                                               cTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGa                                       >  1:4888121/1‑150 (MQ=255)
                                                                                                               cTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGa                                       <  2:4392816/150‑1 (MQ=255)
                                                                                                               cTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGa                                       >  2:590003/1‑150 (MQ=255)
                                                                                                               cTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTCCGGCATTCCAGGCACCAGGCCGa                                       <  2:5559894/150‑1 (MQ=255)
                                                                                                               cTTCGACAACGGTATTATTCTCTGCCATTTGTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGa                                       <  2:1921901/150‑1 (MQ=255)
                                                                                                                ttCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAg                                      <  1:7653938/150‑1 (MQ=255)
                                                                                                                 tCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGt                                     >  2:6371291/1‑150 (MQ=255)
                                                                                                                 tCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGt                                     >  1:8099199/1‑150 (MQ=255)
                                                                                                                  cGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTa                                    >  2:2826844/1‑150 (MQ=255)
                                                                                                                        cGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGa                              <  2:7058534/150‑1 (MQ=255)
                                                                                                                        cGGTATTATTCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGa                              <  2:10377214/150‑1 (MQ=255)
                                                                                                                           tattatTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCa                           >  1:8208153/1‑150 (MQ=255)
                                                                                                                             ttattCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCAcg                         >  2:4682358/1‑150 (MQ=255)
                                                                                                                              tattCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCAcgg                        <  2:1344427/150‑1 (MQ=255)
                                                                                                                                ttCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCAcggcg                      >  2:10693116/1‑150 (MQ=255)
                                                                                                                                ttCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCAcggcg                      >  1:7238457/1‑150 (MQ=255)
                                                                                                                                ttCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCAcggcg                      >  2:10822314/1‑150 (MQ=255)
                                                                                                                                ttCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATAACAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCAcggcg                      >  1:2354789/1‑150 (MQ=255)
                                                                                                                                     tGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCt                 >  2:4024121/1‑150 (MQ=255)
                                                                                                                                      gCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCtt                <  1:10832399/150‑1 (MQ=255)
                                                                                                                                       ccATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCGGAGTAAACCGAGCACGGCGGCGCTTa               <  2:6670196/150‑1 (MQ=255)
                                                                                                                                       ccATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCACGCCAGCCCCAGCAATCCCGCAGTGAACAGCAAAAAAACCAGGTTCGGCATTCACGGCACAAGCCCCAGTAAACCAACCCCGCCGCCGCtaa               >  1:3213289/1‑148 (MQ=255)
                                                                                                                                                ccGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCCGCCCGAGCAATCCGGCAGTGAACCGCAAAAATACCAGGTTCGGCATTCCAGGCACCCGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACaaaa      >  1:8584115/1‑150 (MQ=255)
                                                                                                                                                ccGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGGCAGCGTATACACCAGGCCAGCCCGCGCAAGACGGCAGTGTACAGCCAAAATAACAGGTTCGGCATTCCAGGCACCCGGCCGAGTGAACCGAGCACCGCGGCGCTTAACAACATaa      >  1:10077052/1‑150 (MQ=255)
                                                                                                                                                   gTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACGc   >  2:8870591/1‑150 (MQ=255)
                                                                                                                                                   gTTTGGGTTCGGCAGGCGCTTTTTGTTAGCGTCCGCCTATCCACCCCGCCAGCCCGAGCAATACGGCCGTGAACAGCAAAAATAACACGTTCGGCCATCCCGGCACCAGGCCGAGTAAACCGAGCCCCGCCGCGCGTAACAAAAAAACGc   >  1:5049405/1‑150 (MQ=255)
                                                                                                                                                    tttGGGTTCGGCAGGCGCTTTTTTTTAGCGTCCGCGTATCCACCACGCCCCCCCGAGCAAAACCGCAGTGAAAAGCAAAAAAACCAGGTTCGGCATTACAGGCACCAGGCCCAGTAAACCGAGCACGGCGGCGCTTAACAACAAAACCCt  >  1:2720645/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACGCT  >  NC_000913/1963886‑1964183

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑