Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 646,275 T→A 10.7% intergenic (‑178/‑106) PP_0556 ← / → acoR acetoin catabolism protein/acetoin catabolism regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947646,2750TA10.7% 224.3 / 16.2 85intergenic (‑178/‑106)PP_0556/acoRacetoin catabolism protein/acetoin catabolism regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (34/41);  new base A (6/3);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

GCAGGAGCAAGCGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCATGTGGACTGCGTCAGC  >  NC_002947/646126‑646422
                                                                                                                                                     |                                                                                                                                                   
gCAGGAGCAAGCGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCt                                                                                                                                                     >  1:2178322/1‑150 (MQ=255)
gCAGGAGCAAGCGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCt                                                                                                                                                     <  2:112509/150‑1 (MQ=255)
    gAGCAAGCGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCGGCGt                                                                                                                                                 >  1:357538/1‑150 (MQ=255)
    gAGCAAGCGGGGTGCCAGCTTTCGTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGt                                                                                                                                                 <  2:157746/150‑1 (MQ=255)
     aGCAAGCGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGTCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGtc                                                                                                                                                <  2:832059/150‑1 (MQ=255)
         aGCGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCaa                                                                                                                                            <  1:2008722/150‑1 (MQ=255)
           cGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAagc                                                                                                                                          <  1:2282808/150‑1 (MQ=255)
                  ccAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCGCAAGCCAGCCgg                                                                                                                                   >  1:1932579/1‑150 (MQ=255)
                   cAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGt                                                                                                                                  >  1:336936/1‑150 (MQ=255)
                    aGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAag                                                                                                                                           >  1:1297855/1‑140 (MQ=255)
                    aGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAag                                                                                                                                           <  2:1297855/140‑1 (MQ=255)
                    aGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTa                                                                                                                                 >  2:1055112/1‑150 (MQ=255)
                    aGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTa                                                                                                                                 >  1:1541649/1‑150 (MQ=255)
                                gATGGCTGTTGGGGCCGGGGGAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGc                                                                                                                     <  2:61255/150‑1 (MQ=255)
                                    gCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGATAATCACACAGTGCTTGt                                                                                                                 <  2:2277227/150‑1 (MQ=255)
                                       gAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCgg                                                                                                              <  2:2370976/150‑1 (MQ=255)
                                             gCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATg                                                                                                        >  1:2272671/1‑150 (MQ=255)
                                              ccGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGc                                                                                                       <  2:2272671/150‑1 (MQ=255)
                                              ccGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGc                                                                                                       <  1:1493662/150‑1 (MQ=255)
                                               cGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGAGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATccc                                                                                                      >  1:492737/1‑150 (MQ=255)
                                               cGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTCCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGcc                                                                                                      <  1:1061078/150‑1 (MQ=255)
                                                gggCTAGAGGGCTCTGATCGTAATCGGGCAGGCAATCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCTCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGccc                                                                                                     >  2:869726/1‑150 (MQ=255)
                                                  gCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCag                                                                                                   >  2:1343577/1‑150 (MQ=255)
                                                   gtAGAGGGCTCGGGTCGGGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCaga                                                                                                  <  2:1063139/149‑1 (MQ=255)
                                                   cTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCaga                                                                                                  <  1:742176/150‑1 (MQ=255)
                                                    tAGCGGGCTCTGGTCGTGTTTGGGCAGGTTTTCACGGGTTGGGTGGCGAGACGCCGCGTGCCAGCGCCGAGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCCAGCCAGTCGGGAATCACACAGTGCTTGTCGGGGCATGCCTaggg                                                                                                 <  2:395655/150‑3 (MQ=255)
                                                    tAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCagag                                                                                                 <  2:336936/150‑1 (MQ=255)
                                                    tAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCAGCGTCTCag                                                                                                                                            <  2:1192298/107‑2 (MQ=255)
                                                    tAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCAGCGTCTCag                                                                                                                                            >  1:1192298/1‑106 (MQ=255)
                                                       agGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGgcg                                                                                              >  2:1347200/1‑150 (MQ=255)
                                                          gCTCTGGTCGTGTTCGGGCAGGTTTTCACTGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCtt                                                                                           >  2:1786664/1‑150 (MQ=255)
                                                           ctctGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTa                                                                                          <  1:1347200/150‑1 (MQ=255)
                                                           ctctGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGACCAGAGGCGCTTa                                                                                          >  2:13328/1‑150 (MQ=255)
                                                            tctGGTCGTGTTCGGGCAGGTTTTCACTGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTaa                                                                                         <  1:1786664/150‑1 (MQ=255)
                                                            tctGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTaa                                                                                         >  1:1491344/1‑150 (MQ=255)
                                                            tatGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTaa                                                                                         <  1:2192351/148‑1 (MQ=255)
                                                              tGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATg                                                                                       <  1:1490315/150‑1 (MQ=255)
                                                                 tCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGGGGCGCTTAATGGcc                                                                                    >  1:620256/1‑150 (MQ=255)
                                                                 tCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGcc                                                                                    <  1:13328/150‑1 (MQ=255)
                                                                 tCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTCCCTCTCAAGCCACCCGGTCATCACCCCGTGCTTGTTTGGGCATGCCCAGAAGCGCTCAATGGcc                                                                                    <  2:394565/150‑1 (MQ=255)
                                                                  cGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATAACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCg                                                                                   <  2:462054/150‑1 (MQ=255)
                                                                    tgtTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGa                                                                                 <  1:1343577/150‑1 (MQ=255)
                                                                    tgtTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGa                                                                                 >  2:1766220/1‑150 (MQ=255)
                                                                        cGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAAc                                                                             <  2:1491942/150‑1 (MQ=255)
                                                                               gTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATaaaa                                                                      >  1:1324184/1‑150 (MQ=255)
                                                                               gTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATaaaa                                                                      >  2:2261303/1‑150 (MQ=255)
                                                                               gTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATaaaa                                                                      >  1:135866/1‑150 (MQ=255)
                                                                                ttttCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCAt                                                                                                         >  1:1389624/1‑114 (MQ=255)
                                                                                ttttCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCAt                                                                                                         <  2:1389624/114‑1 (MQ=255)
                                                                                         gTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACga                                                            <  2:1324184/150‑1 (MQ=255)
                                                                                         gTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACga                                                            <  2:135866/150‑1 (MQ=255)
                                                                                          tAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACgag                                                           >  1:89588/1‑150 (MQ=255)
                                                                                             ggTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGgc                                                                                               <  2:1613642/111‑1 (MQ=255)
                                                                                             ggTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGgc                                                                                               >  1:1613642/1‑111 (MQ=255)
                                                                                              gTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAAcc                                                       >  1:631904/1‑150 (MQ=255)
                                                                                              gTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAAcc                                                       >  1:631893/1‑150 (MQ=255)
                                                                                                  cGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAg                                                   <  1:1055112/150‑1 (MQ=255)
                                                                                                   gagaCGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCAt                                                                                                         >  2:1310556/1‑95 (MQ=255)
                                                                                                   gagaCGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCAt                                                                                                         <  1:1310556/95‑1 (MQ=255)
                                                                                                   gagaCGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATaaaa                                                                      <  1:235098/130‑1 (MQ=255)
                                                                                                   gagaCGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATaaaa                                                                      >  2:235098/1‑130 (MQ=255)
                                                                                                   gagaCGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAgg                                                  <  1:2130846/150‑1 (MQ=255)
                                                                                                        gCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTg                                             >  2:318577/1‑150 (MQ=255)
                                                                                                           gggTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAAt                                          <  2:2178322/150‑1 (MQ=255)
                                                                                                           gggTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAAt                                          <  1:2261303/150‑1 (MQ=255)
                                                                                                                ccAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTg                                     >  1:1200856/1‑150 (MQ=255)
                                                                                                                 cAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGc                                    >  2:1822334/1‑150 (MQ=255)
                                                                                                                   gcgcCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGaa                                                         >  1:749827/1‑127 (MQ=255)
                                                                                                                   gcgcCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGaa                                                         <  2:749827/127‑1 (MQ=255)
                                                                                                                            cTGCGTTGCAGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTccc                         <  2:89588/150‑1 (MQ=255)
                                                                                                                                gTTGCCGAGACACTGGGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGc                     <  1:310028/150‑1 (MQ=255)
                                                                                                                                      gagaCACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCAtgt               <  1:855071/150‑1 (MQ=255)
                                                                                                                                      gagaCACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCAtgt               >  2:330434/1‑150 (MQ=255)
                                                                                                                                         acacTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTCACGAGAACCGGAGGGCATGAATGGTTGCCGCGAACTCCCGAGCCCATGTGGa            >  1:1033698/1‑150 (MQ=255)
                                                                                                                                           cttGAGACGCAGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAAc                                                              <  1:1159748/96‑1 (MQ=255)
                                                                                                                                           cttGAGACGCAGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAAc                                                              >  2:1159748/3‑98 (MQ=255)
                                                                                                                                           acTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCATGTGGACt          >  2:2220735/1‑150 (MQ=255)
                                                                                                                                            ttaagagGCAGCGTGTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCATGTGGACTg         >  2:2312194/4‑150 (MQ=255)
                                                                                                                                            ttGAGACGCAGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGGGAACCGGAGGGCCTGAATGCTTGcc                                   >  1:930170/2‑124 (MQ=255)
                                                                                                                                            ttGAGACGCAGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGGGAACCGGAGGGCCTGAATGCTTGcc                                   <  2:930170/123‑1 (MQ=255)
                                                                                                                                            ttGAGACGCAGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAGCTCCCGAGCCCATGTGGACTg         >  1:1132778/2‑150 (MQ=255)
                                                                                                                                            ttGAGACGCAGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCATGTGGACTg         >  1:1500330/2‑150 (MQ=255)
                                                                                                                                            cTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCATGTGGACTg         <  1:2220735/150‑1 (MQ=255)
                                                                                                                                             tGAGACGCTGCGTCTCAAGGCAGCCGGTAATCACACAGTGCTTGGCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTGCCGAGCCCATGTGGACTGc        <  2:631009/150‑1 (MQ=255)
                                                                                                                                                 aCGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAAc                                                              >  1:1718503/1‑92 (MQ=255)
                                                                                                                                                 aCGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAAc                                                              <  2:1718503/92‑1 (MQ=255)
                                                                                                                                                   gCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCATGTGGACTGCGTcagc  >  2:420907/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCAGGAGCAAGCGGGGTGCCAGCTTTCTTCGCGATGGCTGAAGTGGCCGGGCTAGAGGGCTCTGGTCGTGTTCGGGCAGGTTTTCACGGGTAGGGTGGCGAGACGCCGGGTGCCAGCGCCGGGCCTGCGTTGCCGAGACACTGAGACGCTGCGTCTCAAGCCAGCCGGTAATCACACAGTGCTTGTCGGGGCATGCCCAGAGGCGCTTAATGGCCGGACAACGATAAAAACAGTAACGAGAACCGGAGGGCCTGAATGCTTGCCGCGAACTCCCGAGCCCATGTGGACTGCGTCAGC  >  NC_002947/646126‑646422

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