Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,069 T→C 5.1% intergenic (‑110/+18) PP_0386 ← / ← rpoD sensory box protein/RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,0690TC 5.1% 1224.0 / 31.5 372intergenic (‑110/+18)PP_0386/rpoDsensory box protein/RNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base T (174/177);  new base C (9/10);  total (185/187)
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.07e-01

Only 100 of 378 total aligned reads displayed.

AAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTG  >  NC_002947/470920‑471218
                                                                                                                                                     |                                                                                                                                                     
  aTCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTgg                                                                                                                                                     <  2:6939007/150‑1 (MQ=255)
        aTCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCATGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAATACCCCGGCCTGGGCCggg                                                                                                                                               <  1:6729383/150‑1 (MQ=255)
         tCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCt                                                                                                                                                       >  1:1971063/1‑141 (MQ=255)
           aGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtt                                                                                                                                            >  2:882296/1‑150 (MQ=255)
            gCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGttt                                                                                                                                           <  1:2190778/150‑1 (MQ=255)
            gCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGttt                                                                                                                                           >  2:3256132/1‑150 (MQ=255)
            gCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGttt                                                                                                                                           >  2:1748738/1‑150 (MQ=255)
            gCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGCCGGGGttt                                                                                                                                           <  1:577382/150‑1 (MQ=255)
             cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtttt                                                                                                                                          <  1:6787990/150‑1 (MQ=255)
               aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGt                                                                                                                                        <  2:2887518/150‑1 (MQ=255)
               aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGt                                                                                                                                        >  1:4976712/1‑150 (MQ=255)
                cAGAAGGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtc                                                                                                                                       <  1:6233183/150‑1 (MQ=255)
                 agcaTGGCCGGCTTCAAGGACAGGGCCGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtca                                                                                                                                      <  1:294175/147‑1 (MQ=255)
                  gAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtcat                                                                                                                                     <  1:3349117/150‑1 (MQ=255)
                   aaTGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGCCGGGGtttttcttta                                                                                                                                    <  2:3321711/150‑7 (MQ=255)
                    aTGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtcatca                                                                                                                                   <  1:6140103/150‑1 (MQ=255)
                     tGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCAc                                                                                                                                  >  2:1266032/1‑150 (MQ=255)
                       gCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACtc                                                                                                                                >  1:2659418/1‑150 (MQ=255)
                          ggCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACtcgtc                                                                                                                             >  2:2325228/1‑150 (MQ=255)
                          ggCTTCAAGGACATGGCAGCCTCGCACCCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACtcgtc                                                                                                                             >  2:41849/1‑150 (MQ=255)
                             ttCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAATGAAATACCCCGGCCTGGGCCGGGGTGGTGTCATCACTCGTCGag                                                                                                                          >  2:1111764/1‑150 (MQ=255)
                             ttCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGag                                                                                                                          >  2:4999575/1‑150 (MQ=255)
                             ttCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGag                                                                                                                          <  2:1930707/150‑1 (MQ=255)
                             ttCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGag                                                                                                                          >  2:7209624/1‑150 (MQ=255)
                              tCAAGGACATGGCAGCCTCGCACTCAGGTGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAATTCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGagga                                                                                                                        >  1:282906/1‑150 (MQ=255)
                                aaGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggac                                                                                                                       >  1:4537180/1‑149 (MQ=255)
                                aaGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggac                                                                                                                       >  1:5440058/1‑149 (MQ=255)
                                aaGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaa                                                                                                                       >  2:3514938/1‑150 (MQ=255)
                                aaGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCGGGGCCGGGGTGTTGTCATCACTCGTCGaggac                                                                                                                       >  2:3822256/1‑149 (MQ=255)
                                 aGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaag                                                                                                                      <  1:3037256/150‑1 (MQ=255)
                                 aGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaag                                                                                                                      <  1:1023929/150‑1 (MQ=255)
                                        ggCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCa                                                                                                               >  1:799363/1‑150 (MQ=255)
                                           aGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGt                                                                                                            >  1:2356585/1‑150 (MQ=255)
                                            gCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTg                                                                                                           <  2:3877246/150‑1 (MQ=255)
                                               tCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTc                                                                                                        <  2:5464691/150‑1 (MQ=255)
                                                    cTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                                   <  2:5569501/150‑1 (MQ=255)
                                                    cTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                                   <  2:6356817/150‑1 (MQ=255)
                                                     tCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGTGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCg                                                                                                  >  2:6840065/1‑150 (MQ=255)
                                                      cAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGa                                                                                                 <  2:7129877/150‑1 (MQ=255)
                                                       aGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAg                                                                                                >  2:1376619/1‑150 (MQ=255)
                                                               tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTgg                                                                                        <  1:5027394/150‑1 (MQ=255)
                                                               tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTgg                                                                                        >  1:3605028/1‑150 (MQ=255)
                                                               tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTgg                                                                                        >  2:5629197/1‑150 (MQ=255)
                                                                  aGTCGAGTGTAGACGGGCATTGCCGCAGGAGAGCAACGGAAGGCGGGCTTAATGTCACGGGCCACTAAACAAAAACCCCGGCCTGGGCCGGGGAATTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGgcgc                                                                                     <  2:6993308/150‑1 (MQ=255)
                                                                     cGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGGGTCGGGTGGCGCAGc                                                                                  >  1:966884/1‑150 (MQ=255)
                                                                        gtgtAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTg                                                                               <  2:638002/150‑1 (MQ=255)
                                                                        gtgtAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTg                                                                               >  1:7234898/1‑150 (MQ=255)
                                                                          gtAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgcg                                                                             <  2:789728/150‑1 (MQ=255)
                                                                          gtAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgcg                                                                             >  2:4519426/1‑150 (MQ=255)
                                                                              aCGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtcat                                                                                                                                     <  1:350544/90‑1 (MQ=255)
                                                                              aCGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAAc                                                                         <  1:1229023/150‑1 (MQ=255)
                                                                                gggCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGc                                                                       >  1:4926827/1‑150 (MQ=255)
                                                                                    aTTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttc                                                                   >  1:3798801/1‑150 (MQ=255)
                                                                                     ttGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcg                                                                  <  1:3142485/150‑1 (MQ=255)
                                                                                      tGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgc                                                                 <  1:5214082/150‑1 (MQ=255)
                                                                                        ccgccgGTGTGCACCGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgctt                                                               <  2:6319943/150‑1 (MQ=255)
                                                                                          gCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgcttcg                                                             <  2:4070658/150‑1 (MQ=255)
                                                                                                gtgCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGa                                                       >  2:5989263/1‑150 (MQ=255)
                                                                                                   cAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACgg                                                    >  2:4911794/1‑150 (MQ=255)
                                                                                                   cAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACgg                                                    >  2:2250328/1‑150 (MQ=255)
                                                                                                             cGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAc                                          >  1:5855118/1‑150 (MQ=255)
                                                                                                               ggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCg                                                                                             >  1:3934157/1‑97 (MQ=255)
                                                                                                               ggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgc                                        <  2:6692573/150‑1 (MQ=255)
                                                                                                               ggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgc                                        <  2:7346206/150‑1 (MQ=255)
                                                                                                               ggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgc                                        <  1:4226968/150‑1 (MQ=255)
                                                                                                                 cTTTATGTCACGCGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGt                                      <  1:4120492/150‑1 (MQ=255)
                                                                                                                  tttATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGATTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGtt                                     <  1:430687/150‑1 (MQ=255)
                                                                                                                    tATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTAc                                   <  2:5291979/150‑1 (MQ=255)
                                                                                                                        tCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCa                               >  1:1613148/1‑150 (MQ=255)
                                                                                                                         cACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaa                              >  1:1293859/1‑150 (MQ=255)
                                                                                                                          aCGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaaa                             <  1:7138205/150‑1 (MQ=255)
                                                                                                                            gggCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACt                           >  2:4305037/1‑150 (MQ=255)
                                                                                                                              gCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGc                         <  2:2100650/150‑1 (MQ=255)
                                                                                                                              gCCACTAAAGAAAAACCCCGGCCTGGGCCGAGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGc                         <  1:6354857/150‑1 (MQ=255)
                                                                                                                               ccACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGAt                                                           <  2:593989/115‑1 (MQ=255)
                                                                                                                               ccACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCt                        <  1:882296/150‑1 (MQ=255)
                                                                                                                                cACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCtt                       >  1:1352095/1‑150 (MQ=255)
                                                                                                                                 aCTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTg                      >  2:1501309/1‑150 (MQ=255)
                                                                                                                                   tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc                    >  2:4579552/1‑150 (MQ=255)
                                                                                                                                   tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc                    <  2:5198344/150‑1 (MQ=255)
                                                                                                                                   tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc                    >  1:4710106/1‑150 (MQ=255)
                                                                                                                                   tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc                    >  2:1023847/1‑150 (MQ=255)
                                                                                                                                     aaGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCa                  >  1:2258567/1‑150 (MQ=255)
                                                                                                                                     aaGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCa                  >  1:6897083/1‑150 (MQ=255)
                                                                                                                                      aGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAc                 >  1:1436763/1‑150 (MQ=255)
                                                                                                                                      aGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAc                 <  1:5940294/150‑1 (MQ=255)
                                                                                                                                       gAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgcg                                                                             >  2:5830492/1‑89 (MQ=255)
                                                                                                                                       gAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAcc                >  2:36763/1‑150 (MQ=255)
                                                                                                                                       gAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAcc                >  2:2677384/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAcc                <  2:3994975/148‑1 (MQ=255)
                                                                                                                                            aCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTc           >  1:2511229/1‑150 (MQ=255)
                                                                                                                                             ccccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCg          >  1:737086/1‑150 (MQ=255)
                                                                                                                                              cccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa         >  2:4443872/1‑150 (MQ=255)
                                                                                                                                              cccGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa         >  2:3088364/1‑150 (MQ=255)
                                                                                                                                              cccGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa         >  2:6601922/1‑150 (MQ=255)
                                                                                                                                               ccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAg        >  1:555290/1‑150 (MQ=255)
                                                                                                                                                cGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg       <  1:843229/150‑1 (MQ=255)
                                                                                                                                                cGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg       <  1:1266032/150‑1 (MQ=255)
                                                                                                                                                   ccTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGgtg    >  1:6559380/1‑150 (MQ=255)
                                                                                                                                                    cTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGgtgt   >  1:5641573/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTG  >  NC_002947/470920‑471218

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