Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,468,121:1 +A 11.5% coding (819/1047 nt) mltB ← membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,468,1211.A11.5% 294.7 / 37.8 87coding (819/1047 nt)mltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base . (40/37);  new base A (7/3);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

AGTGATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATG  >  NC_002947/5467974‑5468270
                                                                                                                                                    |                                                                                                                                                     
aGTGATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTc                                                                                                                                                     <  1:786885/150‑1 (MQ=255)
  tGATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGa                                                                                                                                                   >  1:95974/1‑150 (MQ=255)
   gATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCGGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGAt                                                                                                                                                  >  1:254311/1‑150 (MQ=255)
         gTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggt                                                                                                                                            <  1:1756709/150‑1 (MQ=255)
          tAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtg                                                                                                                                           <  2:1631483/150‑1 (MQ=255)
          tAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtg                                                                                                                                           >  2:2288370/1‑150 (MQ=255)
             aaGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTc                                                                                                                                        >  1:220314/1‑150 (MQ=255)
              aGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaa                                                                                                                                       <  1:1445819/150‑5 (MQ=255)
               gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                      <  1:1908186/150‑1 (MQ=255)
               gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                      <  1:1889010/150‑1 (MQ=255)
               gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTCCACCGGCATCGAGCCggtgaaaac                                                                                                                                      <  2:1290989/150‑6 (MQ=255)
               gTTCTTCAGGCCCATCCAGTACGCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCACTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                      <  2:1099774/150‑1 (MQ=255)
                 tcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGc                                                                                                                                    >  2:214195/1‑150 (MQ=255)
                 tcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGc                                                                                                                                    <  1:214195/150‑1 (MQ=255)
                        gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtc                                                                                                                             >  2:441990/1‑150 (MQ=255)
                        gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtc                                                                                                                             >  1:825547/1‑150 (MQ=255)
                           cATCCACTTCTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgc                                                                                                                          <  1:1606188/150‑1 (MQ=255)
                            aTCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgcg                                                                                                                         <  2:332869/150‑1 (MQ=255)
                             tCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                      >  2:586329/1‑136 (MQ=255)
                             tCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                      <  1:586329/136‑1 (MQ=255)
                              ccAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgcgcg                                                                                                                       >  2:1763108/1‑150 (MQ=255)
                                aGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCt                                                                                                                     >  1:2405545/1‑150 (MQ=255)
                                     tCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCt                                                                                                                >  1:2110240/1‑150 (MQ=255)
                                           ccGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCAcc                                                                                                          >  1:516375/1‑150 (MQ=255)
                                               tctcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGATGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAgg                                                                                                      >  1:124234/1‑150 (MQ=255)
                                                     ccTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGAc                                                                                                >  1:1674090/1‑150 (MQ=255)
                                                      cTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCACCGAATCATGGACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACt                                                                                               >  1:2256184/1‑150 (MQ=255)
                                                       tCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTg                                                                                              >  1:434414/1‑150 (MQ=255)
                                                        cGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGa                                                                                             <  1:1195625/150‑1 (MQ=255)
                                                         gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGAc                                                                                            <  2:2405545/150‑1 (MQ=255)
                                                           gaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACca                                                                                          <  1:2436478/150‑1 (MQ=255)
                                                                      gTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCg                                                                               <  2:1109299/150‑1 (MQ=255)
                                                                       tAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCgg                                                                              >  1:2013986/1‑150 (MQ=255)
                                                                             gCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCCACTCACCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCCCTCCCCTTCCACCCAGGCGCGACTGACCCCACCCTCGCCGGCCAccc                                                                        >  2:462696/1‑150 (MQ=255)
                                                                                 atcatcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGcc                                                                    >  2:1283687/1‑150 (MQ=255)
                                                                                  tcatcGCGCGGCGAATCATGCACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTAACCGGCTCGCTGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACGGACCACACCCTCGCCGGCCACCCAGCCg                                                                   >  1:584187/1‑150 (MQ=255)
                                                                                       gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGgcg                                                              >  1:259726/1‑150 (MQ=255)
                                                                                           aGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTg                                                          >  1:2078563/1‑150 (MQ=255)
                                                                                            gCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCCCCCAGGCGCGACTGACCacac                                                                                       <  2:931073/120‑1 (MQ=255)
                                                                                            gCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCacac                                                                                       >  1:931073/1‑120 (MQ=255)
                                                                                                aTCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTa                                                     <  2:1148580/150‑1 (MQ=255)
                                                                                                 tCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGGCTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg                                                    <  2:254311/150‑1 (MQ=255)
                                                                                                       aCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGcc                                              >  1:1674746/1‑150 (MQ=255)
                                                                                                             cAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAg                                                       >  1:1522180/1‑135 (MQ=255)
                                                                                                             cAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAg                                                       <  2:1522180/135‑1 (MQ=255)
                                                                                                             cAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTg                                        >  2:1580480/1‑150 (MQ=255)
                                                                                                                ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCTCTTGAAGTAGCTGGCCACGCTGCCg                                     >  1:1818267/1‑150 (MQ=255)
                                                                                                                ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTAAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTGGCTGGCCACGCTGCCg                                     >  1:507245/1‑150 (MQ=255)
                                                                                                                  aaGCGTACGCAACTCCCCTACCGTTTTCGCCGGCTCGATGCCGGTGGTCAGGCCGTCGCCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                   <  2:1876719/150‑1 (MQ=255)
                                                                                                                  aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTGGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                   <  2:879451/150‑1 (MQ=255)
                                                                                                                  aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                   <  1:441990/150‑1 (MQ=255)
                                                                                                                   aGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATg                                  <  2:825547/150‑1 (MQ=255)
                                                                                                                   aGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATg                                  <  2:1674090/150‑1 (MQ=255)
                                                                                                                    gCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATgg                                 >  2:2415027/1‑150 (MQ=255)
                                                                                                                     cGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCt                                         <  2:2302298/141‑1 (MQ=255)
                                                                                                                     cGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCt                                         >  1:2302298/1‑141 (MQ=255)
                                                                                                                        aCGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATc                             <  2:1674746/150‑1 (MQ=255)
                                                                                                                          gCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGt                           >  1:1393950/1‑150 (MQ=255)
                                                                                                                             aCTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTgg                                                >  2:810384/1‑126 (MQ=255)
                                                                                                                             aCTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTgg                                                <  1:810384/126‑1 (MQ=255)
                                                                                                                             aCTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAg                        <  1:1763108/150‑1 (MQ=255)
                                                                                                                              cTCCCCTACCGTTTTCACCGGCTCGATGCCGGTTGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAgg                       <  1:2443515/150‑1 (MQ=255)
                                                                                                                               tCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                       >  1:1521462/1‑133 (MQ=255)
                                                                                                                               tCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                       <  2:1521462/133‑1 (MQ=255)
                                                                                                                               tCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAggg                      >  1:1617352/1‑150 (MQ=255)
                                                                                                                                  ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCa                                                                                                            <  1:2169605/61‑1 (MQ=255)
                                                                                                                                  ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCa                                                                                                            >  2:2169605/1‑61 (MQ=255)
                                                                                                                                   cTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttgt                  >  1:2418566/1‑150 (MQ=255)
                                                                                                                                    tACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCATCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttgtt                 >  2:674773/1‑150 (MQ=255)
                                                                                                                                      ccGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAg                                                       >  1:1701231/1‑110 (MQ=255)
                                                                                                                                      ccGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAg                                                       <  2:1701231/110‑1 (MQ=255)
                                                                                                                                       cGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCa              <  1:2415027/150‑1 (MQ=255)
                                                                                                                                        gTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAg                        >  2:9443/1‑139 (MQ=255)
                                                                                                                                        gTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAg                        <  1:9443/139‑1 (MQ=255)
                                                                                                                                        gTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAg             >  2:2376552/1‑150 (MQ=255)
                                                                                                                                        ctgaccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAg             >  2:2121389/6‑150 (MQ=255)
                                                                                                                                        ctgaccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAg             >  1:1176649/6‑150 (MQ=255)
                                                                                                                                          gaccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCacac                                                                                       <  1:888854/71‑1 (MQ=38)
                                                                                                                                          gaccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCacac                                                                                       >  2:888854/4‑74 (MQ=38)
                                                                                                                                             cACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt        >  2:944429/1‑150 (MQ=255)
                                                                                                                                             cACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt        >  2:2279596/1‑150 (MQ=255)
                                                                                                                                             cACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt        >  2:1630352/1‑150 (MQ=255)
                                                                                                                                             cACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt        >  2:1026052/1‑150 (MQ=255)
                                                                                                                                               ccGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGa      <  2:2418566/150‑1 (MQ=255)
                                                                                                                                                cGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGat     <  1:1580480/150‑1 (MQ=255)
                                                                                                                                                cGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGat     <  2:2110240/150‑1 (MQ=255)
                                                                                                                                                cGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGat     <  2:124234/150‑1 (MQ=255)
                                                                                                                                                   cTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATg  <  1:2288370/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGTGATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATG  >  NC_002947/5467974‑5468270

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: