Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,738,637 T→A 11.6% V265D (GTC→GAC)  purR → transcriptional repressor, hypoxanthine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,738,6370TA11.6% 182.0 / 8.6 59V265D (GTC→GAC) purRtranscriptional repressor, hypoxanthine‑binding
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base A (4/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGT  >  NC_000913/1738492‑1738782
                                                                                                                                                 |                                                                                                                                                 
gtgCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTccc                                                                                                                                                <  1:182838/149‑1 (MQ=255)
      ggTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGa                                                                                                                                          <  2:255629/149‑1 (MQ=255)
       gcgACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCCATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGGGACGCGCatgc                                                                                                                                          >  3:331086/3‑145 (MQ=255)
       gTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCACTGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGAt                                                                                                                                         >  1:550671/1‑149 (MQ=255)
       gTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGGGACGCGCAggc                                                                                                                                          <  4:331086/148‑2 (MQ=255)
        tGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATg                                                                                                                                        <  2:236200/149‑1 (MQ=255)
            tttGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTTCGCGTCCCGCATGATGTTTc                                                                                                                                    >  1:674575/1‑149 (MQ=255)
            tttGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTc                                                                                                                                    >  3:539336/1‑149 (MQ=255)
                aaCCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTg                                                                                                                                >  1:242642/1‑149 (MQ=255)
                  ccTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCa                                                                                                                                             <  4:684700/134‑1 (MQ=255)
                  ccTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCa                                                                                                                                             >  3:684700/1‑134 (MQ=255)
                   cTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATc                                                                                                                             <  4:532981/149‑1 (MQ=255)
                      aaTCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGGGCACTTTGTGCTGCTGATGAAGTGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGt                                                                                                                          >  3:478444/1‑149 (MQ=255)
                            gTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGAt                                                                                                                    >  1:479096/1‑149 (MQ=255)
                               aTCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAAc                                                                                                                 <  1:472959/149‑1 (MQ=255)
                               aTCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAAc                                                                                                                 <  4:478444/149‑1 (MQ=255)
                                tCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACg                                                                                                                >  4:682566/1‑149 (MQ=255)
                                    gcCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTgcg                                                                                                            >  1:338916/1‑149 (MQ=255)
                                        tGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAc                                                                                                        <  2:40563/149‑1 (MQ=255)
                                        tGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAc                                                                                                        <  2:479096/149‑1 (MQ=255)
                                           agcaAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAcgcg                                                                                                     >  1:148175/1‑149 (MQ=255)
                                            gcaAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAcgcgc                                                                                                    <  3:682566/149‑1 (MQ=255)
                                            gcaAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAcgcgc                                                                                                    >  3:282563/1‑149 (MQ=255)
                                              aaaTCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACgcgcgc                                                                                                  <  2:453412/149‑1 (MQ=255)
                                                 tCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATc                                                                                                                             <  2:267769/119‑1 (MQ=255)
                                                 tCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATc                                                                                                                             >  1:267769/1‑119 (MQ=255)
                                                   cTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTAttt                                                                                             >  4:640939/1‑149 (MQ=255)
                                                      tCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTAc                                                                                          >  3:639610/1‑149 (MQ=255)
                                                        gCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGc                                                                                        <  4:131748/149‑1 (MQ=255)
                                                        gCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGc                                                                                        >  3:402076/1‑149 (MQ=255)
                                                         cAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGcc                                                                                       >  3:365637/1‑149 (MQ=255)
                                                          aGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCg                                                                                      >  1:35631/1‑149 (MQ=255)
                                                                   gCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAAAGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGAcc                                                                             >  4:468357/1‑149 (MQ=255)
                                                                       tACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGa                                                                         >  3:413201/1‑149 (MQ=255)
                                                                          tGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATcc                                                                      <  4:539336/149‑1 (MQ=255)
                                                                                 ttctGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGcc                                                               >  4:328274/1‑149 (MQ=255)
                                                                                 ttctGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGcc                                                               >  2:150489/1‑149 (MQ=255)
                                                                                            cGATATCATGGCAATGGGCGCTGTTTGTGCTGCTGATGCAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGTTCCATCAGCCAAAAGATTCGc                                                    <  2:674575/149‑1 (MQ=255)
                                                                                             gATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGGGACGCGCAggc                                                                                                                                          <  4:102334/62‑2 (MQ=255)
                                                                                             gATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGGGACGCGCAggc                                                                                                                                          >  3:102334/1‑61 (MQ=255)
                                                                                             gATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCt                                                   <  3:417777/149‑1 (MQ=255)
                                                                                               tatCATGGCAATGGGCGCTCTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCTAAAAGATTCGCTgg                                                 <  4:255395/149‑1 (MQ=255)
                                                                                                 tCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGt                                               >  3:240166/1‑149 (MQ=255)
                                                                                                   aTGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGa                                             >  1:653494/1‑149 (MQ=255)
                                                                                                     ggCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGaaa                                           <  2:148175/149‑1 (MQ=255)
                                                                                                        aaTGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAg                                        >  3:284271/1‑149 (MQ=255)
                                                                                                          tGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCg                                      <  2:338916/149‑1 (MQ=255)
                                                                                                           gggCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGt                                     <  1:150489/149‑1 (MQ=255)
                                                                                                                     tgtgCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCt                           <  3:113753/149‑1 (MQ=255)
                                                                                                                       tgctgcTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGt                         <  4:639610/149‑1 (MQ=255)
                                                                                                                            ctgATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGat                    <  3:328274/149‑1 (MQ=255)
                                                                                                                             tgatgaAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatc                   <  2:13978/149‑1 (MQ=255)
                                                                                                                               atgaAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatcgt                 >  3:196729/1‑149 (MQ=255)
                                                                                                                                 gaAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGAAAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCAATCAGCCAGAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatcgtat               >  4:563163/1‑149 (MQ=255)
                                                                                                                                    aTGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatcgtatcgt            >  2:164763/1‑149 (MQ=255)
                                                                                                                                       atccTGCGGGACGCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatcgtatc              >  2:481495/3‑144 (MQ=255)
                                                                                                                                       atccTGCGGGACGCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatcgtatc              <  1:481495/142‑1 (MQ=255)
                                                                                                                                       atccTGCGGGACGCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTcaa         >  3:468870/3‑149 (MQ=255)
                                                                                                                                              gcgTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGt  <  4:196729/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGT  >  NC_000913/1738492‑1738782

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: