Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,648 C→T V140I (GTT→ATT)  yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,6480CT97.1% 104.6 / NA 35V140I (GTT→ATT) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/16);  minor base G (0/1);  total (18/17)
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.84e-01

CACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAA  >  NC_000913/4616403‑4616885
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
cACCGAACAGTTGGCGATGAGAGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCCTCAGGGGTGATTTCAGGGTCGTAGATCAAGGTCAGGTTCGGGGAAACCTGATTAAACTCTGAGTCTGACCGTAAGATCGCGCGGGGAATGGGCGAATAAGACGGGCCGATATTGGCGTGCATCAAGGCGTCTGGCAGGGTGCTGTCGAGGTAACGCCAAAAACGTTTTATTCGAATATCg                                                                                                                                                                                                                                           >  2:185897/1‑250 (MQ=255)
      aCCGTTGGCCATGAGCGGTATGCTACATTTACCGATGTTCTTCGACACTTCCGGCAGCAGGCCATAAGGGGTGCTTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTACCGCCAGAAACGTTTTATTCGAATATCGATCTCGt                                                                                                                                                                                                                                    >  2:155431/1‑251 (MQ=255)
        cGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCt                                                                                                                                                                                                                                  >  1:18141/1‑251 (MQ=255)
                aTGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGtt                                                                                                                                                                                                                            >  1:85858/1‑249 (MQ=255)
                    gCGGTTTGCTACATGCACGGAGGTTCTTCGCCACTTCCAGCAGCAGTGCATCAGGGTTGATTTCAGGATCGTGGATAAAGGTCAGGTTCGGGGAAACCTGCGTCATCTCTGCGTCTGCACGTAAGATCTCGCGGGTAAGGGGCGAATCAGGCGTTCCGGTAGTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAAt                                                                                                                                                                                                                       <  1:139449/250‑1 (MQ=255)
                        tttGCTACATTCACAGATGTTATTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTAGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATAAGACGGGCCGATATTGGAGAGCATAAAGGCGTCTGGCAGGGTTCTGTTGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTa                                                                                                                                                                                                                   >  2:91261/1‑250 (MQ=255)
                             tccattcccggATGTTCTTAGCCACTTCCAGCAGCCGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGCAACATGCTTCAACTCTGCGTCTGACCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGCGGTAACGCCCGAAACGTTTTATTCGAATATGGATCTCGTATTGTGTTAGAATTCTAACATat                                                                                                                                                                                                             >  2:20057/11‑250 (MQ=255)
                                   cacaGATGTTCTTCGCCACTTCCCGCAGAAGGTCATCAGGGGTGATTTCAGGCTCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGc                                                                                                                                                                                                       >  2:233310/1‑251 (MQ=255)
                                                    aCTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTg                                                                                                                                                                                      <  2:260820/251‑1 (MQ=255)
                                                                 ggTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                                                                                                                                                                         >  1:85516/1‑251 (MQ=255)
                                                                           ggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg                                                                                                                                                               <  2:85858/251‑1 (MQ=255)
                                                                            gTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGAGCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                                                                                                                                                              <  1:253880/251‑1 (MQ=255)
                                                                            gTGATTTCAGGATCGTAGATAAAGGGCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCGGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGGCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTGGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                                                                                                                                                              <  1:242243/251‑1 (MQ=255)
                                                                               aTTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                                                                                                                                                           >  1:156406/1‑251 (MQ=255)
                                                                                  tCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGATCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGTTCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTCGACCGGCATCGATGTGACCg                                                                                                                                                         >  1:247264/1‑250 (MQ=255)
                                                                                               tAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg                                                                                                                                           >  1:95073/1‑251 (MQ=255)
                                                                                                aaaGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg                                                                                                                                           <  2:85516/250‑1 (MQ=255)
                                                                                                       aGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt                                                                                                                                   <  1:104930/251‑1 (MQ=255)
                                                                                                         gTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt                                                                                                                                   >  1:16105/1‑249 (MQ=255)
                                                                                                              gggAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                                                                                                                                              <  1:82273/199‑1 (MQ=255)
                                                                                                              gggAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                                                                                                                                              >  2:82273/3‑201 (MQ=255)
                                                                                                                       gCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATAGTGGCGTGCATAAAGGCGTCTGGCATGGTTCTGTCGAGGTAACGCCAGAAACGTTTTAGTCGAAGATCGATCGCGTCTTGTGTTAGAATTCTAACATATGGTGGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGt                                                                                                                   <  2:105625/251‑1 (MQ=255)
                                                                                                                                   cgtctgcaCGTAAGATCGCGCGGGTAAGGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCTGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGGAGAGAGAGTGCGTTATCAAGa                                                                                                       <  2:108554/249‑1 (MQ=255)
                                                                                                                                        gcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCt                                                                                                    <  2:18141/249‑1 (MQ=255)
                                                                                                                                                      gcgGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATGCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTtccagct                                                                                    <  2:245602/251‑1 (MQ=255)
                                                                                                                                                                    aaTCAAACCGGCCGATATTGGGGTGAATAAAGGCGACAGCAAGGGCTCTGTCGAGGTAACGACAAAAAAGTTTTATTCGACAATCGATATCGTCTTGCGTTAGCCTTACAACATATGGTTGCAACAACGCAACCAGTTGCACCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCATTCCCGCGCCAGCCATTCggc                                                                       >  2:141661/1‑249 (MQ=255)
                                                                                                                                                                     aTCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAAc                                                                     >  2:1185/1‑251 (MQ=255)
                                                                                                                                                                          aCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTc                                                                >  1:264378/1‑251 (MQ=255)
                                                                                                                                                                                    atTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAAc                                                      <  2:95073/251‑1 (MQ=255)
                                                                                                                                                                                        gCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGc                                                  <  1:1185/251‑1 (MQ=255)
                                                                                                                                                                                            gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTa                                               >  1:237482/1‑250 (MQ=255)
                                                                                                                                                                                                    gCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTa                                       >  1:45030/1‑250 (MQ=255)
                                                                                                                                                                                                                     gtctaGGTAACGCCAGCATCGTTTTATTCGAATATCGATCTCGTCTTGTGGTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGTCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGATCCGTTCGCCAGAAACCGGGCGTAATCGGGTAAGACATAGCGCGGtt                        <  1:45009/244‑1 (MQ=255)
                                                                                                                                                                                                                       cGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATAGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTa                    <  2:16105/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cGTTTTATTCGAATTTCGATCTCGTCTTTTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAATAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCTGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCTCGGTTTGTACGGCGCATGACCTTCaaa  <  1:131828/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
CACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAA  >  NC_000913/4616403‑4616885

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