Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,673 G→C 13.2% intergenic (‑39/‑379) flhA ← / → PP_4345 flagellar export pore protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,936,6730GC13.2% 233.8 / 23.3 83intergenic (‑39/‑379)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base C (8/3);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTAC  >  NC_002947/4936527‑4936821
                                                                                                                                                  |                                                                                                                                                    
catcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGgcgc                                                                                                                                                   >  2:1383217/1‑150 (MQ=35)
catcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGgcgc                                                                                                                                                   <  2:482228/150‑1 (MQ=35)
catcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGgcgc                                                                                                                                                   >  1:2001298/1‑150 (MQ=35)
       atTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCg                                                                                                                                            >  1:2188188/1‑150 (MQ=38)
        ttGCCAACATCACCAGCATCAACAGCGGCACCCCCAGTTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTCCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGt                                                                                                                                           <  1:1127592/150‑1 (MQ=255)
          gCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTgg                                                                                                                                         <  1:1039374/150‑1 (MQ=40)
             aaCATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCtt                                                                                                                                      <  1:570168/150‑1 (MQ=255)
                    ccAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaaga                                                                                                                               >  1:1965379/1‑150 (MQ=255)
                     cagcagCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaagaa                                                                                                                              <  2:1409890/150‑1 (MQ=255)
                     cagcagCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaagaa                                                                                                                              >  1:980554/1‑150 (MQ=255)
                     cagcagCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaagaa                                                                                                                              <  1:58273/150‑1 (MQ=255)
                       ggaggaaCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGc                                                                                                                            >  2:1201687/6‑150 (MQ=255)
                          gcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCttt                                                                                                                         <  2:377731/150‑1 (MQ=255)
                          gcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGt                                                                                                                                           <  2:198085/132‑1 (MQ=255)
                          gcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGt                                                                                                                                           >  1:198085/1‑132 (MQ=255)
                             aCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCa                                                                                                                      >  2:1392796/1‑150 (MQ=255)
                              cAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCaa                                                                                                                     <  1:1121442/150‑1 (MQ=255)
                                   gCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTg                                                                                                                >  1:2182751/1‑150 (MQ=255)
                                     aCCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGc                                                                                                              <  2:980676/150‑1 (MQ=255)
                                     aCCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGc                                                                                                              <  2:1433128/150‑1 (MQ=255)
                                       ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAg                                                                                                            >  1:1711968/1‑150 (MQ=255)
                                       ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAg                                                                                                            >  2:1661756/1‑150 (MQ=255)
                                         ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGac                                                                                                          >  1:1570579/1‑150 (MQ=255)
                                          cAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGaca                                                                                                         <  1:1687056/150‑1 (MQ=255)
                                            ggTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACAcc                                                                                                       <  2:980554/150‑1 (MQ=255)
                                                  cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACg                                                                                                 >  1:1322724/1‑150 (MQ=255)
                                                  cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACg                                                                                                 >  2:464649/1‑150 (MQ=255)
                                                  cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGc                                                                                                                                             >  1:1157404/1‑106 (MQ=255)
                                                  cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGc                                                                                                                                             <  2:1157404/106‑1 (MQ=255)
                                                   ctggCCGGGCTCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGt                                                                                                                                           <  1:2241022/105‑1 (MQ=255)
                                                   ctggCCGGGCTCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGt                                                                                                                                           >  2:2241022/3‑107 (MQ=255)
                                                      gCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaaa                                                                                             <  1:1392796/150‑1 (MQ=255)
                                                        cGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGCAAGAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGa                                                                                           <  1:2332721/150‑1 (MQ=255)
                                                        cGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGa                                                                                           <  1:464649/150‑1 (MQ=255)
                                                        cGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGa                                                                                           <  1:864791/150‑1 (MQ=255)
                                                          aGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATg                                                                                         >  2:1156762/1‑150 (MQ=255)
                                                             ccGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAAc                                                                                      <  2:1322724/150‑1 (MQ=255)
                                                              cGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAAcc                                                                                     >  2:2370875/1‑150 (MQ=255)
                                                               ggCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTACGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCAGGTGGCTTTGCAAGAGGCTTTGCAACTTTGGGAGACACCCACACGGAAAGATGAACCg                                                                                    >  1:775358/1‑150 (MQ=17)
                                                                 ccAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCggg                                                                                  <  2:1570579/150‑1 (MQ=255)
                                                                             gggCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACcg                                                                      >  1:2094306/1‑150 (MQ=255)
                                                                                    gCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGAc                                                               >  2:2057852/1‑150 (MQ=255)
                                                                                          tAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTgg                                                         <  1:1156762/150‑1 (MQ=255)
                                                                                                 gTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGa                                                  <  2:2001298/150‑1 (MQ=255)
                                                                                                  tGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGag                                                 >  1:1011169/1‑150 (MQ=255)
                                                                                                    cGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACg                                                              >  2:1006412/1‑135 (MQ=255)
                                                                                                    cGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACg                                                              <  1:1006412/135‑1 (MQ=255)
                                                                                                     gATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCAt                                                    <  1:1539626/144‑1 (MQ=255)
                                                                                                     gATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCAt                                                    >  2:1539626/1‑144 (MQ=255)
                                                                                                     gATCCACCGCGAGGCCTCTTGATGCAAAACTTTGACGCCCAGGGGCCGCCTGGGCGTGGCTTTGCACAAAGCTTTCTAAATTAGGCAGACACCCACCCGAAAAGTCGAACCGGGCGCGCACCACCGCGTTCCCGGGGAGCTCGTGagcaa                                              >  2:777136/1‑150 (MQ=21)
                                                                                                       tCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGagcaagc                                            >  1:2203641/1‑150 (MQ=255)
                                                                                                          aCCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAAt                                         <  1:1661756/150‑1 (MQ=255)
                                                                                                           ccGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATc                                        >  1:1534281/1‑150 (MQ=255)
                                                                                                              cgAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCAcc                                     <  1:1201687/150‑1 (MQ=255)
                                                                                                                aGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaaa                                                                                             <  1:1552377/92‑1 (MQ=255)
                                                                                                                aGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaaa                                                                                             >  2:1552377/1‑92 (MQ=255)
                                                                                                                  tCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCaa                                   >  2:2114718/1‑148 (MQ=255)
                                                                                                                  tCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCaa                                   <  1:2114718/148‑1 (MQ=255)
                                                                                                                    ctctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGaaa                               >  2:1604059/1‑150 (MQ=255)
                                                                                                                     tctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                                                                                                             >  2:1353765/1‑71 (MQ=255)
                                                                                                                     tctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                                                                                                             <  1:1353765/71‑1 (MQ=255)
                                                                                                                       ttGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAAc                            >  2:167918/1‑150 (MQ=255)
                                                                                                                         gATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga                          <  1:1604059/150‑1 (MQ=255)
                                                                                                                         gATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga                          <  1:2057852/150‑1 (MQ=255)
                                                                                                                         gATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga                          <  1:1781243/150‑1 (MQ=255)
                                                                                                                           tGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACgaga                        <  2:1965379/150‑1 (MQ=255)
                                                                                                                           tGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACgaga                        <  2:2182751/150‑1 (MQ=255)
                                                                                                                              aaaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGacac                                                                                                        >  2:1556786/1‑67 (MQ=255)
                                                                                                                              aaaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGacac                                                                                                        <  1:1556786/67‑1 (MQ=255)
                                                                                                                               aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTc                    <  2:2188188/150‑1 (MQ=255)
                                                                                                                               aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGCGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTc                    >  2:1274762/1‑150 (MQ=255)
                                                                                                                                     ccaCGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                               <  2:971107/51‑1 (MQ=255)
                                                                                                                                     ccaCGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                               >  1:971107/3‑53 (MQ=255)
                                                                                                                                     ccaCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGccc              >  2:1612277/3‑150 (MQ=255)
                                                                                                                                     ccaCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGccc              >  2:407708/3‑150 (MQ=255)
                                                                                                                                     ccaCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTCGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGccc              >  2:49155/3‑150 (MQ=255)
                                                                                                                                      caCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGc                <  2:992101/146‑1 (MQ=255)
                                                                                                                                      caCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGc                >  1:992101/2‑147 (MQ=255)
                                                                                                                                      caCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCa             >  1:232120/2‑150 (MQ=255)
                                                                                                                                         gCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTgg                                                         <  2:14743/103‑1 (MQ=255)
                                                                                                                                         gCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTgg                                                         >  1:14743/1‑103 (MQ=255)
                                                                                                                                           ccAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaaa                                                                                             <  2:320191/65‑1 (MQ=255)
                                                                                                                                           ccAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaaa                                                                                             >  1:320191/1‑65 (MQ=255)
                                                                                                                                           ccAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGc        >  2:1150419/1‑150 (MQ=255)
                                                                                                                                            cAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCACTCGcc       <  2:1011169/150‑1 (MQ=255)
                                                                                                                                             aGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCct      <  1:1383217/150‑1 (MQ=255)
                                                                                                                                                 ggCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTAc  <  2:2094306/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTAC  >  NC_002947/4936527‑4936821

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