Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,982,559 T→A 5.0% intergenic (‑31/+40) PP_4387 ← / ← flgE hypothetical protein/flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,982,5590TA 5.0% 656.0 / 15.5 203intergenic (‑31/+40)PP_4387/flgEhypothetical protein/flagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base T (94/95);  new base A (6/5);  total (103/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 207 total aligned reads displayed.

GAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCA  >  NC_002947/4982410‑4982708
                                                                                                                                                     |                                                                                                                                                     
gAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGTGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTAGTGGCCGCt                                                                                                                                                       >  1:4037652/1‑150 (MQ=255)
  acacAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCttt                                                                                                                                                     <  1:137945/150‑1 (MQ=255)
    acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                   >  2:2028623/1‑150 (MQ=255)
    acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                   <  1:2201078/150‑1 (MQ=255)
    acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                   <  2:4277238/150‑1 (MQ=255)
        gTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGcc                                                                                                                                               >  2:733482/1‑150 (MQ=255)
        gTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGcc                                                                                                                                               >  1:4217901/1‑150 (MQ=255)
              gCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaaga                                                                                                                                         <  1:1297437/150‑3 (MQ=255)
               ccGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaagaa                                                                                                                                        <  2:4309431/150‑4 (MQ=255)
                cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAc                                                                                                                                       >  2:3822088/1‑150 (MQ=255)
                 ggTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACa                                                                                                                                      <  1:581765/150‑1 (MQ=255)
                 ggTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACa                                                                                                                                      <  1:581839/150‑1 (MQ=255)
                 ggTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACa                                                                                                                                      >  1:804793/1‑150 (MQ=255)
                   tGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACaaa                                                                                                                                    <  1:1510820/150‑1 (MQ=255)
                       cAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACAcc                                                                                                                                >  2:434771/1‑150 (MQ=255)
                        aGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCa                                                                                                                               <  2:1569002/150‑1 (MQ=255)
                         ggCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaa                                                                                                                              >  2:3035562/1‑150 (MQ=255)
                          gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCTCCAATTACAAACACCaaa                                                                                                                             <  2:3423180/150‑1 (MQ=255)
                          gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                             <  1:1279024/150‑1 (MQ=255)
                          gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                             <  1:742359/150‑1 (MQ=255)
                            tttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                           >  1:2419096/1‑150 (MQ=255)
                             ttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGcc                                                                                                                          <  1:910003/150‑1 (MQ=255)
                               gCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                                        <  1:2597974/150‑1 (MQ=255)
                               gCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                                        <  2:444943/150‑1 (MQ=255)
                               gCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                                        <  2:3423590/150‑1 (MQ=255)
                                   ttGATAGAACGGGGCGGGGGCAGCGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACaa                                                                                                                   >  1:426031/1‑150 (MQ=255)
                                     gATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAat                                                                                                                  <  1:2060024/150‑1 (MQ=255)
                                      aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatc                                                                                                                 >  2:1052347/1‑150 (MQ=255)
                                      aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatc                                                                                                                 >  1:3473349/1‑150 (MQ=255)
                                          aaCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAtcatca                                                                                                             >  2:178266/1‑150 (MQ=255)
                                           aCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAg                                                                                                            >  1:942155/1‑150 (MQ=255)
                                             ggggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAATCAGGt                                                                                                          >  1:1551737/1‑150 (MQ=255)
                                              gggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTc                                                                                                         >  2:2141476/1‑150 (MQ=255)
                                                 cgggggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATc                                                                                                      <  1:1745663/150‑1 (MQ=255)
                                                    gggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGa                                                                                                   >  1:3829237/1‑150 (MQ=255)
                                                         gcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat                                                                                              >  1:4407814/1‑150 (MQ=255)
                                                          cgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                                             >  2:3374338/1‑150 (MQ=255)
                                                              ggcaggcaGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCt                                                                                         >  1:1048820/1‑150 (MQ=255)
                                                                 aggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa                                                                                      <  1:983924/150‑1 (MQ=255)
                                                                   gcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatg                                                                                    >  2:3228286/1‑150 (MQ=255)
                                                                      gAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                                 <  1:1097762/150‑1 (MQ=255)
                                                                      gAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                                 <  1:4298696/150‑1 (MQ=255)
                                                                        accagcaccagcAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGt                                                                               <  1:1841295/150‑1 (MQ=255)
                                                                           agcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGCTGGTGCt                                                                            >  2:3491829/1‑150 (MQ=255)
                                                                           agcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCt                                                                            >  2:76228/1‑150 (MQ=255)
                                                                             caccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCttt                                                                          >  1:1015316/1‑150 (MQ=255)
                                                                                cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCAGCAGGTCATCTGAATGATGGTCTGCCTGAGGGTGCTTTCAg                                                                       >  2:4105546/1‑150 (MQ=255)
                                                                                cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCGTTGCGGCCCCAATTACAAAAACCAAAGCCAATACAATCATCAGGTAAACTGAATGATGGTCTGCATGATGGTGCTTTCAg                                                                       >  1:3403939/1‑150 (MQ=255)
                                                                                 agcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGt                                                                      <  2:435557/150‑1 (MQ=255)
                                                                                  gcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGc                                                                                                                                                <  2:3784066/75‑1 (MQ=255)
                                                                                  gcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTg                                                                     >  1:4518443/1‑150 (MQ=255)
                                                                                  gcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTg                                                                     <  1:3440366/150‑1 (MQ=255)
                                                                                  gcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTg                                                                     <  2:1434486/150‑1 (MQ=255)
                                                                                       gggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAAt                                                                >  2:2181402/1‑150 (MQ=255)
                                                                                        ggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATg                                                               <  2:3563938/150‑1 (MQ=255)
                                                                                           aCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTc                                                            <  1:173746/150‑1 (MQ=255)
                                                                                           aCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTc                                                            <  1:1054483/150‑1 (MQ=255)
                                                                                             aaGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCtt                                                          <  2:4261391/150‑1 (MQ=255)
                                                                                              aGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                         <  1:3725845/150‑1 (MQ=255)
                                                                                                gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                                       >  2:3242372/1‑150 (MQ=255)
                                                                                                gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCGTTGCGGCCCCAATTACAAACACCAAAGCCCATACAAACAACAGGTCATCTGAATGATGGTCTGCAGGATGGGGCTTTCAGGGGAAATGGGCTTGGc                                                       >  1:2286208/1‑150 (MQ=255)
                                                                                                    gTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTg                                                   <  1:3113189/150‑1 (MQ=255)
                                                                                                     tAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTgg                                                  >  2:2897784/1‑150 (MQ=255)
                                                                                                       aCGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGcc                                                <  2:2937790/150‑1 (MQ=255)
                                                                                                        cGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCt                                               <  1:2649638/150‑1 (MQ=255)
                                                                                                        cGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCt                                               <  2:2419096/150‑1 (MQ=255)
                                                                                                         gACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCATTGCGGCCCCAAATACAAACACCAAAGCCAAAACAAACATCAGGTCAACTGAAAGAAGGTCTGCATGATGGTGCTTTCAGTGGAAATGGACGTGGCATTGGccgg                                              >  2:538233/1‑148 (MQ=255)
                                                                                                          aCAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTgg                                             >  2:4170046/1‑150 (MQ=255)
                                                                                                              ccGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGt                                         >  1:2160330/1‑150 (MQ=255)
                                                                                                                gCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTg                                       <  1:1099826/150‑1 (MQ=255)
                                                                                                                 cTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGc                                      >  2:2837190/1‑150 (MQ=255)
                                                                                                                    cATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGctct                                   <  1:3174964/150‑1 (MQ=255)
                                                                                                                        cgcgTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgcgc                               <  2:4383046/150‑1 (MQ=255)
                                                                                                                         gcgTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc                              <  1:1350432/150‑1 (MQ=255)
                                                                                                                         gcgTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc                              >  1:3224295/1‑150 (MQ=255)
                                                                                                                           gTGAACCTTCACATTCTTGGGGCCGCTTTTCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATTGTCTTGGCATTGGCCTGGTAGTTTCTCTGCGCCtt                            >  1:1997959/1‑150 (MQ=255)
                                                                                                                                ccTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCa                       <  1:2985450/150‑1 (MQ=255)
                                                                                                                                 cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAACGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                      <  1:3877700/150‑1 (MQ=255)
                                                                                                                                 cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                      >  2:192392/1‑150 (MQ=255)
                                                                                                                                 cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGGGGAAATGGTCTTGGCCTTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                      >  2:1619949/1‑150 (MQ=255)
                                                                                                                                   tCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGTTGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGGTCTGCGCCTTGATCAGGt                    >  2:3650981/1‑150 (MQ=255)
                                                                                                                                   tCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAACATCAGGTCATCTGAATGATGGTCTGCATGCTGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGt                    >  2:1924881/1‑150 (MQ=255)
                                                                                                                                   tCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAACATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATGAGGt                    >  2:890238/1‑150 (MQ=255)
                                                                                                                                   tCACATTCTTGGGGCCGCTTTGCGGCCCCAAATACAAACACCAAAGCCAATACAAAAATCAGGTCAAATGAATTATTGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTTGCATTGGCCGGGTAGTTTCTCTGCGCCGTGATCAGGt                    >  1:3546055/1‑150 (MQ=255)
                                                                                                                                   acacaATCTTGGGGCCCGTTTGCGGCCCCAAATACAAACACCAAAGCCAATACAAACAACACGTCAAATCTATCCTGGTATTCATGCTGGTTCTTTCAGATGAAAACGTCGTGGCATTGGCCGTGTAGTTGCTCTGCGCCATGATCATGt                    >  2:2637531/2‑150 (MQ=255)
                                                                                                                                       aTTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAc                >  2:1303809/1‑150 (MQ=255)
                                                                                                                                       aTTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAc                <  2:1209892/150‑1 (MQ=255)
                                                                                                                                        ttCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAcc               >  1:1616674/1‑150 (MQ=255)
                                                                                                                                         tCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCa              <  1:1782343/150‑1 (MQ=255)
                                                                                                                                          cTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAg             <  2:3632542/150‑1 (MQ=255)
                                                                                                                                           ttGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGc            >  2:1478267/1‑150 (MQ=255)
                                                                                                                                           ttGGGGCCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGc            >  2:7713/1‑150 (MQ=255)
                                                                                                                                           ttGGGGCCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGc            >  1:1802952/1‑150 (MQ=255)
                                                                                                                                           ttGGGGCCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGc            >  2:1329209/1‑150 (MQ=255)
                                                                                                                                              gggCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCg         >  2:3898541/1‑150 (MQ=255)
                                                                                                                                              gggCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCg         >  2:3968732/1‑150 (MQ=255)
                                                                                                                                                gCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGcc       >  2:2305262/1‑150 (MQ=255)
                                                                                                                                                  cGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCgg     >  1:3790560/1‑150 (MQ=255)
                                                                                                                                                    cTTTGCGGCCCCAATTACAAACACCAAACCCAATACAAACAACAGCACATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGCTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGcccctc   <  2:2017146/150‑5 (MQ=255)
                                                                                                                                                     tttGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGGCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTCGTTGCTCTGCGCCTTGATCAGGTTTACCCGCTCGCCGGTcc  >  2:3434356/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCA  >  NC_002947/4982410‑4982708

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: