Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_010473	4,658,941	0	A	T	30.8%	64.6 / 43.3	40	intergenic (+196/+224)	yjjV/yjjW	metal‑dependent hydrolase/YjjW family glycine radical enzyme activase
Reads supporting (aligned to +/- strand):  ref base A (18/9);  new base T (2/10);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 5.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACA  >  NC_010473/4658761‑4659130                                                                                                                                                                                                         |                                                                                                                                                                                                             cACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGT‑C‑GAACCGTAGGCCGGATAA‑GGCg                                                                                                                                                                                          >  1:22316S5/1‑200 (MQ=17)  aCGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTACGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAAAGT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCg                                                                                                                                                                                          >  1:12584S25/1‑200 (MQ=17)      cgcATCCGGCAGTCGTGCGCA‑GATGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGC‑GTGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCAcg                                                                                                                                                                                    <  1:32027S6/200‑1 (MQ=2)           tccggcagtcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTA‑AA‑TT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCa                                                                                                                                                                               <  1:1959S76/182‑1 (MQ=2)                 agtcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCg                                                                                                                                                                           <  1:22890S9/188‑1 (MQ=2)                  ttG‑GCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGAATGCCGACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGCGCTAAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATc                                                                                                                                                                             <  1:32573S10/200‑1 (MQ=17)                   tcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGc                                                                                                                                                                         >  1:280S3/11‑200 (MQ=2)                   tcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGc                                                                                                                                                                         <  1:24231S32/190‑1 (MQ=2)                    g‑gCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAg                                                                                                                                                                       >  1:26551S85/1‑200 (MQ=32)                     gtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAg                                                                                                                                                                       <  1:3180S58/192‑1 (MQ=2)                      gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑t                                                                                                                                                                   >  1:6875S15/1‑200 (MQ=32)                        cgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑g                                                                                                                                                               <  1:30953S3/195‑1 (MQ=2)                        cgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GACGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑g                                                                                                                                                               <  1:28620S14/195‑1 (MQ=2)                          acaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑gc                                                                                                                                                            >  1:25996S7/1‑200 (MQ=32)                             aTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑ATCCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑gtgcg                                                                                                                                                           >  1:285S11/1‑196 (MQ=2)                              tGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAA‑AT‑‑TT‑ACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑gtgctggga                                                                                                                                                       <  1:27323S24/200‑9 (MQ=2)                                  tGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATg                                                                                                                                                   >  1:17914S23/1‑200 (MQ=32)                                    aTGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑GTGCGC‑agatgcc                                                                                                                                                  <  1:496S2/200‑7 (MQ=2)                                     tGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑gtgctgggatgggc                                                                                                                                                  >  1:10160S18/1‑187 (MQ=0)                                                aCGCGTCTTATCAGGCCTACAAATT‑C‑GCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAAAGGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑GTGCGC‑A‑GATG‑CCTGATGCGACGCt                                                                                                                                    <  1:19201S149/200‑1 (MQ=2)                                                 cgcgTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCtt                                                                                                                                   >  1:8183S8/1‑200 (MQ=32)                                                               ccTACAGGTT‑T‑GTGC‑CGGACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCa                                                                                                                    <  1:15240S23/200‑1 (MQ=17)                                                                cTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCCT‑ATCAt                                                                                                                   >  1:11331S4/1‑200 (MQ=17)                                                                  aCAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGc                                                                                                                 <  1:8104S94/200‑1 (MQ=32)                                                                   cAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGcc                                                                                                                <  1:27075S58/200‑1 (MQ=32)                                                                                  aaccgttgggCCGGATAAGGCATTAGCGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTCTTGTGGTTTATTATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TCCGG‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGTCTAc                                                                                                             >  1:4332S65/9‑200 (MQ=14)                                                                                                aTAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGcc                                                                                   <  1:4727S8/200‑1 (MQ=32)                                                                                                     gCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATa                                                                            >  1:11548S43/1‑200 (MQ=32)                                                                                                       gTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑AAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAg                                                                          >  1:21091S46/1‑200 (MQ=17)                                                                                                           aCGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGt                                                                      >  1:9452S39/1‑200 (MQ=32)                                                                                                                               gcgcACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAACG‑CCCGACGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGt                                                   <  1:32720S16/200‑1 (MQ=17)                                                                                                                                     aaTGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑A‑GGGCGTTAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGgcg                                              >  1:26832S39/1‑200 (MQ=17)                                                                                                                                      aTGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGca                                            >  1:7900S44/1‑200 (MQ=32)                                                                                                                                       tGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTTTATCACGCC‑GCAGAGGTTGGGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGca                                            <  1:15666S6/200‑1 (MQ=17)                                                                                                                                        gCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGcaca                                          <  1:8188S6/200‑1 (MQ=32)                                                                                                                                                   gACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCCCGAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGa                                  >  1:323S39/1‑200 (MQ=17)                                                                                                                                                             gcgTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTg                       >  1:8761S60/1‑200 (MQ=32)                                                                                                                                                                                         gT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTC‑‑‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAgcgagtct  <  1:764S89/200‑8 (MQ=17)                                                                                                                                                                                          t‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCGTCCGGGAGGAATC‑‑G‑GCGC‑‑CGCTGTAATAATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTCAGGCCCGG‑ATAAGGCATTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGc       >  1:24229S37/1‑200 (MQ=14)                                                                                                                                                                                          t‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa  >  1:19333S17/1‑200 (MQ=32)                                                                                                                                                                                                         |                                                                                                                                                                                                             CACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAG‑‑‑TC‑‑G‑GCGC‑ATAATG‑CCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGT‑AGGCC‑GG‑ATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACA  >  NC_010473/4658761‑4659130

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 49 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑