Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,764,801 A→G 6.7% D261G (GAC→GGC)  ctpA → carboxy‑terminal‑processing protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,764,8010AG 6.7% 212.3 / 4.4 74D261G (GAC→GGC) ctpAcarboxy‑terminal‑processing protease
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base G (3/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

TCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAATT  >  NC_002947/5764654‑5764949
                                                                                                                                                   |                                                                                                                                                    
tCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGAc                                                                                                                                                     <  6:329818/149‑1 (MQ=255)
   cccAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCAc                                                                                                                                                  >  5:150001/1‑149 (MQ=255)
    ccAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCACAACACCCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACt                                                                                                                                                 >  6:179374/1‑149 (MQ=255)
      aGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTc                                                                                                                                               <  5:4327/149‑1 (MQ=255)
           cAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGAc                                                                                                                                          >  6:322927/1‑149 (MQ=255)
               tCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCa                                                                                                                                        >  1:244248/1‑147 (MQ=255)
                 aaGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAg                                                                                                                                     <  6:150001/148‑1 (MQ=255)
                      cGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGa                                                                                                                               >  1:55722/1‑149 (MQ=255)
                       ggTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGAt                                                                                                                              >  7:278281/1‑149 (MQ=255)
                        gTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATc                                                                                                                             >  5:246907/1‑149 (MQ=255)
                           aCGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGt                                                                                                                           <  4:48321/148‑1 (MQ=255)
                           aCGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTc                                                                                                                          >  4:70723/1‑149 (MQ=255)
                             gAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                        >  3:300296/1‑149 (MQ=255)
                              aaGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGCCCAAAGGGCTGATCGTCTa                                                                                                                        >  8:248715/1‑148 (MQ=255)
                               aGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                                      >  7:537/1‑149 (MQ=255)
                               aGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                                      >  5:269079/1‑149 (MQ=255)
                                gTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTacac                                                                                                                     <  6:334794/149‑1 (MQ=255)
                                gTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTacac                                                                                                                     >  1:331721/1‑149 (MQ=255)
                                  gggCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACAcc                                                                                                                    >  6:52844/1‑148 (MQ=255)
                                        ggCCCTGGCCAAGCTGTGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGcc                                                                                                             <  8:11812/149‑1 (MQ=255)
                                        ggCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGcc                                                                                                             >  8:18291/1‑149 (MQ=255)
                                           ccTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGAGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGATGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATAGTCTACACCAAAGGCCGCa                                                                                                          >  4:286454/1‑149 (MQ=255)
                                                ccAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGcc                                                                                                     <  8:537/149‑1 (MQ=255)
                                                 cAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGcc                                                                                                     >  3:27762/1‑148 (MQ=255)
                                                    gCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACt                                                                                                 >  2:190027/1‑149 (MQ=255)
                                                        cgcAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCg                                                                                              <  6:246907/148‑1 (MQ=255)
                                                        cgcAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCg                                                                                              <  2:55722/148‑1 (MQ=255)
                                                             ggACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGc                                                                                        >  4:288705/1‑149 (MQ=255)
                                                               aCAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGTCGCATCGCCAACTCCGAGCTGCGt                                                                                      >  3:74219/1‑149 (MQ=255)
                                                                 aaCGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGtt                                                                                     >  4:237783/1‑148 (MQ=255)
                                                                    ggCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTtctc                                                                                 <  5:322927/149‑1 (MQ=255)
                                                                     gCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCg                                                                                <  2:331721/149‑1 (MQ=255)
                                                                     gCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCg                                                                                >  8:168072/1‑149 (MQ=255)
                                                                           agCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGAc                                                                          <  7:83923/149‑1 (MQ=255)
                                                                            gCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCTGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGAcc                                                                         >  2:347205/1‑149 (MQ=255)
                                                                                    gggTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGAt                                                                 <  3:237783/149‑1 (MQ=255)
                                                                                      gtggtgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGc                                                               <  8:278281/149‑1 (MQ=255)
                                                                                         gtgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCa                                                             <  7:18291/148‑1 (MQ=255)
                                                                                         gtgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAg                                                            <  6:269079/149‑1 (MQ=255)
                                                                                          tgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAgc                                                           >  7:175903/1‑149 (MQ=255)
                                                                                                  cTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGt                                                   >  5:127674/1‑149 (MQ=255)
                                                                                                  cTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGt                                                   >  8:125896/1‑149 (MQ=255)
                                                                                                    gcgcAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGc                                                 >  3:131357/1‑149 (MQ=255)
                                                                                                          caacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtgg                                           <  7:125896/149‑1 (MQ=255)
                                                                                                            acCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                         <  5:52844/149‑1 (MQ=255)
                                                                                                              ccAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGc                                                               >  1:29724/1‑125 (MQ=255)
                                                                                                              ccAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGc                                                               <  2:29724/125‑1 (MQ=255)
                                                                                                               cAGGCGGCGTGCTGCAATCGGCTGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTa                                      <  1:347205/149‑1 (MQ=255)
                                                                                                                   cggcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGa                                  <  3:288705/149‑1 (MQ=255)
                                                                                                                       gTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCaa                              <  4:131357/149‑1 (MQ=255)
                                                                                                                          ctgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACg                            <  4:300296/148‑1 (MQ=255)
                                                                                                                           tgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGt                          >  1:38741/1‑149 (MQ=255)
                                                                                                                            gcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCgg                                                                               <  3:221619/95‑1 (MQ=255)
                                                                                                                            gcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCgg                                                                               >  4:221619/1‑95 (MQ=255)
                                                                                                                            gcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTg                         >  3:96118/1‑149 (MQ=255)
                                                                                                                             cAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTgg                        <  3:70723/149‑1 (MQ=255)
                                                                                                                               aTCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGTCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCa                      <  4:74219/149‑1 (MQ=255)
                                                                                                                               aTCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCAGCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCa                      <  5:179374/149‑1 (MQ=255)
                                                                                                                                 cGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGGTGGTGGTACTGATCAACGGTGGCAGt                    >  4:224421/1‑149 (MQ=255)
                                                                                                                                  ggCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTg                   >  5:226650/1‑149 (MQ=255)
                                                                                                                                  ggCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCCCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGt                    <  2:244248/148‑1 (MQ=255)
                                                                                                                                   gCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCCGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTg                   >  5:20484/1‑148 (MQ=255)
                                                                                                                                     tcaggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcat                >  5:150485/4‑149 (MQ=255)
                                                                                                                                      gTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatc               <  6:169596/149‑1 (MQ=255)
                                                                                                                                       tGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcg              <  6:226650/149‑1 (MQ=255)
                                                                                                                                       aggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcg              >  7:215599/2‑149 (MQ=255)
                                                                                                                                        ggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                         <  6:303656/121‑1 (MQ=255)
                                                                                                                                        ggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                         >  5:303656/1‑121 (MQ=255)
                                                                                                                                         gAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcgg             >  7:191113/1‑148 (MQ=255)
                                                                                                                                           aGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcggcat          <  6:150485/149‑1 (MQ=255)
                                                                                                                                              ggCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcggcatcg        >  7:242227/1‑148 (MQ=255)
                                                                                                                                                ccGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGaa     <  7:168072/149‑1 (MQ=255)
                                                                                                                                                 cGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGaaa    >  6:33112/1‑149 (MQ=255)
                                                                                                                                                 cGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGaaa    >  6:233479/1‑149 (MQ=255)
                                                                                                                                                   aCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGAt                                 >  2:4799/1‑118 (MQ=255)
                                                                                                                                                   aCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAtt  <  3:286454/149‑1 (MQ=255)
                                                                                                                                                   aCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAtt  <  4:27762/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAATT  >  NC_002947/5764654‑5764949

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: