Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,541,825 Δ1 bp 5.7% coding (498/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,8250C. 5.7% 223.7 / 7.7 70coding (498/1503 nt)nhaBNa(+)/H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base . (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

GTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGG  >  NC_002947/4541678‑4541973
                                                                                                                                                    |                                                                                                                                                    
gttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTacca                                                                                                                                                    <  1:1869447/150‑1 (MQ=255)
 ttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTaccac                                                                                                                                                   >  2:777614/1‑150 (MQ=255)
 ttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTaccac                                                                                                                                                   >  2:1583802/1‑150 (MQ=255)
    tcACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACcgc                                                                                                                                                >  2:824243/1‑150 (MQ=255)
    tcACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACcgc                                                                                                                                                >  2:1061114/1‑150 (MQ=255)
       cGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTA‑CACcgcgta                                                                                                                                             <  1:897522/150‑2 (MQ=255)
               gCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCgg                                                                                                                                     <  1:998470/150‑1 (MQ=255)
               gCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCgg                                                                                                                                     <  2:104855/150‑1 (MQ=255)
                cTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCggg                                                                                                                                    <  2:607164/150‑1 (MQ=255)
                 tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCgggg                                                                                                                                   >  1:1515326/1‑150 (MQ=255)
                 tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCgggg                                                                                                                                   <  1:170782/150‑1 (MQ=255)
                 tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCAACGCGTCGCGTCCgggg                                                                                                                                   >  2:1152561/1‑150 (MQ=255)
                  cGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGc                                                                                                                                  <  1:1583802/150‑1 (MQ=255)
                      gTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAAt                                                                                                                              >  2:294824/1‑150 (MQ=255)
                        gCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATcc                                                                                                                            >  1:531900/1‑150 (MQ=255)
                              cAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCAACGCGTCGCGTCCGGGGCCAATCCGCGCGa                                                                                                                      <  1:1152561/150‑1 (MQ=255)
                               aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                     <  1:331867/150‑1 (MQ=255)
                               aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCAGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                     >  1:1485902/1‑150 (MQ=255)
                                  gCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGAc                                                                                                                  >  2:1079585/1‑150 (MQ=255)
                                   ccATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACcgcgtcgc                                                                                                                                           >  2:261871/1‑124 (MQ=255)
                                   ccATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACcgcgtcgc                                                                                                                                           <  1:261871/124‑1 (MQ=255)
                                   ccATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCCAGGAcc                                                                                                                 >  2:1178718/1‑149 (MQ=255)
                                     aTTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCCCCCCCCCCCCCCCCCCCTACGCGGT‑GTACCACCGCGTCGCGTCAGGGGCCAATCGGCGCGAGGACAgc                                                                                                               >  1:924925/1‑150 (MQ=255)
                                           gcgcTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTg                                                                                                         <  2:618460/150‑1 (MQ=255)
                                             gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACcgcgtcgc                                                                                                                                           >  2:72979/1‑114 (MQ=255)
                                             gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACcgcgtcgc                                                                                                                                           <  1:72979/114‑1 (MQ=255)
                                             gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGGGCACTGga                                                                                                       >  1:1298181/1‑150 (MQ=255)
                                             gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGga                                                                                                       <  1:294824/150‑1 (MQ=255)
                                                    ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCa                                                                                                            >  2:1170675/1‑138 (MQ=255)
                                                    ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCa                                                                                                            <  1:1170675/138‑1 (MQ=255)
                                                              tGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATcc                                                                                                                            <  2:1111741/112‑1 (MQ=255)
                                                              tGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATcc                                                                                                                            >  1:1111741/1‑112 (MQ=255)
                                                                  ggCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTAATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCa                                                                                  <  1:209875/150‑1 (MQ=255)
                                                                   gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAAt                                                                                                                              <  2:553687/105‑1 (MQ=255)
                                                                   gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAAt                                                                                                                              >  1:553687/1‑105 (MQ=255)
                                                                    cGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCAcc                                                                                <  2:918583/150‑1 (MQ=255)
                                                                          tATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAcc                                                                          <  1:29916/150‑1 (MQ=255)
                                                                              cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGa                                                                      >  1:688709/1‑150 (MQ=255)
                                                                              cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGAAGAACAGTGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGAGACCAGAAAGTAGCTCACTTGCACCGGGa                                                                      >  1:611900/1‑150 (MQ=255)
                                                                               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                     <  1:1079585/150‑1 (MQ=255)
                                                                               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                     >  1:1244637/1‑150 (MQ=255)
                                                                               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                     <  1:1140968/150‑1 (MQ=255)
                                                                               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                     <  1:824243/150‑1 (MQ=255)
                                                                                 gTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGa                                                                   >  1:1740125/1‑150 (MQ=255)
                                                                                     ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCcgc                                                                                                                          <  1:1348105/91‑1 (MQ=255)
                                                                                     ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCcgc                                                                                                                          >  2:1348105/1‑91 (MQ=255)
                                                                                     ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTc                                                               <  2:924925/150‑1 (MQ=255)
                                                                                             ccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGtt                                                       >  1:902991/1‑150 (MQ=255)
                                                                                             ccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGtt                                                       >  1:843190/1‑150 (MQ=255)
                                                                                              cTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTc                                                      >  1:603246/1‑150 (MQ=255)
                                                                                              cTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTc                                                      >  1:1241156/1‑150 (MQ=255)
                                                                                                 aCCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                                   <  2:688709/150‑1 (MQ=255)
                                                                                                   cGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTGTACGGGGT‑GTAACACGGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCgcgc                                                 >  1:823067/1‑150 (MQ=255)
                                                                                                     tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCt                                               <  1:1865836/150‑1 (MQ=255)
                                                                                                     tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCATCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTc                                                                                       <  1:872446/110‑1 (MQ=255)
                                                                                                     tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCATCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTc                                                                                       >  2:872446/1‑110 (MQ=255)
                                                                                                       aCGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTc                                             >  1:1470173/1‑150 (MQ=255)
                                                                                                        cggcggTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTcc                                            >  2:388511/1‑150 (MQ=255)
                                                                                                                   atcaGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctg                                 <  1:388511/150‑1 (MQ=255)
                                                                                                                      aGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgctg                              >  1:545726/1‑150 (MQ=255)
                                                                                                                              ggTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATg                       >  1:665735/1‑149 (MQ=255)
                                                                                                                              ggTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATg                       <  2:665735/149‑1 (MQ=255)
                                                                                                                                    cTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                     >  2:608603/1‑97 (MQ=255)
                                                                                                                                    cTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                     <  1:608603/97‑1 (MQ=255)
                                                                                                                                    cTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                                >  2:1789026/1‑134 (MQ=255)
                                                                                                                                    cTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                                <  1:1789026/134‑1 (MQ=255)
                                                                                                                                    cTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt                >  2:1604389/1‑150 (MQ=255)
                                                                                                                                    cTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt                <  2:531900/150‑1 (MQ=255)
                                                                                                                                     ttCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg               >  2:829815/1‑150 (MQ=255)
                                                                                                                                        gaCGCGGTGGTA‑CACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                                <  1:1739420/129‑1 (MQ=255)
                                                                                                                                        gaCGCGGTGGTA‑CACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                                >  2:1739420/2‑130 (MQ=255)
                                                                                                                                        gaCGCGGTGGTA‑CACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGt            >  1:389278/2‑150 (MQ=255)
                                                                                                                                          cgcgGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTAc          >  2:1113836/1‑150 (MQ=255)
                                                                                                                                                   accaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCa                                                          >  1:9683/1‑94 (MQ=255)
                                                                                                                                                   accaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCa                                                          <  2:9683/94‑1 (MQ=255)
                                                                                                                                                   accaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTgg  <  2:775127/150‑1 (MQ=255)
                                                                                                                                                    ccaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTC‑ACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg               <  2:915777/135‑1 (MQ=255)
                                                                                                                                                    ccaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTC‑ACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg               >  1:915777/1‑135 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGT‑GTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGG  >  NC_002947/4541678‑4541973

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: ‑