Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_010473	4,658,921	0	T	G	34.2%	68.9 / 48.7	41	intergenic (+176/+244)	yjjV/yjjW	metal‑dependent hydrolase/YjjW family glycine radical enzyme activase
Reads supporting (aligned to +/- strand):  ref base T (18/9);  new base G (4/10);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTG  >  NC_010473/4658730‑4659109                                                                                                                                                                                                                  |                                                                                                                                                                                                           gttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CAGACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACaa                                                                                                                                                                                                    <  1:764S90/200‑1 (MQ=17)  tttAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑g                                                                                                                                                                                                  >  1:24229S36/1‑200 (MQ=32)   ttAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGGCATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑g                                                                                                                                                                                                  >  1:19333S16/1‑200 (MQ=17)         acccGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTGGCATGCTTCGCAGTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACaaattct                                                                                                                                                                                               >  1:31922S16/3‑195 (MQ=2)                                cACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGT‑C‑GAACCGTAGGCCGGATAA‑GGCg                                                                                                                                                              >  1:22316S5/1‑200 (MQ=17)                                 aCGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTACGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAAAGT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCg                                                                                                                                                              >  1:12584S25/1‑200 (MQ=17)                                     cg‑cATCCGGCAGTCGT‑‑GCG‑CA‑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)                                                  ttG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGAATGCCGACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGCGCTAAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATc                                                                                                                                                 <  1:32573S10/200‑1 (MQ=17)                                                    g‑‑‑gCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAg                                                                                                                                           >  1:26551S85/1‑200 (MQ=32)                                                    agtcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCg                                                                                                                                               <  1:22890S9/188‑1 (MQ=2)                                                      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)                                                        gtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAg                                                                                                                                           <  1:3180S58/192‑1 (MQ=2)                                                        gcg‑cACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGt                                                                                                                                          >  1:6875S15/1‑200 (MQ=32)                                                           cgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑g                                                                                                                                      <  1:30953S3/195‑1 (MQ=2)                                                           cgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GACGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑g                                                                                                                                      <  1:28620S14/195‑1 (MQ=2)                                                             acaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑gc                                                                                                                                   >  1:25996S7/1‑200 (MQ=32)                                                                aTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑ATCCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑gtgcg                                                                                                                                  >  1:285S11/1‑196 (MQ=2)                                                                 tGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAA‑AT‑‑TT‑ACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑gtgctggga                                                                                                                              <  1:27323S24/200‑9 (MQ=2)                                                                     tGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATg                                                                                                                          >  1:17914S23/1‑200 (MQ=32)                                                                       aTGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑GTGCGC‑agatgcc                                                                                                                         <  1:496S2/200‑7 (MQ=2)                                                                        tGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑gtgctgggatgggc                                                                                                                         >  1:10160S18/1‑187 (MQ=0)                                                                                   aCGCGTCTTATCAGGCCTACAAATT‑C‑GCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAAAGGCGTTCACGCCGCATCCGGCAGTC‑‑GTGCGC‑A‑GATGCCTGATGCGACGCt                                                                                                            <  1:19201S149/200‑1 (MQ=2)                                                                                    cgcgTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCtt                                                                                                           >  1:8183S8/1‑200 (MQ=32)                                                                                                  ccTACAGGTT‑T‑GTGC‑CGGACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCa                                                                                            <  1:15240S23/200‑1 (MQ=17)                                                                                                   cTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCCT‑ATCAt                                                                                           >  1:11331S4/1‑200 (MQ=17)                                                                                                     aCAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGc                                                                                         <  1:8104S94/200‑1 (MQ=32)                                                                                                      cAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGcc                                                                                        <  1:27075S58/200‑1 (MQ=32)                                                                                                                     aaccgttgggCCGGATAAGGCATTAGCGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTCTTGTGGTTTATTATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTCCGG‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGTCTAc                                                                                     >  1:4332S65/9‑200 (MQ=14)                                                                                                                                   aTAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGcc                                                            <  1:4727S8/200‑1 (MQ=32)                                                                                                                                        gCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATa                                                       >  1:11548S43/1‑200 (MQ=32)                                                                                                                                          gTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑AAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAg                                                     >  1:21091S46/1‑200 (MQ=17)                                                                                                                                              aCGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGt                                                 >  1:9452S39/1‑200 (MQ=32)                                                                                                                                                                  gcgcACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAACGCCCGACGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGt                              <  1:32720S16/200‑1 (MQ=17)                                                                                                                                                                        aaTGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑A‑GGGCGTTAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGgcg                         >  1:26832S39/1‑200 (MQ=17)                                                                                                                                                                         aTGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGca                       >  1:7900S44/1‑200 (MQ=32)                                                                                                                                                                          tGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTTATCACGCC‑GCAGAGGTTGGGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGca                       <  1:15666S6/200‑1 (MQ=17)                                                                                                                                                                           gCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGcaca                     <  1:8188S6/200‑1 (MQ=32)                                                                                                                                                                                      gACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGa             >  1:323S39/1‑200 (MQ=17)                                                                                                                                                                                                gcgTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTg  >  1:8761S60/1‑200 (MQ=32)                                                                                                                                                                                                                  |                                                                                                                                                                                                           GTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTT‑ATCATGCCTACAG‑GTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTG  >  NC_010473/4658730‑4659109

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