Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,183,456 A→T 6.8% I310N (ATT→AAT)  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,4560AT 6.8% 216.9 / 4.9 73I310N (ATT→AAT) mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base T (2/3);  total (35/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 = 9.30e-01

TACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGG  >  NC_002947/1183309‑1183602
                                                                                                                                                      |                                                                                                                                                   
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                                                                                                                                     >  2:291665/1‑149 (MQ=255)
              cTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgc                                                                                                                                     <  5:297778/149‑1 (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                     >  2:75100/1‑149 (MQ=255)
                                                                                                                                cacaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTc                     <  7:231022/149‑1 (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:94751/147‑1 (MQ=255)
                                                                                                                                               gTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgta         <  1:75100/147‑1 (MQ=255)
                                                                                                                                               gTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgta         >  2:94751/1‑147 (MQ=255)
                                                                                                                                                 aGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGtagtag     >  6:292635/1‑149 (MQ=255)
                                                                                                                                                   ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCg   <  5:292635/149‑1 (MQ=255)
                                                                                                                                                    cAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCgg  <  2:210599/149‑1 (MQ=255)
                                                                                                                                                      aTGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTc                     <  7:196932/128‑1 (MQ=255)
                                                                                                                                                      aTGGATGCCAGCACGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTc                     >  8:196932/1‑128 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTC‑‑TTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGG  >  NC_002947/1183309‑1183602

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: