Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,110 T→A 8.1% intergenic (+71/+81) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,1100TA 8.1% 515.8 / 4.7 180intergenic (+71/+81)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (88/68);  new base A (17/2);  total (110/70)
Fisher's exact test for biased strand distribution p-value = 5.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

GCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGG  >  NC_002947/667961‑668258
                                                                                                                                                     |                                                                                                                                                    
     tCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATATGGCCGCAAAGCGGCCCCGACAAGTCAAg                                                                                                                                                 >  2:576410/1‑150 (MQ=255)
     tCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCCTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAg                                                                                                                                                 >  2:6135468/1‑150 (MQ=255)
      cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGTTGGGCAAGGGGAGTTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgc                                                                                                                                                <  1:9016602/150‑1 (MQ=255)
      cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgc                                                                                                                                                <  2:4308065/150‑1 (MQ=255)
         aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcggc                                                                                                                                             <  2:1772045/150‑1 (MQ=255)
         aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGAAAGGGGAGATGACGGTGAAAGCCCGCTAGACACCCCCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcggc                                                                                                                                             <  2:2398949/150‑1 (MQ=38)
         aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGGGCGGTGGGGCGGGCAAGCGGCCCCGCCCAGTCAAtcggg                                                                                                                                             >  1:6936301/1‑149 (MQ=255)
            gATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                                                          <  2:5287628/150‑1 (MQ=255)
            gATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGaa                                                                                                                                          >  2:699450/1‑150 (MQ=255)
             aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCCTTAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAg                                                                                                                                         >  2:1179929/1‑150 (MQ=255)
                 aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGa                                                                                                                                     <  1:5686648/150‑1 (MQ=255)
                  ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCGGGGCGGGGCGGGTGGGGCGCCAAGCGGCCCCGACAAGTCAAGCGGCGAAGCGGac                                                                                                                                    >  1:1601056/1‑149 (MQ=255)
                   cGGCCTAACCTGCAGGCGGTGGGCGATGGGAATGGGTATATGACGGTTAAAGCCCGCTTGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCTATCGGGCCGAAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                                   <  1:5397603/150‑1 (MQ=37)
                   cGGCCGAACCTGCTGGCGGTGGGCGTGGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGGGTAATGTGTGGGCGCGGCCTTGCGTCGCGATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                                   <  1:4833021/150‑1 (MQ=37)
                   cGGCCGAACCTGCAGGCGGTGGGCGATGGGAATGGGGTGATTACGGTTAAAGCCCGCTCGAGACCCGCTCTACGCGGAAAGTGTGGGAGGGGCCTTGGGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                                   <  1:3151362/150‑1 (MQ=255)
                   cGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGCGCGCTAGGGCCGCCAAACCGCCCCGCCCATTCAAGCGGCGAAGcacacc                                                                                                                                   >  1:5123630/1‑145 (MQ=255)
                   cGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCCAAGCGGCCCCCACAATTCAAGCGGCGACGCTGgaa                                                                                                                                   >  2:6912481/1‑150 (MQ=255)
                        gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTg                                                                                                                              >  2:6539462/1‑150 (MQ=255)
                         aaCCTGCTGGCCGTGGGCGATGGGCACGGGGAGTTGACGGTGAACCCCCGCTCGCGACCCGCTCTACGCGGAATGTGTGGGAGGGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                                             <  1:3301644/150‑1 (MQ=255)
                         aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                                             <  2:3246878/150‑1 (MQ=255)
                              gctggCCGTGGGCTTTGGGCATGGGGATATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGGGTCGCGATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctg                                                                                                                        <  1:4179694/147‑1 (MQ=38)
                                aGGCCGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                                      <  2:306475/150‑1 (MQ=255)
                                  gcggggggCGTTGGGCATGGGTAGTTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAATGTGTGGGAGCGGCCTTGCGTCGCTATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgt                                                                                                                    <  1:1104616/145‑1 (MQ=37)
                                  gcggggggCGTTGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAACGGGTGGGCGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgt                                                                                                                    <  1:1104555/145‑1 (MQ=255)
                                   ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCTATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                                                   <  1:6482235/150‑1 (MQ=255)
                                   ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                                                   <  2:7129599/150‑1 (MQ=255)
                                   ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                                                   <  2:1369845/150‑1 (MQ=255)
                                    cGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGc                                                                                                                  >  2:6806727/1‑150 (MQ=255)
                                    cGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGc                                                                                                                  >  2:5871212/1‑150 (MQ=255)
                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                                 <  1:4511945/150‑1 (MQ=255)
                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                                 <  2:2128318/150‑1 (MQ=255)
                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                                 <  2:158963/150‑1 (MQ=255)
                                       gggCGTGGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGCGCGGCCTTGCGAGGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                               <  1:893730/150‑1 (MQ=38)
                                       gggCGAGGGGCATGGGGAGATGAGGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                               <  1:7884099/150‑1 (MQ=255)
                                        ggCGATGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                                              <  2:834539/150‑1 (MQ=255)
                                          cGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCTATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                                            <  1:1917380/150‑1 (MQ=255)
                                           gATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCAGGGGCTGCTGTGCAGCCCTg                                                                                                           >  2:5232135/1‑150 (MQ=255)
                                            ttgGGCATGGGGTTATGACGGTGAAAGCCCGCTAGAGACCCGCTCTCCGCGGAACGTGTGGGAGCGGCGTTGCGTGGCTATCGGGCCCCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGt                                                                                                          >  1:3871822‑M2/2‑149 (MQ=255)
                                             tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGACTGTGTGGGAGCGGCCTTGCGTCGCGATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                                                         >  2:2642898‑M2/1‑148 (MQ=255)
                                              gggCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcg                                                                                                        >  2:7304277‑M2/1‑147 (MQ=255)
                                                 cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                                                     >  1:4602684‑M2/1‑144 (MQ=255)
                                                      ggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                                                                >  1:6804730‑M2/1‑139 (MQ=255)
                                                      ggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                                                                <  2:614634‑M2/150‑12 (MQ=255)
                                                        agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                                              >  1:5540173‑M2/1‑137 (MQ=255)
                                                        agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                                              >  2:8678570‑M2/1‑137 (MQ=255)
                                                             aCGGTGAACGCCCGCTAGAGACCCGATCTACGCGGAATGTGTGGGCGCGGCCTTGCGTCGCTATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgct                                                                                         >  1:2137990‑M2/1‑132 (MQ=255)
                                                               ggTGCAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCCCAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcc                                                                                       <  2:2896398‑M2/150‑21 (MQ=255)
                                                               ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccactcc                                                                                       <  2:4921382‑M2/150‑21 (MQ=255)
                                                               ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCGCCCAGCGGCCCCGACAATTGAAGCGGCGGAGCTGGCAGCCGGGGGCTGCTGTGCCGCCCTGTCGCGGCGCCCGGCCGCGcc                                                                                       >  2:1884649/1‑150 (MQ=255)
                                                                      gCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCaa                                                                                                                                                  <  2:922117/84‑1 (MQ=255)
                                                                      gCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCaa                                                                                                                                                  >  1:922117/1‑84 (MQ=255)
                                                                        ccGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCCCAAAGCGGCCCCGACAAATAAAGCGGCGCAGCGGAAAGCCTGGGGCGGCTGTGCAGCCCAGTCGCGCCGCCACGCCGCGCCCCCCAAGAc                                                                              >  2:8410683/1‑150 (MQ=255)
                                                                           cTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGGGTGGCGATTGGGCCGCGAAGCGGGCCCGACCAATGAAGGGGGGGAGGTGTAAGGCTGGGGGTGGTGTTGAGGCCTGTGGGGTCGGAAGGGCGGTGCCCCAAAGaaacc                                                                           >  2:4031209/1‑146 (MQ=255)
                                                                                 ccccGCTCTCCGGGTAATGTGTGGGTGGGGCGTGGGGTCGCTATCGGCCCGCACCGCGGCCCCGACAATACACGCGGCGACGCTCAATGCCGGGGGCTGCTTTGCCGCCCTGGGGGGACGAACGGCCGCTCCCACAAAGACCGCGCAGAt                                                                     <  2:1236501/149‑1 (MQ=255)
                                                                                   ccGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGCCAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCATGGCCGCTCCCACAAAGACCGCGCAGATCa                                                                   <  2:7477392/150‑1 (MQ=37)
                                                                                   ccGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCCCCAATCGGCCCCCAACAATAAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCCCTACCACAAAGACCGCGCAGAtaa                                                                   >  2:2930525/1‑148 (MQ=25)
                                                                                       tctACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                                            <  2:497841/105‑1 (MQ=255)
                                                                                              ggAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCCCGCAGACCAGTGGATCAccc                                                        >  2:9018184/1‑150 (MQ=255)
                                                                                                    tgtgGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCCGCCCCGCCACTTCAAGCGCCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCGCGCCCCCAAAGAACGCGCAGCTCCGTGGCTCACCCACCtca                                                  >  2:8457428/1‑150 (MQ=21)
                                                                                                            cGGCCTTGCGTCGCGTTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCCAAGCTGAAATCCTGGGGCTGCTGTGAAGCCCTGTCGCGACGCAATGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCcac                                          >  1:8283449/1‑150 (MQ=37)
                                                                                                             ggCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCcaca                                         >  1:8537128/1‑150 (MQ=255)
                                                                                                                 ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc                                     >  2:2038031/1‑150 (MQ=255)
                                                                                                                 ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc                                     >  2:6848268/1‑150 (MQ=255)
                                                                                                                    gggggcgATTGGGCCGCAACGGGGCCCCGAAAATTCCCGCGGCGACGGTTAACGGCTGGGGCGGCTGTGCCGCCCGGGGGCGCCGCATGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAg                                  <  2:8802233/146‑1 (MQ=255)
                                                                                                                      tagcgATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGACTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGAGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCa                                >  1:8415274/3‑150 (MQ=255)
                                                                                                                      tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttg                                <  1:2398213‑M2/150‑76 (MQ=255)
                                                                                                                      tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtct                                                              >  2:4074397‑M2/1‑75 (MQ=255)
                                                                                                                         cgATGGGGCCGCAACGGGGCCCCGACAATACAAGCGGCTAAGCTTACAGCGTGGGGCGGGTGTGCAGCCCGGTCGCGACGAACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAcca                             <  2:7783363/150‑1 (MQ=255)
                                                                                                                          gATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgccca                            <  1:5781005‑M2/150‑80 (MQ=255)
                                                                                                                           aTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagc                           <  2:7697732/150‑1 (MQ=255)
                                                                                                                             aGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCTAGGCGGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCCGCTCATCACCCACAGCCCCAGCAccagccc                         >  2:2939635/1‑149 (MQ=255)
                                                                                                                              gggCCGCTTTGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagccag                        >  1:4582209/1‑150 (MQ=255)
                                                                                                                               ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga                       >  1:2967482/1‑150 (MQ=255)
                                                                                                                               ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga                       >  2:8525673/1‑150 (MQ=255)
                                                                                                                                gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgat                      >  1:6972531/1‑150 (MQ=255)
                                                                                                                                gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgat                      >  2:339924/1‑150 (MQ=255)
                                                                                                                                 ccGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatg                     >  2:8851337/1‑150 (MQ=255)
                                                                                                                                  cGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatga                    >  1:6981529/1‑150 (MQ=255)
                                                                                                                                  cGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatga                    >  1:6981423/1‑150 (MQ=255)
                                                                                                                                  cGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatga                    >  1:7480891/1‑150 (MQ=255)
                                                                                                                                    cAAAGCGGCCCCTACAATCCAAGCGGCTATGCTGAACGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGAtt                  <  2:1010554/150‑1 (MQ=255)
                                                                                                                                    cAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGAtt                  >  2:5351194/1‑150 (MQ=255)
                                                                                                                                     aaaGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTc                 >  1:7535978/1‑150 (MQ=255)
                                                                                                                                     aaaGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTc                 >  2:8104421/1‑150 (MQ=255)
                                                                                                                                        gCGGCCCCGACAATTCCAGCGGCGACGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGCCGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGa              <  2:5214671/150‑1 (MQ=255)
                                                                                                                                         cgcccccCACAATTCACGGGGCGACGCTGAATGCGTGGGGGTGGTGTGCCGCCCTGCGGCGCCGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAg             <  2:4801621/147‑1 (MQ=255)
                                                                                                                                           gCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTTTGCGGCCCTAtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccatcgcccacggcctgcag           <  1:3763588‑M2/150‑97 (MQ=255)
                                                                                                                                            ccccGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgat          >  1:4197360/1‑150 (MQ=255)
                                                                                                                                             cccTACAATTCAAGCCGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGACCAGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGAGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgatg         >  2:8434608/1‑150 (MQ=39)
                                                                                                                                             cccGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgatg         >  2:8701958/1‑150 (MQ=255)
                                                                                                                                             accgACAATTCAATCGGCTAAGCAGAAAGCATGGGGCTTCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacaaagaccgcgcagagagggtctatagcgggctttcaccgacatctccccttgcccatcgcccacggcctgcaggt         <  2:5099012‑M2/149‑99 (MQ=255)
                                                                                                                                              ccGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgatga        >  1:7109293/1‑150 (MQ=255)
                                                                                                                                              ccGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgatga        >  2:5292402/1‑150 (MQ=255)
                                                                                                                                              ccGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgatga        >  2:8020039/1‑150 (MQ=255)
                                                                                                                                                gACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATc      >  1:2263741/1‑150 (MQ=255)
                                                                                                                                                gACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATc      >  2:3535232/1‑150 (MQ=255)
                                                                                                                                                gACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATc      >  2:7886553/1‑150 (MQ=255)
                                                                                                                                                gACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATc      >  1:7811833/1‑150 (MQ=255)
                                                                                                                                                   aaTTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAg   >  2:2179723/1‑150 (MQ=255)
                                                                                                                                                   aaTTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCACCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGTTGATTCCGAGGATGATagag   >  2:7611528/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGG  >  NC_002947/667961‑668258

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