Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,223	T→C	N207N (AAT→AAC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,223	0	T	C	90.9%	81.3 / ‑0.6	35	N207N (AAT→AAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (16/14);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

GCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAA‑GAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAA  >  NC_000913/4619986‑4620462                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                gctgGATAAACTGGTCGAACGCTCTAATCTGCCGGTTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTa                                                                                                                                                                                                                                    <  1:179060/251‑1 (MQ=255)    ggATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAA‑GAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATcc                                                                                                                                                                                                                                 >  1:132260/1‑250 (MQ=255)     gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCCGGCGTGTTATCGc                                                                                                                                                                                                                                >  1:101726/1‑251 (MQ=255)        aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCCGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCCCGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGACGGCTACACAATCGGTCGGGTTATTGCTCg                                                                                                                                                                                                                             >  2:76168/1‑251 (MQ=255)              gTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTGCAACATCGGtagggtttcctcccggtcgtt                                                                                                                                                                                                                       >  2:143519/1‑231 (MQ=255)                 aaccgcctaaattgtcggggttctcTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGAACAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAA‑GAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTAt                                                                                                                                                                                                                    <  2:132260/227‑1 (MQ=255)                               gCCGGTTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCg                                                                                                                                                                                                      <  1:143519/251‑1 (MQ=255)                                         cTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTcc                                                                                                                                                                                             >  1:314005/1‑250 (MQ=255)                                                                   ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTgg                                                                                                                                                                  >  1:92808/1‑251 (MQ=255)                                                                   ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTAGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCCGCGTCCCGGGAACCGTCACGACCTgg                                                                                                                                                                  >  1:285591/1‑251 (MQ=255)                                                                    gTCATTCTGGAGCAGCTGTGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGACCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGCAGAGCTGACCAACGGCGGCTACACCATCGGTCGTGTTATCGCTAGGCCGTTTATCGGTGACCAAGCACGTAACTTCCAGCGTAACGGTAACCGTCaagcaccggc                                                                                                                                                                 >  2:67843/1‑242 (MQ=255)                                                                                         gAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCaccgac                                                                                                                                            <  1:333415/251‑1 (MQ=255)                                                                                                acaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTTAGCCGCCAGCCCCGACCGtgctg                                                                                                                                     >  1:350041/1‑251 (MQ=255)                                                                                                       aaCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCACAACATCGGTCGTGGTATCGCCCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTCACGGTAACCGTCACGACCTGGCTGTTGGGCCACCAGCACCAGCCGGGCTGCaaaacc                                                                                                                              >  1:98013/1‑246 (MQ=255)                                                                                                          cGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGATCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                                                                                                           <  1:350669/251‑1 (MQ=255)                                                                                                           ggCAAGCCTATTTTCTATCCCTCCGCTGTCTCCGTGTTCCAGATTGCCTGCCACGAATGATCCTTT‑GGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACATCACCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCTGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGt                                                                                                                          <  1:256993/251‑1 (MQ=255)                                                                                                                      tttCTATACCTCCGCTGCCTCCGTGTACCAGATTGCCTGCCACGAC‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCCGAAACTGGTTGATGAAAAac                                                                                                               >  1:303974/1‑251 (MQ=255)                                                                                                                       ttCTATACCTCCGCTGACTCCGTGTTCAAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAac                                                                                                               <  1:151379/250‑1 (MQ=255)                                                                                                                           atcccTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCGCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAacac                                                                                                             <  1:360013/245‑1 (MQ=255)                                                                                                                                   gCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtggt                                                                                                  <  1:158973/251‑1 (MQ=255)                                                                                                                                         tCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGAGGCTACAAAAACGGTCGTGTTATCGCTCGTCAGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGGGCTGCAGAAACTGGTTGATGTAAAACACGGTCAGGGGGTTTATGt                                                                                            >  1:22104/1‑251 (MQ=255)                                                                                                                                          ccGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                                                                                           >  1:306294/1‑251 (MQ=255)                                                                                                                                          ccGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                                                                                           <  1:319598/251‑1 (MQ=255)                                                                                                                                           cGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCg                                                                                          >  1:313362/1‑251 (MQ=255)                                                                                                                                                    gATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCAGCTACAACATCTGTCGCGTTATCGCACCTCCGTTTATCGGTGACCAAGCCGGTGACTTCCAGCGTACCGGGAACCGTCACGCCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGGAACTGGTTGATGCAAAACACGGTTAGGTGGGTTCTTGCGGTAAAAttt                                                                                 >  1:132679/1‑250 (MQ=255)                                                                                                                                                                   aaggaaccttaggccaggATAAACCCTACGAACTGTGCGAAATCGCCCGTGAAGAGCGGACCACCGGCGGCGACAACATCGGTCGGGTTATCGCTCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGcc                                                                   <  2:100362/235‑1 (MQ=255)                                                                                                                                                                              gtgccggATAAACTCTACGAACTGCTCGAAATCGCTCTTTGAGTGCTGACCCACGGCGGCTACAATTTCGGTCGTGTTATCGCTCGTCCGTTTATGGGCGACAAAGCCGGTAACTCCCATCGTACCGGTAACCGTCATGACCTGGCTGTTGAGCCGCGAGCACCGACCGTGCGGCAGAAGCTGGTTGATGAAAAACACGGGCAGGTGGTTTCTGTCGGGAAAATTGCGGACATCTACGCCAACTGCGGTAt                                                        <  2:165880/246‑1 (MQ=255)                                                                                                                                                                               gccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTCAAATTGCAGACATCTACGCCAACTGCGGTATc                                                       >  1:311224/3‑251 (MQ=255)                                                                                                                                                                                   ggATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCAcc                                                    >  1:339207/1‑250 (MQ=255)                                                                                                                                                                                      tAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGCCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAATTTCCGACATCTACGCCCACTGCGGTATCAccc                                                   >  1:159809/1‑247 (MQ=255)                                                                                                                                                                                        aaccctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaagaaa                                              <  1:374797/247‑1 (MQ=255)                                                                                                                                                                                         aCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAg                                             <  1:161927/251‑1 (MQ=255)                                                                                                                                                                                              aCGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGaaa                                        <  1:76168/251‑1 (MQ=255)                                                                                                                                                                                                                                cAACGGCGGCTACAACATCTGTCGTGTTATCGCTCGTCCGTTTTTCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGTCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGtcacca      <  1:69963/251‑1 (MQ=255)                                                                                                                                                                                                                                       ggcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCaa  >  1:53314/1‑248 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                GCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAA‑GAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAA  >  NC_000913/4619986‑4620462

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: ‑