Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	3,382,603	0	A	T	27.8%	27.3 / 11.1	18	L193Q (CTG→CAG)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base A (7/6);  new base T (3/2);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGATCACCTGCTTCGAGAAGGCTTCCCAGTCGAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGACCAGGCGGGCGTGCACGCGCTT  >  AM260479/3382472‑3382733                                                                                                                                    |                                                                                                                                   cGATCACCTGCTTCGAGAAGGCTTCCCAGTCGAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCTgg                                                                                                                                  >  1:832311/1‑134 (MQ=255)           cTTCGAGAAGGCTTCCCAGTCGAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCTGGcgcgcgcga                                                                                                                         <  2:832311/133‑4 (MQ=255)                agaAGGCTTCCCAGTCGAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTg                                                                                                                   <  1:222178/134‑1 (MQ=255)                agaAGGCTTCCCAGTCGAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTg                                                                                                                   <  2:494192/134‑1 (MQ=255)                    ggCTTCCCAGTCGAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCa                                                                                                               >  2:322011/1‑134 (MQ=255)                                gAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGccc                                                                                                    <  1:882105/133‑1 (MQ=255)                                   gCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCg                                                                                                >  1:837777/1‑134 (MQ=255)                                                        cgcgTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGggctggccggctgg                                                                            <  2:837777/133‑1 (MQ=255)                                                          cgTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCTGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGgctggccggctggccg                                                                         >  1:582917/1‑134 (MQ=255)                                                            ttGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTg                                                                       <  1:234058/134‑1 (MQ=255)                                                                             gcCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCg                                                      >  1:1001244/1‑134 (MQ=255)                                                                                                 gcTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCTGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTg                                  >  2:792032/1‑134 (MQ=255)                                                                                                           cTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGa                        >  1:767453/1‑134 (MQ=255)                                                                                                               aTGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGACCAg                    >  2:477023/1‑134 (MQ=255)                                                                                                                         tCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGACCAGGCGGGCGTg           >  2:462328/1‑133 (MQ=255)                                                                                                                                cgcCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGAc                       <  1:680261/114‑1 (MQ=255)                                                                                                                                cgcCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGAc                       >  2:680261/1‑114 (MQ=255)                                                                                                                                 gcCTGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGACCAGGCGGGCGTGCACGCGCtt  <  1:792032/134‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   CGATCACCTGCTTCGAGAAGGCTTCCCAGTCGAAGCTCGGGTAGGCCGACAGGTGCGCGTTGTAGTTGCCCACCGCGCCGTTCATCTTGCCCAGCAGCTCGACCTGCTCGATGCGCTGGATCGCGCGCGCCAGGCGCGCGGCCACGTTGGCCATTTCCTTGCCCAGCGTGGTCGGGCTGGCCGGCTGGCCGTGGGTGCGCGACAGCATCGGCTGCGCGGCATTGAGCTTGGCCAGCTCGACCAGGCGGGCGTGCACGCGCTT  >  AM260479/3382472‑3382733

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