Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,468,758 T→G 6.9% D45A (GAC→GCC)  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,7580TG 6.9% 201.4 / 4.9 69D45A (GAC→GCC) mltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base G (3/2);  total (37/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.14e-01

GGTTTGACCCGCTCGGCCGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCGGTGCGTAT  >  NC_002947/5468611‑5468904
                                                                                                                                                    |                                                                                                                                                  
ggTTTGACCCGCTCGGCCGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCcgtc                                                                                                                                                   <  8:134684/149‑1 (MQ=255)
    tGACCCGCTCGGCCGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTa                                                                                                                                                <  5:203472/148‑1 (MQ=255)
       cccGCTCGGCCGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTc                                                                                                                                             <  7:206796/148‑1 (MQ=255)
        ccGCTCGGCCGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGc                                                                                                                                           >  6:340126/1‑149 (MQ=255)
          gCTCGGCCGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCg                                                                                                                                         <  6:240075/149‑1 (MQ=255)
                ccGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGGCGAATACGACGGctcgc                                                                                                                                   <  8:77586/149‑5 (MQ=255)
                    gcgcGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCTGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCa                                                                                                                               <  8:74051/149‑1 (MQ=255)
                         agaTGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGt                                                                                                                          <  6:25096/149‑1 (MQ=255)
                          gaTGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTa                                                                                                                         >  5:135564/1‑149 (MQ=255)
                            tGGCGTCGATAATCGACTGCGTGCGCGGGACCTCACCATCCTACTCCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTAc                                                                                                                        <  1:287508/149‑1 (MQ=255)
                              gCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAgc                                                                                                                     <  3:213850/149‑1 (MQ=255)
                                gTCGAGAATCGACTGATTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAgcgc                                                                                                                   >  8:171192/1‑149 (MQ=255)
                                 tCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGcc                                                                                                                  >  1:59328/1‑149 (MQ=255)
                                     gaATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAAc                                                                                                              <  7:299201/149‑1 (MQ=255)
                                       aTCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACa                                                                                                             <  8:302732/148‑1 (MQ=255)
                                           aCTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGccc                                                                                                        <  2:30935/149‑1 (MQ=255)
                                             tGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCAGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCcac                                                                                                      >  7:112291/1‑149 (MQ=255)
                                               cTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCcaccg                                                                                                    >  5:234608/1‑149 (MQ=255)
                                                           cTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAg                                                                                        >  6:94900/1‑149 (MQ=255)
                                                           cTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAg                                                                                        >  6:152561/1‑149 (MQ=255)
                                                            tCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAgg                                                                                       >  2:40560/1‑149 (MQ=255)
                                                              aCGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGgcg                                                                                     >  3:274157/1‑149 (MQ=255)
                                                                gAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGca                                                                                   >  8:129889/1‑149 (MQ=255)
                                                                   ac‑acacCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACaa                                                                                >  1:117006/1‑149 (MQ=255)
                                                                      acacCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCcaccgca                                                                                                  >  7:25922/1‑129 (MQ=255)
                                                                      acacCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCcaccgca                                                                                                  <  8:25922/129‑1 (MQ=255)
                                                                        acCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACggg                                                                            >  1:228740/1‑149 (MQ=255)
                                                                             tCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGcc                                                                       <  1:40560/149‑1 (MQ=255)
                                                                                  tGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGt                                                                   >  3:158056/1‑148 (MQ=255)
                                                                                   gCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTa                                                                 >  7:88755/1‑149 (MQ=255)
                                                                                   gCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTa                                                                 >  6:247093/1‑149 (MQ=255)
                                                                                    cTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTAc                                                                <  5:94900/149‑1 (MQ=255)
                                                                                    cTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTAc                                                                <  8:88755/149‑1 (MQ=255)
                                                                                       aGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACg                                                               <  7:129889/147‑1 (MQ=255)
                                                                                            cAAAAGCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCt                                                        >  7:310697/1‑149 (MQ=255)
                                                                                             aaaaCCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCTCCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCCCTGCtt                                                       >  3:86583/1‑149 (MQ=255)
                                                                                             aaaaCCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCtt                                                       >  8:280451/1‑149 (MQ=255)
                                                                                             aaaaCCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCtt                                                       >  6:262107/1‑149 (MQ=255)
                                                                                              aaaCCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCa                                                                     <  3:86071/134‑1 (MQ=255)
                                                                                              aaaCCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCa                                                                     >  4:86071/1‑134 (MQ=255)
                                                                                                 ccATAGTAGCTGCTCATTTCGACAACAAACGCGGCTACCTCGGGCGAGCCGGCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACCGTCAGCACAGTTACGCACTGCTTGCAt                                                   <  6:8206/149‑1 (MQ=255)
                                                                                                      gTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGaaa                                               <  5:247093/148‑1 (MQ=255)
                                                                                                       tCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAAtt                                             >  5:169414/1‑149 (MQ=255)
                                                                                                           ggCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTc                                         <  5:340126/149‑1 (MQ=255)
                                                                                                            gCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCa                                        <  6:135564/149‑1 (MQ=255)
                                                                                                            gCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCa                                        >  2:69410/1‑149 (MQ=255)
                                                                                                                aTTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCtt                                    >  1:194864/1‑149 (MQ=255)
                                                                                                                aTTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCt                                     >  6:134854/1‑148 (MQ=255)
                                                                                                                 tttCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTa                                   >  8:289173/1‑149 (MQ=255)
                                                                                                                   tCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTAt                                  <  5:262107/148‑1 (MQ=255)
                                                                                                                   tCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTAt                                  <  7:280451/148‑1 (MQ=255)
                                                                                                                       caacaaACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGc                                                         >  1:297707/1‑121 (MQ=255)
                                                                                                                       caacaaACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGc                                                         <  2:297707/121‑1 (MQ=255)
                                                                                                                        aacaaACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGtt                                                                  <  3:212633/111‑1 (MQ=255)
                                                                                                                        aacaaACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGtt                                                                  >  4:212633/1‑111 (MQ=255)
                                                                                                                            aaCGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGTCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTgc                        <  7:212989/149‑1 (MQ=255)
                                                                                                                             aCGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTgcg                       >  7:119202/1‑149 (MQ=255)
                                                                                                                             aCGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTgcg                       <  8:310697/149‑1 (MQ=255)
                                                                                                                             aCGCGGCTACCTCGGCCGAGCCGTCGAATTCGCCGGCATGGGCCAGCTGTACAGAGCCGAACAGGCCCACCCCACCGAACCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTgc                        >  3:202941/1‑148 (MQ=255)
                                                                                                                                   cTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTgcg                       >  7:89990/1‑143 (MQ=255)
                                                                                                                                   cTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTgcg                       <  8:89990/143‑1 (MQ=255)
                                                                                                                                      ccTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCAt               <  2:117006/148‑1 (MQ=255)
                                                                                                                                         cGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATACAAGGCGCACAACGGGCAGCCCAGTTACGAACTGCTTGCATGAAATTGTTAACCTAATTCAAACCTGCGCGATCCAttt             >  4:132134/1‑147 (MQ=255)
                                                                                                                                          gggCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTg            <  2:194864/147‑1 (MQ=255)
                                                                                                                                          gcccatgCCGGCGAATACGACGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCAGCGCGATCCAATTGc           >  5:129250/7‑148 (MQ=255)
                                                                                                                                           cccatgCCGGCGAATACGACGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCgg         >  4:322541/6‑149 (MQ=255)
                                                                                                                                            gCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCGGt        <  4:274157/149‑1 (MQ=255)
                                                                                                                                              gAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCGGTGc      >  6:135395/1‑149 (MQ=255)
                                                                                                                                              atgCCGGCGAATACGACGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCGGTGc      >  6:339398/3‑149 (MQ=255)
                                                                                                                                                gCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCGGTGCGt    <  7:171192/149‑1 (MQ=255)
                                                                                                                                                  cgtcgtATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCGGTGCGTAt  >  5:134176/1‑149 (MQ=255)
                                                                                                                                                   gtcgtATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTa                                   >  3:69485/1‑115 (MQ=255)
                                                                                                                                                   gtcgtATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTa                                   <  4:69485/115‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTTTGACCCGCTCGGCCGGGCGCGAGATGGCGTCGAGAATCGACTGCTTGCGCTGCACCTCACGAAAC‑ACACCCATCAGTTGCTCAGGCGCAAAACCATAGTCGCGGCTCATTTCGCCAACAAACGCGGCTACCTCGGGCGAGCCGTCGTATTCGCCGGCATGGGCCAGCTGTACAGCGCCGAACAGGCCCACCGCACCGATCCAAGGCGCACAACGGGCAGCCCAGTTACGCACTGCTTGCATGAAATTGTTCACCTTATTCAAACCTGCGCGATCCATTTGCGGTGCGTAT  >  NC_002947/5468611‑5468904

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