Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,673 G→C 8.1% 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,6730GC 8.1% 223.4 / 9.6 74intergenic (‑39/‑379)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base C (3/3);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTA  >  NC_002947/4936524‑4936820
                                                                                                                                                     |                                                                                                                                                   
cAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGAACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAggggg                                                                                                                                                     <  2:1551763/150‑1 (MQ=21)
 aaCATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGgc                                                                                                                                                    >  1:1433667/1‑150 (MQ=35)
   catcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGgcgc                                                                                                                                                  >  2:1266092/1‑150 (MQ=35)
        tcatTGCCAACATCACCAGCAGCACCAGCGGAACCCCAAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTggg                                                                                                                                             >  1:42177/1‑150 (MQ=255)
          atTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCg                                                                                                                                           <  2:1572259/150‑1 (MQ=255)
              ccAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGAGGGCGCCGGGGCGTGGc                                                                                                                                       >  2:606928/1‑150 (MQ=255)
               cAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGtcaaa                                                                                                                                      <  2:1021340/150‑5 (MQ=11)
                   aTCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGc                                                                                                                                  >  2:1415597/1‑150 (MQ=255)
                        cagcagCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaagaa                                                                                                                             <  2:667230/150‑1 (MQ=255)
                        cagcagCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCCGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaagaa                                                                                                                             >  1:234496/1‑150 (MQ=255)
                         agcagcAACAGCGGCACCCCCAGCTTGCCCCGGCCGATCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAGCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaagaag                                                                                                                            <  2:714568/150‑1 (MQ=255)
                             gggactgCGGCACCCCCACCTTGCCCCGGCCCAGACAGGCCCGGTTGGAGCGGGTGTTGCTGAGGAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCttt                                                                                                                        <  2:1755387/144‑1 (MQ=17)
                             gcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCttt                                                                                                                        >  1:1577943/1‑150 (MQ=255)
                              caACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTc                                                                                                                       >  1:810297/1‑150 (MQ=255)
                               aaCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTcc                                                                                                                      <  1:324823/150‑1 (MQ=255)
                                  aGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAAc                                                                                                                   <  2:1628520/150‑1 (MQ=255)
                                  aGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAAc                                                                                                                   <  2:1433667/150‑1 (MQ=255)
                                       cACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                              <  1:560301/150‑1 (MQ=255)
                                       cACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                              <  1:41514/150‑1 (MQ=255)
                                        aCCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGc                                                                                                             <  2:1564405/150‑1 (MQ=255)
                                          ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAg                                                                                                           >  2:966240/1‑150 (MQ=255)
                                          ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAg                                                                                                           >  2:968979/1‑150 (MQ=255)
                                                gTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGAGATCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACAccc                                                                                                     >  1:680719/1‑150 (MQ=255)
                                                cttGCCCCCGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGACCCACCGCGAGTCCTCTTGATGAAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTAGGCAGACAccc                                                                                                     >  2:652613/2‑150 (MQ=255)
                                                   gCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCcaca                                                                                                  <  1:1266092/150‑1 (MQ=255)
                                                    ccccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACttt                                                                                                                <  1:437176/135‑1 (MQ=255)
                                                    ccccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACttt                                                                                                                >  2:437176/1‑135 (MQ=255)
                                                     cctggCCGGGCTCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACg                                                                                                >  1:1389637/4‑150 (MQ=255)
                                                        ggCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaa                                                                                             <  1:898420/150‑1 (MQ=255)
                                                         gCCGAGCCCGAAAAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaaa                                                                                            <  1:390716/150‑1 (MQ=255)
                                                               cccGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGaa                                                                                      >  1:835555/1‑150 (MQ=255)
                                                                ccGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAAc                                                                                     <  1:722620/150‑1 (MQ=255)
                                                                  ggCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCg                                                                                   <  2:141168/150‑1 (MQ=255)
                                                                         ttgttgCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGc                                                                            >  1:1029057/1‑150 (MQ=255)
                                                                                 ggCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACcgc                                                                    <  2:1902290/150‑1 (MQ=255)
                                                                                   cGTTGCTGATTAACTGAGTGCGATGCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGTGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAACGATGAACCGGGCGAGCACCACCGCGt                                                                  <  2:965439/150‑1 (MQ=255)
                                                                                   cGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGt                                                                  <  2:1029057/150‑1 (MQ=255)
                                                                                    gTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAATCGGGCGAGCACCACCGCGtt                                                                 <  1:966240/150‑1 (MQ=255)
                                                                                    gTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGtt                                                                 <  1:968979/150‑1 (MQ=255)
                                                                                     ttGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaa                                                                                             <  1:754555/121‑1 (MQ=255)
                                                                                     ttGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaa                                                                                             >  2:754555/1‑121 (MQ=255)
                                                                                         tGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACg                                                             <  2:1798689/149‑1 (MQ=255)
                                                                                         tGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACg                                                             >  1:1798689/1‑149 (MQ=255)
                                                                                          gATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACggg                                                           >  1:999000/1‑150 (MQ=255)
                                                                                                 tGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCa                                                    >  1:1718320/1‑150 (MQ=255)
                                                                                                  gAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCAt                                                   <  2:999000/150‑1 (MQ=255)
                                                                                                        gATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGCACCGGGCGAGCACCACCGCGCTGACGGGTGGCTCATGagcaa                                             <  1:1030659/150‑1 (MQ=255)
                                                                                                        gATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                                          <  2:12593/137‑1 (MQ=255)
                                                                                                        gATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                                          >  1:12593/1‑137 (MQ=255)
                                                                                                        gATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGagcaa                                             <  2:1718320/150‑1 (MQ=255)
                                                                                                        gATACAACGAGAGTCTTCTTGTTCCAAAACTTTGACGCCCAGGGGGCGTCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACAAAACTATTGATACATGATCCGAGCGAGCACCACCGCGTTGACTTGTGGCTCATGagcat                                             <  2:266428/150‑2 (MQ=16)
                                                                                                         aTCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCAAAACCGCGTTGACGGGTGGCTCATGagcaag                                            <  1:817084/150‑1 (MQ=255)
                                                                                                             aCCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAAt                                        >  2:1288988/1‑150 (MQ=255)
                                                                                                                    gTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTgg                                                        >  1:1739514/1‑127 (MQ=255)
                                                                                                                    gTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTgg                                                        <  2:1739514/127‑1 (MQ=255)
                                                                                                                       ctctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGc                                                                                >  2:564315/1‑100 (MQ=255)
                                                                                                                       ctctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGc                                                                                <  1:564315/100‑1 (MQ=255)
                                                                                                                          ttGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAAc                           <  2:1389637/150‑1 (MQ=255)
                                                                                                                            gATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga                         <  1:664742/150‑1 (MQ=255)
                                                                                                                             aTGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACtag                        <  1:1508268/150‑1 (MQ=255)
                                                                                                                                cAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                              <  1:1130711/61‑1 (MQ=255)
                                                                                                                                cAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                              >  2:1130711/1‑61 (MQ=255)
                                                                                                                                cAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGAcc                     >  1:915133/1‑150 (MQ=255)
                                                                                                                                cAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGAcc                     >  2:570099/1‑150 (MQ=255)
                                                                                                                                 aaaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCt                    >  2:1618588/1‑150 (MQ=255)
                                                                                                                                 aaaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCt                    >  2:5740/1‑150 (MQ=255)
                                                                                                                                  aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGc                                                                                <  2:1388656/89‑1 (MQ=255)
                                                                                                                                  aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGc                                                                                >  1:1388656/1‑89 (MQ=255)
                                                                                                                                   aaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTcc                  >  1:1721068/1‑150 (MQ=255)
                                                                                                                                      tttGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGc               >  1:1393288/1‑150 (MQ=255)
                                                                                                                                        tGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGccc             <  2:810297/150‑1 (MQ=255)
                                                                                                                                          aCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGagcaa                                             >  2:353651/1‑116 (MQ=255)
                                                                                                                                          aCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGagcaa                                             <  1:353651/116‑1 (MQ=255)
                                                                                                                                          aCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCaa           >  1:156040/1‑150 (MQ=255)
                                                                                                                                          aCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCaa           >  2:554720/1‑150 (MQ=255)
                                                                                                                                           cGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAAt          <  1:1099255/150‑1 (MQ=255)
                                                                                                                                               cAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGcc      >  1:1816731/1‑150 (MQ=255)
                                                                                                                                                   gggCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTa  >  1:1378116/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTA  >  NC_002947/4936524‑4936820

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