Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 24,781 C→T G797G (GGC→GGT ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091324,7810CT94.1% 86.5 / ‑3.5 35G797G (GGC→GGTileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (14/18);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGTGTGGCGCGTCGTAGCTGCCAGACTGCGCTATATCACATCGCAGAAGCGCTGGTGCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAAC  >  NC_000913/24535‑25012
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
gATGCGCTTCTGCCCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACCTCGCAGAAGCGCTGGTTCGCTGGATGGAACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTgg                                                                                                                                                                                                                                      >  2:353166/1‑251 (MQ=255)
        tCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGt                                                                                                                                                                                                                              >  1:155458/1‑251 (MQ=255)
                gTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGcaat                                                                                                                                                                                                                      >  1:312345/1‑248 (MQ=255)
                    agaTGGGTTCCTTCTACCTCGACATCAGCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTATCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAAc                                                                                                                                                                                                                  <  1:234293/251‑1 (MQ=255)
                               ttCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATAAACGATGCGTTCTg                                                                                                                                                                                                       >  1:29320/1‑251 (MQ=255)
                                       cGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAg                                                                                                                                                                                                <  1:329928/250‑1 (MQ=255)
                                           atcagcAAAGCCCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGtt                                                                                                                                                                                           <  1:48118/251‑1 (MQ=255)
                                                      ccGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTg                                                                                                                                                                                <  1:70206/251‑1 (MQ=255)
                                                            gtgcacCGCCAAGGGGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCTCAGAATCGCTGGTTCGCGGGATGGCACCAAACCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAg                                                                                                                                                                          <  2:288619/248‑1 (MQ=255)
                                                             tACACAGCCAAAGCGGACAGCGTGGAGCGGCGTAGCTGCCAGCCTGCGCTGTCTCACATCGCAGAAGCGCTGGTACGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAGAAAACGTCATCACCGTGGAGTGGTACGCAGGCCGGTACGGTCAGGCAGCCAGTGAACCCAGGAACGATGCGTTATGGGACGGGCTGTTAAAAATGCGGGGCGAAgg                                                                                                                                                                         >  2:267587/1‑250 (MQ=255)
                                                                           ggACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTg                                                                                                                                                           <  1:238202/251‑1 (MQ=255)
                                                                               agcgtgGCGCGTCGTAGCTGCCAGAATGCGCTGTATCACATCGCAGAAGCGCTGGCGCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGGGGGGGGCTACCTGCCGGGCGAACGTGAAAATTACTTCTTCACCGGTGAGAGGTACGAAGGCCTGTTTGGGCTGGCAGACAGTGAAGCAATGAACTAAGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGTCGAAGAGAACAAAGACATAGAGc                                                                                                                                                        >  2:76316/4‑250 (MQ=255)
                                                                                  gtgGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTCCGCTGGCTGGCACCAATCCTCTACTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGCGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGCGCAAgc                                                                                                                                                   >  1:228072/1‑251 (MQ=255)
                                                                                        cgtcgtAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGc                                                                                                                                             >  1:252436/1‑251 (MQ=255)
                                                                                         gtcgtAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGCGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGAGCGTGcc                                                                                                                                            >  1:42673/1‑251 (MQ=255)
                                                                                            gtAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGGGCTGGCTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGAc                                                                                                                                         >  1:310186/1‑251 (MQ=255)
                                                                                                   ccAGACTGCGCTGTATCACATCGGAGAAGCGATGGTACGCAGGAAGGCACCAATCCTCTCCTTCACCGCTGATGAAGAGTGGGGCAAACCGCCAGGCGAAAGTGAAAAATACGTCTTCCCCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGTGCTGTTGAAAGTGCGGAGCGACATGAACAAAGGCATTGAGAAGACGCGTGGCGACAAGaa                                                                                                                                    >  2:226816/1‑249 (MQ=255)
                                                                                                     aGACTGAGCTGTATCACATTGTAGAAGCGCTGGGTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGTGCTACCTGCCGAGCGACCGTGACAAATTCGTCTTCACCGGTGAGTGATACGAAGGCCTTTGAGGTGTGGACGACAGTGACGCAATGTACGAGTCGTGCTGGGACGAGCTGATGCAAGTGCGTGGCGAAGCGAACAAAGTCATTGAGCAAGCGCGTGCCTATAAGaacgt                                                                                                                                 >  2:33408/1‑247 (MQ=255)
                                                                                                                   cacgtcggataaGCGCTGGGGCGCTGGATGGCACCAAGCGGCTCCTGCACCGCTGATGAAGTGTGGGGCTACCGGCCGGGTGAACGTGAAAAATACGTCTGCACCGGTGAGTGGTACGCAGGCCTGTTTGGCCTGGCAGTCAGTGAAGCGAGGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCAATGAGCAAGCTCGTGCCGACAATAAAGTGGGTGGCGCGCTGGa                                                                                                                  <  2:153608/241‑1 (MQ=255)
                                                                                                                   cacaTCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTgg                                                                                                                   <  1:171151/250‑1 (MQ=255)
                                                                                                                    acaccGCAGAAGCGCCGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAg                                                                                                                                                                                                                          <  2:282876/143‑1 (MQ=255)
                                                                                                                    acaTCGCAGAAGCGATGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAg                                                                                                                                                                                                                          >  1:282876/1‑147 (MQ=255)
                                                                                                                                                     cAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTg                                                                                <  1:230040/251‑1 (MQ=255)
                                                                                                                                                                  cACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGGTAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGGTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGa                                                                   <  2:312345/251‑1 (MQ=255)
                                                                                                                                                                         gatgaAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCt                                                            <  1:23790/251‑1 (MQ=255)
                                                                                                                                                                                   ggggCTACCTGCCGGGCGAACGTGAAAAATACGTCTGCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGa                                                   <  1:232520/250‑1 (MQ=255)
                                                                                                                                                                                             gCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTg                                        <  1:353166/251‑1 (MQ=255)
                                                                                                                                                                                                 ggCGAACGTGAAAAATACGTTTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGGAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCt                                    <  1:77037/251‑1 (MQ=255)
                                                                                                                                                                                                          gAAAAATACGTCTGCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCATGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTTGACCt                            <  1:88053/250‑1 (MQ=255)
                                                                                                                                                                                                                 aCGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTTGACCTCCGGCGCt                    <  1:230468/251‑1 (MQ=255)
                                                                                                                                                                                                                    tCTTCACCGGTGAGTAGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGAGATGGGCGATGAATGACGATTTTGCTTGTTGACATCCTGCGCtccc                 >  2:287914/1‑248 (MQ=255)
                                                                                                                                                                                                                         cccGGTGAGTGGGACGAAGGGCTTGTTGGTCTGGCAATCAGTGAAGCAATGAACTAGGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGGATAAAGTCATTGTGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAGCCTTGTATGCAGAACCGGAACTGGCGGCGAACCTGAGCTCCCTGGGCGATGAATTACGATTTGTCTTGTTGACCTCCGGCGCTACCGTTGc            <  2:311396/250‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATACGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTTGACCTCCGGCGCTACCGTTGc            >  1:287105/1‑248 (MQ=255)
                                                                                                                                                                                                                                   ggTACGAAGGCCTGGTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTTGACCTCCGGCGCTACCGTTGCAGACTATAac  <  1:337561/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGaa                                                  <  1:205991/199‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
GATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGTGTGGCGCGTCGTAGCTGCCAGACTGCGCTATATCACATCGCAGAAGCGCTGGTGCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAAC  >  NC_000913/24535‑25012

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