Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,964,035 T→G 5.2% 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 5.2% 1524.4 / 3.0 555K339T (AAA→ACA) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base T (258/245);  new base G (3/47);  total (261/294)
Fisher's exact test for biased strand distribution p-value = 7.28e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 572 total aligned reads displayed.

TGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACGCTT  >  NC_000913/1963886‑1964184
                                                                                                                                                     |                                                                                                                                                     
tGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGtt                                                                                                                                                       >  1:9053531/1‑150 (MQ=255)
  cGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTg                                                                                                                                                     >  2:4294420/1‑150 (MQ=255)
   gTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTTCCGGTTTgg                                                                                                                                                    <  1:6841019/150‑1 (MQ=255)
   gTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTgg                                                                                                                                                    <  1:9505928/150‑1 (MQ=255)
   gTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTgg                                                                                                                                                    <  1:7217788/150‑1 (MQ=255)
      taCGGCCCAACAACTCACCATCCTGCTTGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGTTTTGGGtt                                                                                                                                                 <  2:6472284/150‑1 (MQ=255)
      taCGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGtt                                                                                                                                                 <  1:6482472/150‑1 (MQ=255)
      taCGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGtt                                                                                                                                                 >  2:8112802/1‑150 (MQ=255)
         ggCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCgg                                                                                                                                              <  1:1028474/150‑1 (MQ=255)
           cccAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCa                                                                                                                                            >  1:9387516/1‑150 (MQ=255)
              aacaacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGgc                                                                                                                                         <  2:8392833/150‑1 (MQ=255)
               acaacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGgcg                                                                                                                                        >  2:2963007/1‑150 (MQ=255)
                caacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGgcgc                                                                                                                                       <  1:3704579/150‑1 (MQ=255)
                caacTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGgcgc                                                                                                                                       <  2:6257553/150‑1 (MQ=255)
                     cACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCttttt                                                                                                                                  <  1:1294557/150‑1 (MQ=255)
                      aCCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTg                                                                                                                                 <  2:914177/150‑1 (MQ=255)
                        cATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACCTCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGtt                                                                                                                               <  2:9066122/150‑1 (MQ=255)
                        cATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGtt                                                                                                                               <  1:764090/150‑1 (MQ=255)
                         aTCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTc                                                                                                                              <  2:4779784/150‑1 (MQ=255)
                           ccTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTcgc                                                                                                                            <  1:2449295/150‑1 (MQ=255)
                            ctgctgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTcgcg                                                                                                                           >  2:6863207/1‑150 (MQ=255)
                                tgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGACTCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTGTTACCGGGTTGGGTTCGGCAGGCGCTGTTTGTTCGCGTCcg                                                                                                                       >  1:8047141/1‑150 (MQ=255)
                                tgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcg                                                                                                                       >  2:9632307/1‑150 (MQ=255)
                                tgGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcg                                                                                                                       <  2:3087311/150‑1 (MQ=255)
                                tgGAAATCGACCATCGGGATCAGGCGATAACCCACTTCCATTCACAGAGAATCTTCCAGTTGTACATCGGTCCACGTCGCTTCGACAACGGTATTATACTCTGCCATTTTTACCGGTTAGGGTTCGGCAGGCGCTTGTTGTTCGCGTCcg                                                                                                                       >  2:4394047/1‑150 (MQ=255)
                                 ggAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgc                                                                                                                      >  1:5646197/1‑150 (MQ=255)
                                  gAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCcgcg                                                                                                                     >  1:4709115/1‑150 (MQ=255)
                                      tCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATc                                                                                                                 <  1:5154192/150‑1 (MQ=255)
                                      tCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATc                                                                                                                 >  1:5003095/1‑150 (MQ=255)
                                         aCCATCGGGATCAGTCGATAACCCACTTCCATTCCCCGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATccac                                                                                                              <  2:2815836/150‑1 (MQ=255)
                                          ccATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATccacc                                                                                                             <  2:1354083/150‑1 (MQ=255)
                                           cATCGGGATCATTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGAGTCCGCGTATccacca                                                                                                            <  2:3629815/150‑1 (MQ=255)
                                           cATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCTGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATccacca                                                                                                            <  1:8154196/150‑1 (MQ=255)
                                            aTCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAg                                                                                                           >  1:9746776/1‑150 (MQ=255)
                                            aTCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAg                                                                                                           <  1:1853332/150‑1 (MQ=255)
                                             tCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAgg                                                                                                          >  2:10053248/1‑150 (MQ=255)
                                                   tCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGcc                                                                                                    >  2:9449518/1‑150 (MQ=255)
                                                    cAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGccc                                                                                                   <  1:8831129/150‑1 (MQ=255)
                                                      gTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGa                                                                                                 <  1:8436614/150‑1 (MQ=255)
                                                      gTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGa                                                                                                 <  1:1148912/150‑1 (MQ=255)
                                                      gTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGa                                                                                                 <  2:8561438/150‑1 (MQ=255)
                                                           tAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAAt                                                                                            <  2:6098061/150‑1 (MQ=255)
                                                               ccACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTa                                                                                                                   >  2:9897145/1‑123 (MQ=255)
                                                               ccACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTa                                                                                                                   <  1:9897145/123‑1 (MQ=255)
                                                               ccACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCgg                                                                                        <  2:6012737/150‑1 (MQ=255)
                                                               ccACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCgg                                                                                        >  2:4980189/1‑150 (MQ=255)
                                                                   ttCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGt                                                                                    <  2:3926137/150‑1 (MQ=255)
                                                                   ttCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGt                                                                                    >  1:7007270/1‑150 (MQ=255)
                                                                      cATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGaa                                                                                 <  2:7731378/150‑1 (MQ=255)
                                                                       aTTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGGATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCACTCCGGCAGTGacc                                                                                >  2:489162/1‑148 (MQ=255)
                                                                        ttCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACa                                                                               <  1:6222964/150‑1 (MQ=255)
                                                                         tCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAAcgg                                                                              <  2:9754482/150‑3 (MQ=255)
                                                                           ccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCa                                                                            >  1:2786589/1‑150 (MQ=255)
                                                                             agagAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCaaa                                                                          <  2:7435235/150‑1 (MQ=255)
                                                                              gagaATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTTGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATACACCAGGCCAGCCCGAGCAATCCGGCAGGGAACAGCaaaa                                                                         >  1:4551397/1‑150 (MQ=255)
                                                                              gagaATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCaaaa                                                                         >  2:2020670/1‑150 (MQ=255)
                                                                                   tCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATAcc                                                                    <  2:1723228/150‑1 (MQ=255)
                                                                                   tCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATAcc                                                                    <  2:3462663/150‑1 (MQ=255)
                                                                                      tCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAgg                                                                 >  1:7417299/1‑150 (MQ=255)
                                                                                        cAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCCCCCCGAGCACACCGGCAGTGAACAGCAAAAATACCAGGtt                                                               >  1:1984433/1‑150 (MQ=255)
                                                                                         aGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTc                                                              >  2:7144341/1‑150 (MQ=255)
                                                                                          gTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCg                                                             >  2:8930689/1‑150 (MQ=255)
                                                                                           ttGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCgg                                                            >  1:5310689/1‑150 (MQ=255)
                                                                                           ttGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCgg                                                            >  2:6701505/1‑150 (MQ=255)
                                                                                           ttGTACATCGTTCCACGTCGCTTCGACAACGGTATTAGTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTGTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCgg                                                            >  2:8780727/1‑150 (MQ=255)
                                                                                           ttGTACATCGGTGCTCGTCGCTTCGACAACGGTATTAGTCTCTGCCATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCGGGCAGTGAACAGCAAAAATACCAGGTTCgg                                                            <  1:7848021/150‑1 (MQ=255)
                                                                                              tACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCAt                                                         <  2:8251971/150‑1 (MQ=255)
                                                                                               aCATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCCGCCCGAGCAATCCGGCAGTGACCCGCACACATACCAGGGTCGGCAtg                                                        >  1:2322971/1‑149 (MQ=255)
                                                                                                  tCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCa                                                     <  1:4375092/150‑1 (MQ=255)
                                                                                                  tCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTGGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAACAATACCAGGTTCGGCATTCCa                                                     <  2:9345333/150‑1 (MQ=255)
                                                                                                          cGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAg                                             >  2:4117453/1‑150 (MQ=255)
                                                                                                            tCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATAACAGGTTCGGCATTCCAGGCACCAGGc                                           >  1:2225771/1‑150 (MQ=255)
                                                                                                              gCTTCGACAACGGTAGTATTCTCTGCCATTTTTCCCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACATCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCg                                         <  1:2987185/150‑1 (MQ=255)
                                                                                                                ttCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAg                                       >  2:926050/1‑150 (MQ=255)
                                                                                                                ttCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAg                                       >  1:10108235/1‑150 (MQ=255)
                                                                                                                ttCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAg                                       >  1:8604262/1‑150 (MQ=255)
                                                                                                                        cGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGa                               <  1:788951/150‑1 (MQ=255)
                                                                                                                        cGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGa                               <  2:6618970/150‑1 (MQ=255)
                                                                                                                         ggTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCTGCATTCCAGGCACCAGGCCGAGTAAACCGAg                              >  2:146641/1‑150 (MQ=255)
                                                                                                                          gTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGc                             <  2:2580982/150‑1 (MQ=255)
                                                                                                                               attCTCTGCCATTTTGACCGGTGGGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCAcggc                        <  2:2228821/150‑1 (MQ=255)
                                                                                                                                    ctGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGgcgc                   <  2:458560/150‑1 (MQ=255)
                                                                                                                                       ccATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTa                >  2:7295898/1‑150 (MQ=255)
                                                                                                                                       ccATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTa                >  2:5146376/1‑150 (MQ=255)
                                                                                                                                        cATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTaa               <  1:3325165/150‑1 (MQ=255)
                                                                                                                                        cATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTaa               <  2:2372881/150‑1 (MQ=255)
                                                                                                                                        cATTTTTACCGGTGTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTaa               <  1:6015916/150‑1 (MQ=255)
                                                                                                                                         aTTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCCGGCCAGCCCCAGCAAACCGGCAGTGAACACCAAAAATACCACGTTCGGCAATCCAGGCACCACGCCGACTAAACCGAGCACCGCGGCGCGTaac              >  1:3477699/1‑150 (MQ=255)
                                                                                                                                         aTTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTaac              >  2:10262312/1‑150 (MQ=255)
                                                                                                                                          tttttACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTAGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATACCGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTaaca             >  1:3217764/1‑150 (MQ=255)
                                                                                                                                           ttttACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTaacaa            <  1:10033457/150‑1 (MQ=255)
                                                                                                                                           ttttACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTaacaa            <  2:5897566/150‑1 (MQ=255)
                                                                                                                                               cccGGTGGGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATa        <  2:5717959/149‑1 (MQ=255)
                                                                                                                                               aCCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATa        <  1:9138365/150‑1 (MQ=255)
                                                                                                                                                ccGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATaa       <  2:7926914/150‑1 (MQ=255)
                                                                                                                                                 cGGTTTGGGTTCGGCAGGCGCTTTTTTTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGGAAACCGAGCACGGCGGCGCGTAACAACAAAAc      >  1:9317339/1‑150 (MQ=255)
                                                                                                                                                 cGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAAc      <  1:3974332/150‑1 (MQ=255)
                                                                                                                                                  ggTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACg     <  2:3301995/150‑1 (MQ=255)
                                                                                                                                                  ggTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACg     <  2:541675/150‑1 (MQ=255)
                                                                                                                                                   gTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGCTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACGc    <  1:4164799/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAGGCACCAGGCCGAGTAAACCGAGCACGGCGGCGCTTAACAACATAACGCTT  >  NC_000913/1963886‑1964184

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: