Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,183,450 T→C 6.7% Y312C (TAC→TGC)  mltF ← membrane‑bound lytic murein transglycosylase F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,183,4500TC 6.7% 216.9 / 4.3 74Y312C (TAC→TGC) mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base C (2/3);  total (36/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 = 7.89e-01

ATGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTA  >  NC_002947/1183303‑1183595
                                                                                                                                                      |                                                                                                                                                  
aTGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTa                                                                                                                                                   <  1:330264/149‑1 (MQ=255)
 tGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAg                                                                                                                                                  >  3:330759/1‑149 (MQ=255)
 tGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAg                                                                                                                                                  <  4:40132/149‑1 (MQ=255)
   aaGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGc                                                                                                                                                 >  2:15587/1‑148 (MQ=255)
      tACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAAt                                                                                                                                             <  5:121589/149‑1 (MQ=255)
       aCTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTg                                                                                                                                            >  1:112885/1‑149 (MQ=255)
            gcgcCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCtac                                                                                                                                       <  2:112885/149‑3 (MQ=255)
              gccgccCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                    <  6:96222/149‑1 (MQ=255)
                cgccCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCcag                                                                                                                                   >  8:290571/1‑148 (MQ=255)
                 gccCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACcagga                                                                                                                                 <  1:340042/149‑3 (MQ=255)
                 gccCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACcagg                                                                                                                                  >  2:340042/1‑147 (MQ=255)
                  cccTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCcagcag                                                                                                                                >  1:159894/1‑149 (MQ=255)
                    cTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgc                                                                                                                              <  5:297778/149‑1 (MQ=255)
                    cTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgc                                                                                                                              >  2:291665/1‑149 (MQ=255)
                         tGCTCTGCCTCGGGGCCAGTCGGGTCCACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAg                                                                                                                         <  2:333646/149‑1 (MQ=255)
                           ctctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTc                                                                                                                       <  3:154064/149‑1 (MQ=255)
                            tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgc                                                                                                                              >  6:284332/1‑141 (MQ=255)
                            tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgc                                                                                                                              <  5:284332/141‑1 (MQ=255)
                            tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTc                                                                                                                       >  1:157556/1‑148 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCCGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCca                                                                                                                                    >  6:162838/1‑131 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCca                                                                                                                                    <  5:162838/131‑1 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGa                                                                                                                     <  2:99232/146‑1 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCt                                                                                                                  >  6:144952/1‑149 (MQ=255)
                                   cGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCt                                                                                                               >  6:148597/1‑149 (MQ=255)
                                    gggTCCAGTCGGTTCGACACGCACATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGTTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCttgca                                                                                                              >  8:114434/1‑145 (MQ=255)
                                      gTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGa                                                                                                            >  6:10760/1‑149 (MQ=255)
                                         cAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCCAGCAGGCGCCAGTCGACCTGCTCGAcc                                                                                                          <  1:15587/149‑1 (MQ=255)
                                           gTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                       <  4:330759/149‑1 (MQ=255)
                                            tCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTgg                                                                                                      >  3:10944/1‑149 (MQ=255)
                                              ggTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCg                                                                                                    <  4:10944/149‑1 (MQ=255)
                                                ttCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCAGGTGGAGGTGACTTCCGGCTGCCACTTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCg                                                                                                    >  1:107564/1‑149 (MQ=255)
                                                  cGACACGCCGATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAg                                                                                                 <  7:240377/148‑1 (MQ=255)
                                                   gcaccGCCCAGGCCATGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGct                                                                                               >  1:237007/5‑149 (MQ=255)
                                                   gACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGc                                                                                                <  5:10760/148‑1 (MQ=255)
                                                    acacGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGctc                                                                                              >  4:301799/1‑149 (MQ=255)
                                                      acGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCg                                                                                                                      <  1:328773/123‑1 (MQ=255)
                                                      acGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCg                                                                                                                      >  2:328773/1‑123 (MQ=255)
                                                      acGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGctct                                                                                             <  8:41139/148‑1 (MQ=255)
                                                        gCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTg                                                                                           >  7:85346/1‑148 (MQ=255)
                                                                    ggCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTtctc                                                                                <  1:291665/147‑1 (MQ=255)
                                                                        gTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACtt                                                                          <  2:237007/149‑1 (MQ=255)
                                                                         tGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTg                                                                         <  5:144952/149‑1 (MQ=255)
                                                                           cgcTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTggg                                                                       >  2:219641/1‑149 (MQ=255)
                                                                           cgcTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCAGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTggg                                                                       >  4:90118/1‑149 (MQ=255)
                                                                            gcTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCTATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGc                                                                      <  4:324703/149‑1 (MQ=255)
                                                                                      catcatCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGtt                                                            <  1:219641/149‑1 (MQ=255)
                                                                                          atcaGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCTCATCTATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCa                                                         <  4:262020/148‑1 (MQ=255)
                                                                                          atcaGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAg                                                        <  5:148597/149‑1 (MQ=255)
                                                                                           tcaGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGAAGGCGCTGTTGCAgg                                                       >  4:227355/1‑149 (MQ=255)
                                                                                             aGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTg                                                     >  5:246676/1‑149 (MQ=255)
                                                                                                     cACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCg                                              >  3:215346/1‑148 (MQ=255)
                                                                                                      aCGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGaa                                            <  2:159894/149‑1 (MQ=255)
                                                                                                       cGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAg                                           <  4:195869/149‑1 (MQ=255)
                                                                                                             tC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGGATTCCTGTTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGgtg                                        >  1:210599/1‑147 (MQ=255)
                                                                                                                   ggAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGgcgc                                 >  8:231022/1‑149 (MQ=255)
                                                                                                                         gACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGa                              <  3:339323/146‑1 (MQ=255)
                                                                                                                         gACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGa                              >  4:339323/1‑146 (MQ=255)
                                                                                                                            ttCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAgg                                                       >  8:246322/1‑118 (MQ=255)
                                                                                                                            ttCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAgg                                                       <  7:246322/118‑1 (MQ=255)
                                                                                                                               cGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCaa                     <  2:157556/149‑1 (MQ=255)
                                                                                                                                ggCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAg                    <  8:85346/149‑1 (MQ=255)
                                                                                                                                     ccACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGt               >  3:265882/1‑149 (MQ=255)
                                                                                                                                      cacaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTc              <  7:231022/149‑1 (MQ=255)
                                                                                                                                      cacaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTc              >  2:75100/1‑149 (MQ=255)
                                                                                                                                       acaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAg                                     >  6:362060/1‑125 (MQ=255)
                                                                                                                                       acaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAg                                     <  5:362060/125‑1 (MQ=255)
                                                                                                                                       acaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCg             >  7:281275/1‑149 (MQ=255)
                                                                                                                                         aTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGa            >  6:27322/1‑148 (MQ=255)
                                                                                                                                             gATTCCTTGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATg        <  5:27322/149‑1 (MQ=255)
                                                                                                                                               ttCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGc      >  5:96257/1‑149 (MQ=255)
                                                                                                                                                 ccTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGcc     <  6:246676/148‑1 (MQ=255)
                                                                                                                                                     gTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgta  <  1:75100/147‑1 (MQ=255)
                                                                                                                                                     gTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgta  <  1:94751/147‑1 (MQ=255)
                                                                                                                                                     gTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgta  >  2:94751/1‑147 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ATGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTA  >  NC_002947/1183303‑1183595

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: