Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,234,479	C→T	A241V (GCG→GTG)	pgi →	glucose‑6‑phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,234,479	0	C	T	95.7%	136.3 / ‑6.4	46	A241V (GCG→GTG)	pgi	glucose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/20);  minor base A (2/0);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTCCACCACGCCTGCCGAGAAAAACCTGCCTGTACTGCTGGCGCTGA  >  NC_000913/4234237‑4234715                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             cGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTAGGACCAACGCTCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACGCGGTGCAAAACGCTTTTTGGCGCt                                                                                                                                                                                                                                        <  2:240114/249‑1 (MQ=255)     ccATACATGGTGACCGAAGCTCTGCGTCCGTACAAACACCACCTAAACATGCCCTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCCGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTccac                                                                                                                                                                                                                                  >  2:41441/1‑251 (MQ=255)         aCATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAAt                                                                                                                                                                                                                              <  1:43554/251‑1 (MQ=255)             ggTGACCGAAGCTCTGCGTCCGTTCAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGcc                                                                                                                                                                                                                           <  1:373708/250‑1 (MQ=255)                            gCGTCCGTACAAAAACCACCTGAACATGCACTTTTTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCTCCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGTTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCGACCAATGCCAAAGCCGTTGGCGAg                                                                                                                                                                                                            <  2:228715/250‑1 (MQ=255)                            gCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGGCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCCCCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTGGCAAAACACTTTGTGGCGCTGTCCACCAATGCCATAGCCGTTGGCg                                                                                                                                                                                                              <  2:71143/248‑1 (MQ=255)                             cGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGAGATGAAAAACACGTTGCAAAACACATTGAGGCGCTTTGCCCCAATGCCAAAGCCGTTGGCGAGtt                                                                                                                                                                                                          >  2:259922/1‑251 (MQ=255)                             cGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTAACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGATTCCTGAAAGCGGCAGGTGATGACAAACACGTTGCAAAACACTTTGTTGCGCTTTCCACCAATGCCCCAGCCGTTGGAGAGtt                                                                                                                                                                                                          >  2:19202/1‑251 (MQ=255)                                ccGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTTATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCTTTGGCGAGTTTgg                                                                                                                                                                                                       <  1:367788/251‑1 (MQ=255)                                ccGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTg                                                                                                                                                                                                        <  2:384996/250‑1 (MQ=255)                                                                ccgatcggtagATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCCCGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACAAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACAAGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTgg                                                                                                                                                                        <  1:328686/240‑1 (MQ=255)                                                                    aaCGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTg                                                                                                                                                                   >  1:110244/1‑251 (MQ=255)                                                                       gTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGCAACCACGCTGGTCTTGATAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCATGACTGATTCCTGAAAGCGGCAGGAGATGAAAAACACGTTGCAAAACACTTTGTTGCGCTTTCCACCAATGCCAAAGCCGTTGGCGCGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTggagatcggg                                                                                                                                                                >  2:328686/1‑243 (MQ=255)                                                                       gTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGt                                                                                                                                                                >  1:330734/1‑251 (MQ=255)                                                                               gACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGcc                                                                                                                                                        >  1:219192/1‑251 (MQ=255)                                                                                   cacaTCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTa                                                                                                                                                    <  1:102288/251‑1 (MQ=255)                                                                                           ggAAGTGCTTAAAAAAGTAAACCCGGAACCCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACGTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTtgt                                                                                                                                            <  2:320676/251‑1 (MQ=255)                                                                                           ggAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTtgg                                                                                                                                            >  2:358667/1‑250 (MQ=255)                                                                                           ggAAGTGCTGAAAAAAGTAAAACCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGAGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGACGTTGGCGAGTTTGGTATTGATACTGCCAACATGGTCGAGTTCTGGGACTGGGTTGGGGGCCGTTACTCTTtgg                                                                                                                                            >  2:450302/1‑250 (MQ=255)                                                                                              aGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGt                                                                                                                                         <  1:394296/251‑1 (MQ=255)                                                                                                tGCTGAAAAAAGTAAACCCGGAACCCACGCTGTTCTTGGTAGCATCTAAAAACTTCACCACACAGGAAACTATGACCAACGACCATAGCGCGCGTGACTGGCTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACAATTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTc                                                                                                                                        >  2:166710/1‑250 (MQ=255)                                                                                                  cTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATATAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGc                                                                                                                                     >  1:355234/1‑251 (MQ=255)                                                                                                           gtaaaaccggaacacaagctgataatGGTAGCCTCTAAAACCTTCACCCCTCACGAAACTCTGACAACCGCCCATAGAGCCCGTGACTCGTTCCTGAACGCGGCACGTGCTGAAAAACACGTTGCAAAACAATTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGCGTTCTGGGCCTGGGTTGGCGGCCGTTACTCTTTGTGGTCCGCGGTTGGCCt                                                                                                                            >  2:242224/26‑251 (MQ=255)                                                                                                            aaaaCCCTGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTg                                                                                                                           <  2:161254/250‑1 (MQ=255)                                                                                                             aaaCCCGGAAACCACGCTGTTATTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGt                                                                                                                          >  1:78856/1‑251 (MQ=255)                                                                                                                       aCCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTc                                                                                                                                                                           <  2:2986/192‑1 (MQ=255)                                                                                                                       aCCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTc                                                                                                                                                                           >  1:2986/1‑192 (MQ=255)                                                                                                                              tGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGCCTGGGTTGGCGGCAGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGCTTGTTCTCTCCATc                                                                                                         >  2:22934/1‑251 (MQ=255)                                                                                                                                    tGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCttt                                                                                                   >  1:350725/1‑251 (MQ=255)                                                                                                                                         gCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGCTTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATa                                                                                               >  2:287656/1‑250 (MQ=255)                                                                                                                                          cATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAAc                                                                                             >  2:266669/1‑251 (MQ=255)                                                                                                                                                    cTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCACCATGTTCGGGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAAc                                                                                   >  2:428583/1‑251 (MQ=255)                                                                                                                                                      ttccctcTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCTGCAGGTGATGAAAAACACGTTGCAAATCATTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTTCCAACATGTTCGAGTTCTGGGACTGGGTTTGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAAc                                                                                   <  1:107170/243‑1 (MQ=255)                                                                                                                                                         ccacTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGCCTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCtt                                                                              >  2:68014/1‑251 (MQ=255)                                                                                                                                                           acTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTc                                                                            <  2:372420/251‑1 (MQ=255)                                                                                                                                                              cAGGAAACTATGACCAACTCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCATAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGTGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGGCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTGCCg                                                                          <  1:264960/250‑1 (MQ=255)                                                                                                                                                              cAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTGCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTTGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCg                                                                          <  1:264937/250‑1 (MQ=255)                                                                                                                                                               aGGAAACTATGACCACCGCCCATAGCGCGCGTGACTGGTTCCTGAACGCGGCGGGTGATGCAAAACCCGTTGCAACACACTTTGTGGCGCTTTCCACCAATGCCACAGCCGTCGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGGTCTCTCCTTCGGCTTTGATACCTTCGTTGACCTGCTTTCCGgc                                                                        >  2:58594/1‑251 (MQ=255)                                                                                                                                                                                        cgcgcgTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCGGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTtct                                               >  1:464024/1‑251 (MQ=255)                                                                                                                                                                                          cgtgTGACTGGTGCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACGGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTcc                                             <  2:156960/248‑1 (MQ=255)                                                                                                                                                                                                        cTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTCCACCACGCCTGCCg                                >  1:147270/1‑250 (MQ=255)                                                                                                                                                                                                               cttcAGGTGATGAAAAACACGATGCAAAACACTTTGAGGCGATTTCCACCAATGCCAAAGCCGATGGAGAGCTTGGTATTGATACTGCCAACATCTTCGAGTTAAGGGACTGGGTTGGCGGACGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCGATAGGCTTTGATAACTTCGTTGAACTGCTTTCCGACTCGCACGCGATGGACACGCCTTTCTCCACCACGCCTGCCgaggtcaac                        >  2:115108/4‑245 (MQ=255)                                                                                                                                                                                                                          gaAAAACACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTCCACCACGCCTGCCGAGAAAAACCTGCCTGTAct             >  1:320264/1‑251 (MQ=255)                                                                                                                                                                                                                           aaaaaCACGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTCCACCACGCCTGCCGAGAAAAACCTGCCTGTActg            >  1:270280/1‑251 (MQ=255)                                                                                                                                                                                                                               acacGTTGCAAAACACTTTGTGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTCCACCACGCCTGCCGAGAAAAACCTGCCTGTACTGCTgg        <  1:50999/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggCAAATCACTTTGTGGCGCTTCCCAACAATGCTAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGTCTGGGTTGGCGGCCGTTACTCTTGGGGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTCCACCACGCCTGCCGAGAAAAACCTGCCTGTACTGCTGGCGCTGa  <  1:58594/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             CGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCCATCGGCTTTGATAACTTCGTTGAACTGCTTTCCGGCGCACACGCGATGGACAAGCATTTCTCCACCACGCCTGCCGAGAAAAACCTGCCTGTACTGCTGGCGCTGA  >  NC_000913/4234237‑4234715

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