Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,719,614 A→G 6.9% Y133H (TAC→CAC)  PP_1515 ← Transcriptional regulator, tetR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,719,6140AG 6.9% 164.2 / 2.9 58Y133H (TAC→CAC) PP_1515Transcriptional regulator, tetR family
Reads supporting (aligned to +/- strand):  ref base A (23/31);  new base G (2/2);  total (25/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 = 8.82e-01

GGTAATCTGGCGCCCCCTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTACCCCCTCGGGATAC  >  NC_002947/1719467‑1719753
                                                                                                                                                   |                                                                                                                                           
ggTAATCTGGCGCCCCCTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTAc                                                                                                                                            >  8:41989/1‑149 (MQ=255)
  tAATCTGGCGCCCCCTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTc                                                                                                                                          <  7:41989/149‑1 (MQ=255)
    aTCTGGCGCCCCCTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCg                                                                                                                                         <  4:317732/148‑1 (MQ=255)
    aTCTGGCGCCCCCTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGc                                                                                                                                        <  7:59317/149‑1 (MQ=255)
         gcgcCCCCTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCaaaa                                                                                                                                   >  3:255883/1‑149 (MQ=255)
              cccTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCgg                                                                                                                              <  7:165760/149‑1 (MQ=255)
               ccTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCggg                                                                                                                             >  1:324307/1‑149 (MQ=255)
                   gACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTc                                                                                                                           >  8:17446/1‑147 (MQ=255)
                      cAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGc                                                                                                                      >  6:119041/1‑149 (MQ=255)
                      cAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGc                                                                                                                      >  8:261802/1‑149 (MQ=255)
                         gCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGc                                                                                                                   <  3:289740/149‑1 (MQ=255)
                          cAGAAGAAGTGCTCCGCCGCCCTGAATGGGTTGTCGTTGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGcc                                                                                                                  <  8:175700/149‑1 (MQ=255)
                           agaagaAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGccc                                                                                                                 >  3:123839/1‑149 (MQ=255)
                             aagaagTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTg                                                                                                               >  3:280183/1‑149 (MQ=255)
                              agaagTGCTCCGCCGCCCTCAAAGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTgg                                                                                                              >  7:201383/1‑149 (MQ=255)
                               gaagTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTggg                                                                                                             <  3:306052/149‑1 (MQ=255)
                                aagTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGc                                                                                                            <  8:23817/149‑1 (MQ=255)
                                 agTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGc                                                                                                            >  5:97735/1‑148 (MQ=255)
                                    gCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCAt                                                                                                         >  3:272560/1‑148 (MQ=255)
                                    gCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATg                                                                                                        >  1:233084/1‑149 (MQ=255)
                                      tCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTGCTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGa                                                                                                       >  2:78577/1‑148 (MQ=255)
                                      tCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTGCTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGAt                                                                                                      >  2:82862/1‑149 (MQ=255)
                                             cccTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAggcgg                                                                                               <  2:105397/149‑1 (MQ=255)
                                             cccTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAggcgg                                                                                               >  6:67192/1‑149 (MQ=255)
                                               cTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCt                                                                                             <  4:280183/149‑1 (MQ=255)
                                                  aaTGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGa                                                                                           <  7:261802/148‑1 (MQ=255)
                                                  aaTGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAg                                                                                          >  7:116891/1‑149 (MQ=255)
                                                      ggTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCtt                                                                                       >  7:101200/1‑148 (MQ=255)
                                                      ggTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTg                                                                                      <  3:122363/149‑1 (MQ=255)
                                                        ttGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTGCTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGAc                                                                                    <  1:82862/149‑1 (MQ=255)
                                                        ttGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTGCTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGAc                                                                                    <  1:78577/149‑1 (MQ=255)
                                                        ttGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGAc                                                                                    <  3:60102/149‑1 (MQ=255)
                                                          gTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCg                                                                                  >  2:325913/1‑149 (MQ=255)
                                                             gATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAATACTCgg                                                                                                                              >  1:267480/1‑102 (MQ=255)
                                                             gATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCgg                                                                                                                              <  2:267480/102‑1 (MQ=255)
                                                              aTGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTc                                                                              <  1:325913/149‑1 (MQ=255)
                                                                 cgcAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGt                                                                            <  5:119041/148‑1 (MQ=255)
                                                                   cagcagCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCAc                                                                         <  1:116066/149‑1 (MQ=255)
                                                                                       cgccgcGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGaa                                                      <  8:201383/148‑1 (MQ=255)
                                                                                         ccgcGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGtt                                                   <  4:255883/149‑1 (MQ=255)
                                                                                                 tAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCaa                                           <  5:67192/149‑1 (MQ=255)
                                                                                                   gcgcTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAAt                                          <  5:75183/148‑1 (MQ=255)
                                                                                                         ccATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGCTGCGGGCAATGGTcagc                                   <  4:123839/149‑1 (MQ=255)
                                                                                                          cATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTcag                                    >  7:105138/1‑147 (MQ=255)
                                                                                                           aTGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGc                                                                                         <  1:232803/93‑1 (MQ=255)
                                                                                                           aTGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGc                                                                                         >  2:232803/1‑93 (MQ=255)
                                                                                                           aTGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTcagcag                                 >  7:23211/1‑149 (MQ=255)
                                                                                                                       aCGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCAc                     >  7:60584/1‑149 (MQ=255)
                                                                                                                        cgcgCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCAcc                    <  2:233084/149‑1 (MQ=255)
                                                                                                                             ttGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGt                 <  6:198271/147‑1 (MQ=255)
                                                                                                                             ttGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGt                 >  5:198271/1‑147 (MQ=255)
                                                                                                                               gggCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTAc               <  4:272560/147‑1 (MQ=255)
                                                                                                                               gggCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTAc               <  8:23211/147‑1 (MQ=255)
                                                                                                                                   ccGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTACCCCCTc         >  4:47273/1‑149 (MQ=255)
                                                                                                                                    cGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGGCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTACCCCCTCg        <  7:17446/149‑1 (MQ=255)
                                                                                                                                    cGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTACCCCCTCg        <  6:97735/149‑1 (MQ=255)
                                                                                                                                        cTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTACCCCCTCGGGAt    <  3:47273/149‑1 (MQ=255)
                                                                                                                                          cGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTACCCCCTCGGGATAc  <  4:181250/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GGTAATCTGGCGCCCCCTTGACCAGGCAGAAGAAGTGCTCCGCCGCCCTCAATGGGTTGTCGATGCGCAGCAGCCCGCGCTCATGTGCGCCGCGCAGTAGCGCTTCCATGCCTGCGAGCACGCGCTTGGGCCCGGCCTCGTAGAAGTACTCGCCAAAACTCGGGTCCAGGCTGCCCTGGGCCATGATCAGGCGGCTGAGCTTGACCGCTTCGTCACTGCTGATAAGCGCCTGGAAGTTGCGGGCAATGGTCAGCAGCACTTCTTCCACCGGTACCCCCTCGGGATAC  >  NC_002947/1719467‑1719753

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