Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,453,971 G→T 9.1% S202I (AGC→ATC)  PP_1272 → multidrug MFS transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,453,9710GT 9.1% 171.6 / 6.4 55S202I (AGC→ATC) PP_1272multidrug MFS transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base T (3/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CCGAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGCC  >  NC_002947/1453830‑1454117
                                                                                                                                              |                                                                                                                                                  
ccGAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGAt                                                                                                                                             <  5:202217/149‑1 (MQ=255)
 cGAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATg                                                                                                                                            >  4:77639/1‑149 (MQ=255)
   aGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGc                                                                                                                                          <  8:266177/149‑1 (MQ=255)
       gggCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGc                                                                                                                                      <  3:176230/149‑1 (MQ=255)
        ggCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCt                                                                                                                                     >  2:136580/1‑149 (MQ=255)
           gcgGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTc                                                                                                                                    >  4:104701/1‑147 (MQ=255)
           gcgGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCaa                                                                                                                                  <  4:240908/149‑1 (MQ=255)
             gggTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGAGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGc                                                                                                                                <  8:133031/149‑1 (MQ=255)
               gTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGTGCCGATGGCCAGCTCAAGCGt                                                                                                                              >  6:1462/1‑149 (MQ=255)
                   cGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATGTCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGc                                                                                                                            <  1:218935/148‑1 (MQ=255)
                   cGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGc                                                                                                                            >  2:218935/1‑147 (MQ=255)
                          aaCGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAgg                                                                                                                   >  6:169692/1‑149 (MQ=255)
                           aCGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGc                                                                                                                  <  3:104701/149‑1 (MQ=255)
                             gCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCt                                                                                                                >  3:42814/1‑149 (MQ=255)
                             gCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCt                                                                                                                <  3:77639/149‑1 (MQ=255)
                              cccAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTg                                                                                                               >  1:293658/1‑149 (MQ=255)
                                  gcagGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGaaa                                                                                                           <  1:136580/149‑1 (MQ=255)
                                    agGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGc                                                                                                                  >  1:195546/1‑140 (MQ=255)
                                    agGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGc                                                                                                                  <  2:195546/140‑1 (MQ=255)
                                    agGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAg                                                                                                         <  2:64473/149‑1 (MQ=255)
                                      gcgcgcAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGTGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGc                                                                                                       <  5:1462/149‑1 (MQ=255)
                                      gcgcgcAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGc                                                                                                       >  4:85048/1‑149 (MQ=255)
                                        gcgcAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCATCGGCGCTGGCCactt                                                                                                                                     >  6:245261/1‑114 (MQ=255)
                                              cgcTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTc                                                                                                >  1:185893/1‑148 (MQ=255)
                                              cgcTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTc                                                                                                >  8:119562/1‑148 (MQ=255)
                                                cTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGAt                                                                                             <  2:185893/149‑1 (MQ=255)
                                                cTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGAt                                                                                             <  3:85048/149‑1 (MQ=255)
                                                   tCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCt                                                                                                                >  5:212017/1‑127 (MQ=255)
                                                   tCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCt                                                                                                                <  6:212017/127‑1 (MQ=255)
                                                      aTCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTcc                                                                                       <  7:164954/149‑1 (MQ=255)
                                                       tCAAGCCCGTGCTCGACTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCt                                                                                      <  7:281683/149‑1 (MQ=255)
                                                             ccGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACAccc                                                                                <  7:119562/149‑1 (MQ=255)
                                                               gTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGc                                                                              <  2:50641/149‑1 (MQ=255)
                                                                   tCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTac                                                                                    <  4:282029/139‑1 (MQ=255)
                                                                   tCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTac                                                                                    >  3:282029/1‑139 (MQ=255)
                                                                     gCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATcac                                                                         <  5:169692/148‑1 (MQ=255)
                                                                       cTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAgc                                                                      >  4:198909/1‑149 (MQ=255)
                                                                             aaCGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGt                                                                >  7:133548/1‑149 (MQ=255)
                                                                                gCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGa                                                             >  2:32468/1‑149 (MQ=255)
                                                                                                   cggcgCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAAcgc                                          >  2:44469/1‑149 (MQ=255)
                                                                                                                ggTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTgg                             <  6:183972/149‑1 (MQ=255)
                                                                                                                            tCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGt                   >  7:120615/1‑148 (MQ=255)
                                                                                                                            tCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGt                   >  8:123521/1‑148 (MQ=255)
                                                                                                                                   tCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg           >  6:35763/1‑149 (MQ=255)
                                                                                                                                     aaGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGgc         <  1:44469/149‑1 (MQ=255)
                                                                                                                                     aaGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGAGGACGGTATGGTCGGCGAAc                                            <  4:129303/114‑1 (MQ=255)
                                                                                                                                     aaGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGAGGACGGTATGGTCGGCGAAc                                            >  3:129303/1‑114 (MQ=255)
                                                                                                                                      gctGGCCATCGGCGCTGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg           >  2:109833/3‑146 (MQ=255)
                                                                                                                                      gctGGCCATCGGCGCTGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg           >  4:255252/3‑146 (MQ=255)
                                                                                                                                      gctGGCCATCGGCGCTGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg           <  1:109833/144‑1 (MQ=255)
                                                                                                                                      gctGGCCATCGGCGCTGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg           <  3:255252/144‑1 (MQ=255)
                                                                                                                                          gCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTc     <  3:198909/148‑1 (MQ=255)
                                                                                                                                          gCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTc     <  2:293658/148‑1 (MQ=255)
                                                                                                                                           ccAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGc   >  4:206119/1‑149 (MQ=255)
                                                                                                                                            cAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGcc  >  2:16811/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CCGAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACAT‑CCTCACGGCTCAAGTGGCCAGCGCCGATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGCC  >  NC_002947/1453830‑1454117

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