Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,915 C→T A71T (GCG→ACG)  yjcO ← Sel1 repeat‑containing protein YjcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,9150CT97.0% 106.7 / NA 34A71T (GCG→ACG) yjcOSel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (16/16);  minor base G (1/0);  total (18/16)
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.98e-01

AACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGCGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGT  >  NC_000913/4296678‑4297144
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
aaCCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGgcatcagc                                                                                                                                                                                                                           >  1:375234/1‑250 (MQ=255)
  ccAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGgcatcagcatc                                                                                                                                                                                                                        >  1:388214/1‑251 (MQ=255)
         gcaTTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCtt                                                                                                                                                                                                                 <  2:9071/251‑2 (MQ=255)
           aTTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGCGGCGTTTTCTAAAAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACAAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGCACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTc                                                                                                                                                                                                               >  1:21499/1‑251 (MQ=255)
                 gCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCggcggc                                                                                                                                                                                                         >  1:42621/1‑251 (MQ=255)
                                     ctcCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCgg                                                                                                                                                                                     >  1:391218/1‑251 (MQ=255)
                                                           tCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGGCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGc                                                                                                                                                                <  2:418730/250‑1 (MQ=255)
                                                                                                 cGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGt                                                                                                                         >  1:216735/1‑251 (MQ=255)
                                                                                                     ggAGTATTTAACAGAATATGCGCCAGCGTTACTTCCCCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGAGAGTCCA‑GACTGACCTGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCGGTGCCCAATATTCGGCGTTGCTCAAAGCGGCGGAGCTAAACCAGCTGGCGGCAAGAAAAGATTGTGCGCTTCGTTcccc                                                                                                                     <  2:157284/250‑1 (MQ=255)
                                                                                                      gAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTcccc                                                                                                                     <  1:6792/250‑1 (MQ=255)
                                                                                                             tACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAGCTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGCTCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGATGTCGGCGGTTTCCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCctct                                                                                                              <  1:487783/250‑1 (MQ=255)
                                                                                                               ccAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GGCTGACCGGATTGGTGAGTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCAGCGCCGCTGGCGGCGGCTGTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTctgct                                                                                                           <  1:465713/251‑1 (MQ=255)
                                                                                                                    atatGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTcccc                                                                                                                     >  1:248159/1‑236 (MQ=255)
                                                                                                                    atatGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTcccc                                                                                                                     <  2:248159/236‑1 (MQ=255)
                                                                                                                          gcCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAAt                                                                                                 >  2:101823/1‑250 (MQ=255)
                                                                                                                                  tACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTg                                                                                         <  1:89050/250‑1 (MQ=255)
                                                                                                                                        aCCTTATTTACTGCCCGCTTGCGCCGCTTTCTATGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCAAGCAACGCGCAGGCATCAGAATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAACTGTCGGCAAGAAAATATTGTGCGCGTAGGTCCCCGGCCTCTGCTGCATTTCAATACTGGCTGCCAggtc                                                                                  >  1:141691/1‑249 (MQ=255)
                                                                                                                                                tactactcGCTTGAGCCGTTGTCTTTGCAAGAACTTTGTCTTGTGGATAGTCCATGA‑TGACCGGATTTGTGATTTTGATCTGCGCCAGCAGCGTTCAGGCATCAGCATCACCGCTGTCGGCGGCTTGCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGCGCTAAACCAGCTTGCGGCAAGAAAATATTGTGCGCGTGGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCATGCTCGGAGTc                                                                           <  2:346834/243‑1 (MQ=255)
                                                                                                                                                        gCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCGTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGGCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTg                                                                  <  2:388214/251‑1 (MQ=255)
                                                                                                                                                         cTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGCTTGGCGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAGTCGCCGGCGCTAAACCAGCTGTCGGCAAGAACATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGt                                                                   >  2:221124/1‑249 (MQ=255)
                                                                                                                                                         cTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAGTCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGGCGGGCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGt                                                                   <  1:221124/249‑1 (MQ=255)
                                                                                                                                                           tGCGCAGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTACACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGc                                                               >  1:331249/1‑251 (MQ=255)
                                                                                                                                                                cgcTTTCTCTGCAAGACCTTTTGCTTGTGGATAGTCCA‑GACTGACCGGAGTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAATCCAGCTGGCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGa                                                          <  2:155204/251‑1 (MQ=255)
                                                                                                                                                                 gcTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTTATCTGCGCCAGCAGCGTGCAGGCATCAGCCTCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTTCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCTCGTCGGTCCCCGGCCTCTGCTGCCTTTCAATACTGGCTGCCAGGCTCGGAGTCGTTTGCGTGGGCAAAGa                                                          >  2:250658/1‑250 (MQ=255)
                                                                                                                                                                       tctGCAAGAGCCTTTACTTGTGGCTAGACCA‑GACTGACCCGATTGGACATTTTGATCTGCGCCCGCACAGGGCCGGCATCCACCTCAACGCTGTCGGTGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACAAGCTGTCGGAAAGAACAAATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGGCAGACTAGGAGTCGTTGGCGTGGGGAAAGAGTGTCa                                                    >  2:481699/1‑250 (MQ=255)
                                                                                                                                                                            aaGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCACGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAaac                                               >  1:422695/1‑250 (MQ=255)
                                                                                                                                                                                          gtgGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGAt                                 >  1:343156/1‑250 (MQ=255)
                                                                                                                                                                                          gtgGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATa                                <  2:343156/251‑1 (MQ=255)
                                                                                                                                                                                           tgGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATa                                >  1:478620/1‑250 (MQ=255)
                                                                                                                                                                                                     a‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTGCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTTGTTTGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCAt                     <  1:283982/251‑1 (MQ=255)
                                                                                                                                                                                                     a‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCAt                     <  2:112830/251‑1 (MQ=255)
                                                                                                                                                                                                                 aTTGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCCTCAGCATCACCGCTGACGGCGGCTTTCTGTGCACAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGCAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTAAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATAc          >  1:483236/1‑251 (MQ=255)
                                                                                                                                                                                                                   tGGTGATTTTGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACg         <  1:94274/250‑1 (MQ=255)
                                                                                                                                                                                                                           ttGATCTGCGCCAGCAGCGTGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGt  >  2:455460/1‑249 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
AACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGAATATGCGCCAGCGTTACTTCACCTTCTTTACTGCCCGCTTGCGCCGCTTTCTCTGCAAGAACTTTTGCTTGTGGATAGTCCA‑GACTGACCGGATTGGTGATTTTGATCTGCGCCAGCAGCGCGCAGGCATCAGCATCACCGCTGTCGGCGGCTTTCTGTGCCCAATATTCGGCTTTGCTCAAATCGCCGGAGCTAAACCAGCTGTCGGCAAGAAAATATTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGT  >  NC_000913/4296678‑4297144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑