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

CGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGCC  >  NC_002947/5467932‑5468225
                                                                                                                                                   |                                                                                                                                                    
cgcgTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGc                                                                                                                                                    >  8:253625/1‑149 (MQ=255)
    tAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggt                                                                                                                                                <  3:295148/149‑1 (MQ=255)
     aGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAGCTCCCCTCCCGTTTTCACCGGCTCGAGGCCggtg                                                                                                                                               >  8:56502/1‑149 (MQ=255)
     aGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtg                                                                                                                                               >  4:127188/1‑149 (MQ=255)
         gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtggt                                                                                                                                             >  8:94541/1‑147 (MQ=255)
         gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTATCACCGGCATCGAGCCggtgaaaa                                                                                                                                           <  5:54637/149‑5 (MQ=255)
          ttcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                          <  5:25763/149‑1 (MQ=255)
                aGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCCTCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtc                                                                                                                                    <  4:279554/149‑1 (MQ=255)
                aGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACAGTTTTAACCGGCTCGATGCCGGTGGTCAGGCCTtc                                                                                                                                    <  2:137899/149‑1 (MQ=255)
                    ccATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtcg                                                                                                                                >  7:233612/1‑149 (MQ=255)
                     cATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCgg                                                                                                                               <  6:105070/149‑1 (MQ=255)
                          aGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgcgcgc                                                                                                                         >  5:185673/1‑149 (MQ=255)
                           gTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCAGGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCt                                                                                                                        >  6:201737/1‑149 (MQ=255)
                                   ggCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCa                                                                                                                                           >  6:82390/1‑123 (MQ=255)
                                   ggCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCa                                                                                                                                           <  5:82390/123‑1 (MQ=255)
                                    gCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCa                                                                                                               <  7:253625/149‑1 (MQ=255)
                                              cccTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGgcgc                                                                                                      <  4:123058/148‑1 (MQ=255)
                                              cccTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGgcgc                                                                                                      >  3:123058/1‑148 (MQ=255)
                                                   gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACt                                                                                                  <  8:233612/147‑1 (MQ=255)
                                                     gaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCAGCGCTCACCTTCCACCCAGGCGCGACTGAc                                                                                               <  1:196755/149‑1 (MQ=255)
                                                     gaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGACCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGAcc                                                                                              >  4:308870/1‑149 (MQ=255)
                                                        ggAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACcaca                                                                                           >  6:7879/1‑149 (MQ=255)
                                                                gTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGcc                                                                                   >  4:233507/1‑149 (MQ=255)
                                                                   aCCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGc                                                                                <  5:201737/149‑1 (MQ=255)
                                                                              atcGCGCAGCGAAACATGAACCGACCAGCCAAGCGTACGCAACTCGCCTACCGTTTTCACCGGCTCGATGCCGGAGGTCAGGCCTTCGTCGGCGCCCTCACCTTCCACAAAGGCCCGACTGACCACACCCTCGCAGGCCACCCAGACGt                                                                     >  1:233619/1‑149 (MQ=255)
                                                                                   gcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCAc                                                                                                              >  4:274740/1‑103 (MQ=255)
                                                                                   gcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCAc                                                                                                              <  3:274740/103‑1 (MQ=255)
                                                                                   gcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGgcgc                                                                >  7:300897/1‑149 (MQ=255)
                                                                                     aGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCtt                                                              >  3:128884/1‑149 (MQ=255)
                                                                                     aGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCtt                                                              <  3:127188/149‑1 (MQ=255)
                                                                                        gAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTg                                                             >  1:233893/1‑147 (MQ=255)
                                                                                             aTGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGc                                                      >  3:103896/1‑149 (MQ=255)
                                                                                              tGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCt                                                     >  8:80298/1‑149 (MQ=255)
                                                                                                acccGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTg                                                    <  4:2984/146‑1 (MQ=255)
                                                                                                aaCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTgg                                                   >  3:213667/1‑149 (MQ=255)
                                                                                                aaCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTgg                                                   >  3:107544/1‑149 (MQ=255)
                                                                                                  ccGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGc                                                  >  5:49407/1‑148 (MQ=255)
                                                                                                          ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGcc                                         <  6:49407/149‑1 (MQ=255)
                                                                                                           cAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGcc                                         >  4:164434/1‑148 (MQ=255)
                                                                                                            aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                          <  1:254154/146‑1 (MQ=255)
                                                                                                            aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                          >  2:254154/1‑146 (MQ=255)
                                                                                                            aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGa                                       <  3:233507/149‑1 (MQ=255)
                                                                                                              gCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATg                                     <  4:107544/149‑1 (MQ=255)
                                                                                                                   cGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCAt                                 <  8:300897/148‑1 (MQ=255)
                                                                                                                      aaCTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACt                                                                                                  <  2:257987/80‑1 (MQ=255)
                                                                                                                      aaCTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACt                                                                                                  >  1:257987/1‑80 (MQ=255)
                                                                                                                         tCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAgg                          >  4:210959/1‑149 (MQ=255)
                                                                                                                           cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAgg                          <  5:7879/147‑1 (MQ=255)
                                                                                                                           cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgt                        >  5:172473/1‑149 (MQ=255)
                                                                                                                           cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGACCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAggg                         <  3:308870/148‑1 (MQ=255)
                                                                                                                            ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                          <  5:69481/130‑1 (MQ=255)
                                                                                                                            ccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                          >  6:69481/1‑130 (MQ=255)
                                                                                                                             cTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttg                      >  8:244888/1‑149 (MQ=255)
                                                                                                                             cTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttg                      >  6:93192/1‑149 (MQ=255)
                                                                                                                              tACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGCTGGCAGCGTCAGGgttgg                     >  6:230711/1‑148 (MQ=255)
                                                                                                                              tACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttgt                     >  3:174410/1‑149 (MQ=255)
                                                                                                                               aCCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttgtt                    >  7:11529/1‑149 (MQ=255)
                                                                                                                                  gTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCa                 <  4:128884/149‑1 (MQ=255)
                                                                                                                                      ccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTGCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGa               <  1:281512/146‑1 (MQ=255)
                                                                                                                                      ccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATc                                <  2:318103/129‑1 (MQ=255)
                                                                                                                                      ccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATc                                >  1:318103/2‑130 (MQ=255)
                                                                                                                                      ccACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGa               >  2:281512/2‑147 (MQ=255)
                                                                                                                                       cACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGc                                   >  2:78535/1‑126 (MQ=255)
                                                                                                                                       cACCGGCATCGAGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGc                                   <  1:78535/126‑1 (MQ=255)
                                                                                                                                        aCCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt           <  3:210959/149‑1 (MQ=255)
                                                                                                                                         ccGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTg          <  5:230711/149‑1 (MQ=255)
                                                                                                                                               tCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGc                                                      >  6:26819/1‑100 (MQ=255)
                                                                                                                                               tCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGc                                                      <  5:26819/100‑1 (MQ=255)
                                                                                                                                                 gATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGGTATgg    <  4:103896/148‑1 (MQ=255)
                                                                                                                                                  aTGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGcc  >  7:308767/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGCC  >  NC_002947/5467932‑5468225

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