Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,541,821:1 +G 5.4% coding (494/1503 nt) nhaB → Na(+)/H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,541,8211.G 5.4% 308.0 / 10.5 74coding (494/1503 nt)nhaBNa(+)/H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base . (38/32);  new base G (2/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

CCTGTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTAC  >  NC_002947/4541675‑4541965
                                                                                                                                                   |                                                                                                                                                
cctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTa                                                                                                                                               <  2:801815/150‑1 (MQ=255)
   gttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTacca                                                                                                                                            <  1:1869447/150‑1 (MQ=255)
    ttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTaccac                                                                                                                                           >  2:777614/1‑150 (MQ=255)
    ttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTaccac                                                                                                                                           >  2:1583802/1‑150 (MQ=255)
       tcACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgc                                                                                                                                        >  2:1061114/1‑150 (MQ=255)
       tcACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgc                                                                                                                                        >  2:824243/1‑150 (MQ=255)
          cGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgta                                                                                                                                     <  1:897522/150‑2 (MQ=255)
                  gCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgg                                                                                                                             <  2:104855/150‑1 (MQ=255)
                  gCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgg                                                                                                                             <  1:998470/150‑1 (MQ=255)
                   cTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCggg                                                                                                                            <  2:607164/150‑1 (MQ=255)
                    tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                           >  1:1515326/1‑150 (MQ=255)
                    tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                           <  1:170782/150‑1 (MQ=255)
                    tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCAACGCGTCGCGTCCgggg                                                                                                                           >  2:1152561/1‑150 (MQ=255)
                     cGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGc                                                                                                                          <  1:1583802/150‑1 (MQ=255)
                         gTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAAt                                                                                                                      >  2:294824/1‑150 (MQ=255)
                           gCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATcc                                                                                                                    >  1:531900/1‑150 (MQ=255)
                                 cAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCAACGCGTCGCGTCCGGGGCCAATCCGCGCGa                                                                                                              <  1:1152561/150‑1 (MQ=255)
                                  aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                             <  1:331867/150‑1 (MQ=255)
                                  aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCAGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                             >  1:1485902/1‑150 (MQ=255)
                                     gCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGAc                                                                                                          >  2:1079585/1‑150 (MQ=255)
                                      ccATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgc                                                                                                                                   <  1:261871/124‑1 (MQ=255)
                                      ccATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgc                                                                                                                                   >  2:261871/1‑124 (MQ=255)
                                      ccATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCCAGGAcc                                                                                                         >  2:1178718/1‑149 (MQ=255)
                                        aTTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCCCCCCCCCCCCCCCCCCCTACGCGGTGTACCACCGCGTCGCGTCAGGGGCCAATCGGCGCGAGGACAgc                                                                                                       >  1:924925/1‑150 (MQ=255)
                                              gcgcTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTg                                                                                                 <  2:618460/150‑1 (MQ=255)
                                                gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgc                                                                                                                                   <  1:72979/114‑1 (MQ=255)
                                                gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgc                                                                                                                                   >  2:72979/1‑114 (MQ=255)
                                                gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGGGCACTGga                                                                                               >  1:1298181/1‑150 (MQ=255)
                                                gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGga                                                                                               <  1:294824/150‑1 (MQ=255)
                                                       ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCa                                                                                                    <  1:1170675/138‑1 (MQ=255)
                                                       ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCa                                                                                                    >  2:1170675/1‑138 (MQ=255)
                                                                 tGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATcc                                                                                                                    >  1:1111741/1‑112 (MQ=255)
                                                                 tGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATcc                                                                                                                    <  2:1111741/112‑1 (MQ=255)
                                                                     ggCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTAATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCa                                                                          <  1:209875/150‑1 (MQ=255)
                                                                      gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAAt                                                                                                                      >  1:553687/1‑105 (MQ=255)
                                                                      gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAAt                                                                                                                      <  2:553687/105‑1 (MQ=255)
                                                                       cGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCAcc                                                                        <  2:918583/150‑1 (MQ=255)
                                                                             tATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAcc                                                                  <  1:29916/150‑1 (MQ=255)
                                                                                 cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGa                                                              >  1:688709/1‑150 (MQ=255)
                                                                                 cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGAAGAACAGTGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGAGACCAGAAAGTAGCTCACTTGCACCGGGa                                                              >  1:611900/1‑150 (MQ=255)
                                                                                  gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                             <  1:824243/150‑1 (MQ=255)
                                                                                  gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                             >  1:1244637/1‑150 (MQ=255)
                                                                                  gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                             <  1:1079585/150‑1 (MQ=255)
                                                                                  gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                             <  1:1140968/150‑1 (MQ=255)
                                                                                    gTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGa                                                           >  1:1740125/1‑150 (MQ=255)
                                                                                        ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCcgc                                                                                                                  >  2:1348105/1‑91 (MQ=255)
                                                                                        ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCcgc                                                                                                                  <  1:1348105/91‑1 (MQ=255)
                                                                                        ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTc                                                       <  2:924925/150‑1 (MQ=255)
                                                                                                ccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGtt                                               >  1:902991/1‑150 (MQ=255)
                                                                                                ccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGtt                                               >  1:843190/1‑150 (MQ=255)
                                                                                                 cTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTc                                              >  1:1241156/1‑150 (MQ=255)
                                                                                                 cTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTc                                              >  1:603246/1‑150 (MQ=255)
                                                                                                    aCCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                           <  2:688709/150‑1 (MQ=255)
                                                                                                      cGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTGTACGGGGTGTAACACGGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCgcgc                                         >  1:823067/1‑150 (MQ=255)
                                                                                                        tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCt                                       <  1:1865836/150‑1 (MQ=255)
                                                                                                        tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCATCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTc                                                                               >  2:872446/1‑110 (MQ=255)
                                                                                                        tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCATCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTc                                                                               <  1:872446/110‑1 (MQ=255)
                                                                                                          aCGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTc                                     >  1:1470173/1‑150 (MQ=255)
                                                                                                           cggcggTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTcc                                    >  2:388511/1‑150 (MQ=255)
                                                                                                                      atcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctg                         <  1:388511/150‑1 (MQ=255)
                                                                                                                         aGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgctg                      >  1:545726/1‑150 (MQ=255)
                                                                                                                                 ggTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATg               >  1:665735/1‑149 (MQ=255)
                                                                                                                                 ggTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATg               <  2:665735/149‑1 (MQ=255)
                                                                                                                                       cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                             >  2:608603/1‑97 (MQ=255)
                                                                                                                                       cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                             <  1:608603/97‑1 (MQ=255)
                                                                                                                                       cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                        >  2:1789026/1‑134 (MQ=255)
                                                                                                                                       cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                        <  1:1789026/134‑1 (MQ=255)
                                                                                                                                       cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt        <  2:531900/150‑1 (MQ=255)
                                                                                                                                       cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt        >  2:1604389/1‑150 (MQ=255)
                                                                                                                                        ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg       >  2:829815/1‑150 (MQ=255)
                                                                                                                                           gaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                        <  1:1739420/129‑1 (MQ=255)
                                                                                                                                           gaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                        >  2:1739420/2‑130 (MQ=255)
                                                                                                                                           gaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGt    >  1:389278/2‑150 (MQ=255)
                                                                                                                                             cgcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTAc  >  2:1113836/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCTGTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTAC  >  NC_002947/4541675‑4541965

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