Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,252,646	0	T	G	41.8%	86.9 / 12.7	56	V316G (GTA→GGA)	yeiI	putative kinase
Reads supporting (aligned to +/- strand):  ref base T (3/29);  new base G (23/0);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAAT  >  NC_000913/2252506‑2252792                                                                                                                                              |                                                                                                                                                     tgtTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCCGGG‑GTATAGCttt                                                                                                                                           >  2:1007871/1‑151 (MQ=255)     cAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCtt                                                                                                                                        <  1:34597/150‑1 (MQ=255)        cAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCtg                                                                                                                                        <  1:712578/147‑2 (MQ=255)         aaCTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAgg                                                                                                                                   <  1:799280/151‑1 (MQ=255)              tttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑ATATAACtttttt                                                                                                                                        >  2:1075779/1‑141 (MQ=255)                 gtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGt                                                                                                                            <  1:719451/150‑1 (MQ=255)                  tttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCAGGTAATTTAACTTTTTTGatgataa                                                                                                                                >  2:529699/1‑140 (MQ=255)                   ttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGG‑GTATAGCTTTTTTTAAGGAGACAGtttt                                                                                                                         >  2:841309/1‑150 (MQ=255)                   ttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGttt                                                                                                                          <  1:282871/150‑1 (MQ=255)                     aTTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGG‑GTATATTTTTTTTTAAAGGTAAAGttt                                                                                                                          >  1:459281/1‑148 (MQ=255)                     aTTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCcg                                                                                                                       <  1:532839/151‑1 (MQ=255)                       ttGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGTCCGGG‑ATATAGGttttttttaagggacagttt                                                                                                                          >  1:203536/1‑132 (MQ=255)                        tGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCcgcg                                                                                                                     <  2:80541/150‑1 (MQ=255)                          cccGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑GTATAGCTTTTTTGTAGGGAACAAttttccccgc                                                                                                                   >  1:1064144/1‑143 (MQ=255)                           ccGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAg                                                                                                                 <  2:52266/151‑1 (MQ=255)                            cGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCTGG‑‑ATATAGCttttttttagggtaaaagtt                                                                                                                          >  1:736119/1‑127 (MQ=255)                            cGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAgc                                                                                                                <  2:221401/151‑1 (MQ=255)                               tGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑ATTTAGGTTTTTTTTAGGATACAGTTTCCGCGACAGCGcc                                                                                                             >  2:630573/1‑151 (MQ=255)                                 aGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGG‑GTATAGCTTTTTTTTAGGGTAAAGATTTCGCGACAGGGccc                                                                                                            >  1:1105795/1‑150 (MQ=255)                                   ttaGTTTATTGCAGGGAAAAGGATGGAGAACAATTGTTGCTGATCGCGCTAGCGCATACGACAGTAGACAGTTTCTGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCcgcgttagc                                                                                                                <  1:936441/142‑6 (MQ=255)                                   tCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGG‑GTATATTTTTTTTTTAAGGTACAGGTTTCGCGACAGCGccc                                                                                                            >  1:182682/1‑148 (MQ=255)                                    cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGttt                                                                                                        <  1:309488/151‑1 (MQ=255)                                    cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGttt                                                                                                        <  1:2919/151‑1 (MQ=255)                                     aGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCCGGG‑ATATATCTTTTTTTTAGGGTACAGATTTCGCGACAGCGcccccttt                                                                                                       >  2:11906/1‑146 (MQ=255)                                            tGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGG‑ATATAAATTTCTTTTAGGATACAGTTTTCCCGGCAGCGCCCCGTTTgtg                                                                                                    >  1:264595/1‑145 (MQ=255)                                              cAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑GTATAGCTTTCTTGTAGGATACAGTTTTCGCGACAGCGCCCGTTTTTCGGGAGcc                                                                                              >  2:31519/1‑151 (MQ=255)                                              cAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑CTATTGCTTTTTTTAAGGATAAAGTTTTCCCGACA‑CGccccgcttttgggt                                                                                                 >  1:460305/1‑136 (MQ=255)                                                 cGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑ATATAACTTTTTTTTAGGATACAGGTTTCCCGGCCGCGCCCGGTTTTCGGGAGcc                                                                                              >  2:463148/1‑148 (MQ=255)                                                    aaaGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGATTTTTTTTAGGATAAAGTTTTCGCGACAGCGCCCCtttt                                                                                                       >  2:250142/1‑136 (MQ=255)                                                            gagaACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑ATATAGCTTTTTTTTAGGATAAAGTTTTCCCGACAGCGCCCCTTTTTCGGGAGCCCTCGCGGCAAtttt                                                                                >  2:910133/1‑150 (MQ=255)                                                               aaCAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTcgcg                                                                             <  2:731236/151‑1 (MQ=255)                                                                    tttttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑ATATATCTTTTTTTAAGGATACAATTTTCGCGACAGCGCCCGCTTTGCGGGTGCCTGCGCGGCAACTTCTCGCGcc                                                                         >  1:533513/1‑150 (MQ=255)                                                                     ttttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCGGGG‑ATAAAAGTTTTTTGAAGGAGACAGATTTCCCCCCACCGCCCCCTTTTCTGTTGGCTGCGCGGCAATTTTGCTCGcc                                                                         >  2:831910/1‑149 (MQ=255)                                                                              ccGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑Ctt                                                             <  1:988794/151‑1 (MQ=255)                                                                                gcgcCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGG‑ATGTAGATTTCTTTGAGGGGATAATTTTCCCGACAGCGCCCCCTTTTTTGGAGCGTGTGCGGCAATTTTGCGCGCCAGCGGCGG‑Ctt                                                             >  2:309310/1‑149 (MQ=255)                                                                                 cgcCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTaaa                                                          <  2:403675/151‑1 (MQ=255)                                                                                  gcCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAaac                                                         <  2:692014/151‑1 (MQ=255)                                                                                    cAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAaacaa                                                       <  2:1025084/151‑1 (MQ=255)                                                                                    cAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAaacaa                                                       <  2:415986/151‑1 (MQ=255)                                                                                        gcATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGTAG‑CTTAAACAaccct                                                   <  1:354736/151‑1 (MQ=255)                                                                                          aTACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGG‑ATATAGCTTTTTTTTAGGATAAAGTTTTCGCGGCAGCGCCCGTTTTTCGGTAGCCTGCGCGGGCATTTTGCGCGCCCGCGGCAG‑CTTAAAAAaaccccc                                                 >  2:671150/1‑148 (MQ=255)                                                                                             cGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGG‑ATATATAGTTCTTTTAGGGTACAGTTTTCGCGACAGCGCCCGGTTTGCGGGAGCCTGCGCGGCAATTTTGCGCGCCCGCCGCAG‑CCTAAACAaccccccccc                                              >  2:622889/1‑145 (MQ=255)                                                                                                   tAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGG‑ATATAGCTTTTTTTTAGGATAAAGTTTTCGCGGCAGCGCCCGTTTTTCGGTAGCCTGCGCGGCAATTTTGCGCGCCCGCGGCCG‑CCTAAAAAAACCCCCCCCGTCt                                          >  2:552870/1‑149 (MQ=255)                                                                                                    aGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGc                                        <  2:974976/150‑1 (MQ=255)                                                                                                        aGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGAt                                    <  1:993404/150‑1 (MQ=255)                                                                                                         gTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATaa                                  <  1:121511/151‑1 (MQ=255)                                                                                                          ttttGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATaa                                  <  1:588330/150‑1 (MQ=255)                                                                                                           tttGGTGCTGACGATGGTTTTATGGC‑GGGGCCGGG‑GTATAACTTTTTTTTAGGGTAAAGTTTTCGCGCCAGCGCCCGTTTTTCGGTAGCCTGCGCGGCGATTTCGCGCGCCCGCGGCCGCCTTAAAAAACCC‑ACCCCGTCTTCCGCTAAc                                 >  2:59495/1‑150 (MQ=255)                                                                                                                 gCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCtt                           <  2:915071/150‑1 (MQ=255)                                                                                                                    gACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCa                        <  1:525549/150‑1 (MQ=255)                                                                                                                     aCGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCAtt                      <  2:622731/151‑1 (MQ=255)                                                                                                                                aTGGC‑GGGCCCGGG‑ATATAACTTTTTTGAATGATAAAGTTTTCGCGACAGCGCCCGTTTTGCGGTAACCTGCGCGGCAATTTCGCGCGCCAGCGGGCG‑CTTAAACAACCCCACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATgg           >  1:830793/1‑151 (MQ=255)                                                                                                                                      gggCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTaa      <  1:823313/151‑1 (MQ=255)                                                                                                                                         ccTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACaa   <  1:342063/151‑1 (MQ=255)                                                                                                                                         ccTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGCAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACa    >  2:880653/1‑150 (MQ=255)                                                                                                                                            ggT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAAt  <  2:735769/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     TGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGT‑ATATAGCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAAT  >  NC_000913/2252506‑2252792

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