Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,468,078 Δ1 bp 5.9% coding (814/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,0780T. 5.9% 218.5 / 6.2 68coding (814/1047 nt)mltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (28/36);  new base . (2/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGCC  >  NC_002947/5467930‑5468225
                                                                                                                                                      |                                                                                                                                                   
aTCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGAt                                                                                                                                                     >  5:301691/1‑149 (MQ=255)
   gcgTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGcc                                                                                                                                                  >  6:96804/1‑149 (MQ=255)
    cgTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCg                                                                                                                                                 <  3:168741/149‑1 (MQ=255)
        gAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtg                                                                                                                                              <  3:27876/148‑1 (MQ=255)
         aaGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTTACCGGCTCGATGCCggtggt                                                                                                                                            <  1:12680/149‑1 (MQ=255)
         aaGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtggt                                                                                                                                            >  7:150442/1‑149 (MQ=255)
            ttcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                         <  8:82699/149‑1 (MQ=255)
            ttcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                         >  5:164025/1‑149 (MQ=255)
            ttcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                         <  6:138805/149‑1 (MQ=255)
            ttcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                         <  2:252408/149‑1 (MQ=255)
              cttcAGGCCCATGCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACGAGCCAAGCGTACGCTACTCCCCTACCGTTTTCACCGGCTCGATTCCGGTGGTCAgg                                                                                                                                        <  4:112678/148‑1 (MQ=255)
                  aGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGGCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCtt                                                                                                                                    <  2:19818/148‑1 (MQ=255)
                  aGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtc                                                                                                                                   <  8:201150/149‑1 (MQ=255)
                            aGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgcgcgc                                                                                                                         >  1:316157/1‑149 (MQ=255)
                              tACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaac                                                                                                                                         <  2:63540/131‑6 (MQ=255)
                              tACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaac                                                                                                                                         >  1:63540/1‑126 (MQ=255)
                                 tCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCAcc                                                                                                                    >  7:82281/1‑149 (MQ=255)
                                           tctcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAg                                                                                                          >  2:266096/1‑149 (MQ=255)
                                               gCCCTCGAGACGGAATGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGgcgc                                                                                                      <  5:96710/149‑1 (MQ=255)
                                               gCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGgc                                                                                                        <  3:285641/147‑1 (MQ=255)
                                                 ccTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGa                                                                                                    <  6:301691/149‑1 (MQ=255)
                                                  cTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGAc                                                                                                   >  8:65431/1‑149 (MQ=255)
                                                     gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGa                                                                                                <  7:65431/149‑1 (MQ=255)
                                                     gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGa                                                                                                >  8:11685/1‑149 (MQ=255)
                                                     gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGAGCGCTCACCTTCCACCCAGGCGCGACTGa                                                                                                >  2:201961/1‑149 (MQ=255)
                                                      agaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGAc                                                                                               <  6:170064/149‑1 (MQ=255)
                                                       gaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCGCACCTTCCACCCAGGCGCGACTGAc                                                                                               >  8:261489/1‑148 (MQ=255)
                                                          ggAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCa                                                                                                           <  8:281403/133‑1 (MQ=255)
                                                          ggAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCa                                                                                                           >  7:281403/1‑133 (MQ=255)
                                                                 ggTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACAccc                                                                                        >  6:251379/1‑145 (MQ=255)
                                                                 ggTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACAccc                                                                                        <  5:251379/145‑1 (MQ=255)
                                                                  gTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCgtc                                                                                   >  4:299096/1‑147 (MQ=255)
                                                                      ccGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCg                                                                                                                                                 <  1:58862/83‑1 (MQ=255)
                                                                      ccGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCg                                                                                                                                                 >  2:58862/1‑83 (MQ=255)
                                                                         gCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCAcc                                                                            <  1:266096/149‑1 (MQ=255)
                                                                           aGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCa                                                                          <  8:82281/149‑1 (MQ=255)
                                                                            gATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaa                                                                                                                                          >  5:91197/1‑80 (MQ=38)
                                                                            gATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaa                                                                                                                                          <  6:91197/84‑5 (MQ=38)
                                                                                 tcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCggg                                                                    <  6:124795/149‑3 (MQ=255)
                                                                                   gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTAACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTgg                                                                   >  1:229588/1‑149 (MQ=255)
                                                                                     gcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGgcgc                                                                <  2:229588/149‑1 (MQ=255)
                                                                                         cGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGa                                                            <  8:152830/149‑1 (MQ=255)
                                                                                                   aCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGTCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGc                                                  <  3:299096/149‑1 (MQ=255)
                                                                                                   aCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGc                                                  <  7:11685/149‑1 (MQ=255)
                                                                                                       aCCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACg                                              >  6:174896/1‑149 (MQ=255)
                                                                                                         cAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCt                                            >  3:99293/1‑149 (MQ=255)
                                                                                                            ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGcc                                         >  8:212049/1‑149 (MQ=255)
                                                                                                            ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                          >  1:159671/1‑148 (MQ=255)
                                                                                                            ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                          <  2:159671/148‑1 (MQ=255)
                                                                                                            ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGCTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGAGGCGCTTGAAGTAGCTCGCCACGCTGcc                                         >  5:325413/1‑149 (MQ=255)
                                                                                                                  gTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGc                                   <  5:96804/149‑1 (MQ=255)
                                                                                                                   tACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCGGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCa                                  >  5:325903/1‑149 (MQ=255)
                                                                                                                   tACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCa                                  <  5:174896/149‑1 (MQ=255)
                                                                                                                     cGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATc                                <  7:261489/149‑1 (MQ=255)
                                                                                                                      gCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCg                               >  1:81085/1‑149 (MQ=255)
                                                                                                                       cAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGt                              <  7:212049/149‑1 (MQ=255)
                                                                                                                             cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAgg                          <  6:164025/147‑1 (MQ=255)
                                                                                                                             cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgt                        >  3:108464/1‑149 (MQ=255)
                                                                                                                              ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgtt                       >  4:163888/1‑149 (MQ=255)
                                                                                                                              ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgt                        <  3:223831/148‑1 (MQ=255)
                                                                                                                              ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgt                        >  4:223831/1‑148 (MQ=255)
                                                                                                                                  aCCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttgtt                    >  8:16551/1‑149 (MQ=255)
                                                                                                                                    cGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTcc                  <  2:81085/149‑1 (MQ=255)
                                                                                                                                           aCCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGatat            <  4:99293/148‑1 (MQ=255)
                                                                                                                                            ccGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTg          <  1:201961/149‑1 (MQ=255)
                                                                                                                                             cGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGa         <  3:163888/149‑1 (MQ=255)
                                                                                                                                              ggCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGa         >  8:225353/1‑148 (MQ=255)
                                                                                                                                                    gATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGc   >  4:240286/1‑149 (MQ=255)
                                                                                                                                                     aTGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGcc  <  8:150442/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ATCGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGCC  >  NC_002947/5467930‑5468225

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