Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,414,829	A→G	E28G (GAA→GGA)	pta →	phosphate acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,414,829	0	A	G	97.1%	96.8 / ‑3.9	34	E28G (GAA→GGA)	pta	phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/10);  total (24/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

AAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCT  >  NC_000913/2414589‑2415075                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       aaaGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGc                                                                                                                                                                                                                                              <  2:59884/251‑1 (MQ=255)                ccTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGCAAGAGGATAAACCGGTTCCCCGATTTTTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGAACAGCGT‑TAGCCTTGG‑CGTGATCCGTGCAATGGGCCGCCAAGGCGTTCGTCTGAGCGtt                                                                                                                                                                                                                               >  2:208450/1‑251 (MQ=255)                        ggtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGGATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAAcc                                                                                                                                                                                                                       >  1:142552/1‑251 (MQ=255)                                    caacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGcc                                                                                                                                                                                                           >  1:40902/1‑251 (MQ=255)                                        gaagaaTTGGTTATCACGAAAGACGCGAGCAGCATGACAGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGAAGCTGTTTTGTAACCAGCCAAAACGGCGGGAACGCAAAAGGATACACCGAGTCCCGAGTTATTATGATGATACCTACCGGCACCAGCCTCGGGTAAACCAGCATACA‑CCTTGGTCG‑GAGGAGCGAAATGGGACGCAAAGGGGGTCGTCTGAGCGTTTTCCAAACTATCGCTCAGCcgcgg                                                                                                                                                                                                       >  2:225206/1‑250 (MQ=255)                                             aCTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCgg                                                                                                                                                                                                  >  1:202393/1‑251 (MQ=255)                                                ggttaccgggcaggccgcgACCCGCCTGACCGCCT‑ATGTTCACTCCCCCAGCTCAGCTGGCGGTTCTGTTTTGTATCCCGCCAATGCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTg                                                                                                                                                                                                <  2:189084/235‑1 (MQ=255)                                                      cgcgCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATgg                                                                                                                                                                                         >  1:179365/1‑250 (MQ=255)                                                          aaaGACACGAGCCGCCTTACTGGCTGAT‑TTCACACCGCCAGCGCAGCGGGCGGTGCGGGTTTGTAACCCGCCAAGTCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCTAAGACGTTGTTCTGAGCGTTTTCAAACGTATCGCTCAGCCGCGTACCGGTGGCGATGCGccc                                                                                                                                                                                     <  2:112764/250‑1 (MQ=255)                                                           aaGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCg                                                                                                                                                                                    >  1:50018/1‑251 (MQ=255)                                                                            aCTGCCTGAT‑TTCACCCCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATc                                                                                                                                                                   <  2:141434/251‑1 (MQ=255)                                                                                     atgtccaacCGCCAGCTCAGCTGGCGGTGCTGGTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATcgtgcgtg                                                                                                                                                            <  2:222125/243‑1 (MQ=255)                                                                                                     tcagctGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATTGGACGCAAAGGCGTCCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCacg                                                                                                                                           >  1:25769/1‑251 (MQ=255)                                                                                                       agctGGCGGTGATGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGGGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCacgac                                                                                                                                         <  2:3267/251‑1 (MQ=255)                                                                                                          tGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGGGTCCCGTTTTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGGCGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCCCGTACCGGGGGCGATGCGCCCGATCAGACTACGACAATCGTGCGTGCGAACTCTTCCACCACGACGGc                                                                                                                                      >  1:161257/1‑251 (MQ=255)                                                                                                           ggCGGTGTTGTTTTGTAACCTGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGcc                                                                                                                                     <  1:105072/251‑1 (MQ=255)                                                                                                              ggTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGccgcc                                                                                                                                  >  2:194666/1‑250 (MQ=255)                                                                                                                tGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGccgctg                                                                                                                                 >  1:212550/1‑250 (MQ=255)                                                                                                                                      tCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGtt                                                                                                          >  1:16095/1‑251 (MQ=255)                                                                                                                                      tCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCAGATCCATACCGGAACCAGCGTCGGTCTGACCAGAGT‑CAGCCTTGG‑CGTGCTCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGCTCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGtt                                                                                                          >  1:167950/1‑251 (MQ=255)                                                                                                                                          cggTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGa                                                                                                        >  2:187359/1‑249 (MQ=255)                                                                                                                                                    agagGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCGGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTAGCGTGCGTGCGAACTCTTCCACCACTACGGCGGCTTAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTc                                                                                            <  1:15796/251‑1 (MQ=255)                                                                                                                                                      agGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGACGGTCTGCTTTcc                                                                                           >  2:203095/1‑250 (MQ=255)                                                                                                                                                        gATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGc                                                                                        >  2:70748/1‑251 (MQ=255)                                                                                                                                                        gATAAACCCAGTCCCGTATTAATATGCTGTTCCCTACCAGCACCAGCGGCGGCCCGACCAGCGT‑CACCCTTGG‑CGTCCTCCGTGCCATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCCGGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGGGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGCAGGTCTGCTTTCCAGc                                                                                        >  2:112210/1‑251 (MQ=255)                                                                                                                                                                  gtCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGaa                                                                                >  2:10006/1‑249 (MQ=255)                                                                                                                                                                  gtCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTATGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGCCGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGaa                                                                                >  1:228571/1‑249 (MQ=255)                                                                                                                                                                    cccGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CCGCCTTGG‑CGTTATCATTGCAATGGGACTCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGATACGTTGAAGGGCTGCTTTCCAGCAATCAGaa                                                                                >  2:31629/1‑247 (MQ=255)                                                                                                                                                                                         ctcccGGAAACAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGCCGCAAAGGCGCTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCGGATCAGACTCCGACTATCGCGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGCAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGAtcgtc                                                       >  2:3940/4‑251 (MQ=255)                                                                                                                                                                                            ccGGAACTAGCGTCGGTATGACCAGCGT‑CAGCCCTCG‑CGGGATCCGTGAAATGCGGCGCAAAGGCGTTCGTCTGAACGTTTTCAAACCTATCGCTCCGCCGAGTACCGGTGGCGAAGCGCCCGATCAGACTACGACTATCGAGCGTGCGAACTCTTCAACCACGACGGACGCTGAACCGCTCACAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCACAAAGATGTGCTGCTGGAAGAGATCGTCGc                                                     >  2:124723/1‑250 (MQ=255)                                                                                                                                                                                                                      gT‑CAGCCTTGG‑CGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGc                          <  2:140677/251‑1 (MQ=255)                                                                                                                                                                                                                                    gTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCt              >  1:202084/1‑251 (MQ=255)                                                                                                                                                                                                                                      gATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTgg            <  1:70748/251‑1 (MQ=255)                                                                                                                                                                                                                                                aTGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCt  >  1:87603/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       AAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAT‑TTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGT‑CAGCCTTGG‑CGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCT  >  NC_000913/2414589‑2415075

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