Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,670	0	T	A	38.1%	17.9 / 17.3	50	pseudogene (934/1617 nt)	glvC	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (8/18);  new base A (13/3);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 3.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGG‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTG  >  NC_000913/3862456‑3862907                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑GGGT‑TGGGG‑AAACAGCGAG‑CCCGCCCCCTTTA‑ACCg                                                                                                                                                                                                                            >  2:222021/1‑239 (MQ=255)  cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGT‑GCGGG‑GAACAGCAAA‑CCCCCCACTTTTA‑ccc                                                                                                                                                                                                                             >  1:282506/1‑237 (MQ=255)      gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑AGAT‑GCGGG‑AATAAGAAAA‑CCCGCCACTTTTT‑CCCCGtttt                                                                                                                                                                                                                       >  1:1018197/1‑239 (MQ=255)               cgtggggtcccAGTTTGGCGGAAAACTCTAGACAAGCACACCACCCCCGTTTTCGCCCACCACAGAAAAGAAAAACGTAACGGGCGACATGGAGGCCGCCAGCAGGGCGTGCACCAAAAACAGCAGGGGTAAAAGGAAGAGGAAGAAAAATCCCG‑GGGGTTCGGTAATCCCCACCAGGTGGGG‑GTTTAG‑GGGG‑GGGGG‑AACCAGCAAC‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCgg                                                                                                                                                                                                              <  2:282506/230‑1 (MQ=255)                                                          ccGCCAATGTTCCCCACACCACAAAAGAAAAACATTGCGGCCGAAATTGAGGCCGCCGGCACGGCGTGTACCACAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCG‑GCGGTTCGGTAATCCCCCCCAGCA‑GGG‑GGTTAA‑GGGT‑GGGGGAAATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCgg                                                                                                                                                                                                              <  2:236472/196‑1 (MQ=255)                                                                atgtcccccACCACCCAAAAGAAATACATTCCGGCCAACATGGGGCCCCCCAGCCCGGCGTCCCCCGCAAACAGCAACGGGAAAATGAACAGAAAGGTAAATCCCA‑GCGTTCCGTTATTCCCCCCCAGCA‑GGG‑GGTTAAGGGTT‑GGGGA‑AATAAGAAAC‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTa                                                                                                                                                                                                   <  2:214413/194‑1 (MQ=255)                                                                aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGT‑GCGGG‑GATCAGCAAA‑ACCCCCCCCTTTT‑CCCCCTTTTCCGGGGCGGGCG‑GGAAAGAC‑AACG‑CTAAAGAAAAACCCACCGGA‑CC‑AAAAAatttttt                                                                                                                                                        >  2:317404/1‑237 (MQ=255)                                                                    tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑AGGT‑GCGGG‑GAACAGCCAG‑CCCGCCCCCTTTA‑CCCGATTTTTCGGGGCGGCAG‑TGGAgga                                                                                                                                                                                                   >  1:236472/1‑194 (MQ=255)                                                                     gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCCACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGCGGGTTTA‑GGTT‑‑CGGG‑GGACCGCCAA‑CCCCCCACTTTTT‑CCCCCTTTTTCCGGGCGGCCG‑GGGAGGACAACCC‑CTAAAAAAATCCCCCCGGCC‑CC‑AAAAAttttttggta                                                                                                                                                     >  2:172302/1‑233 (MQ=255)                                                                          ccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGT‑GCGGG‑GAACAGCAAA‑CCCCCCACCTTTA‑CCCGATTTTTCCGGGCGGCAG‑GGAAGGAA‑AACG‑CCaaa                                                                                                                                                                                       >  1:214413/1‑200 (MQ=255)                                                                             caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑AGTT‑GTGGG‑AAGAAGCGAG‑CCCCCCCCCTTTT‑CCCGGGTTTTCGGTGCGGCAG‑GGAAGTAC‑AACG‑CCAAACAAATAACCACGGGC‑CC‑CAAGATTTTTTAGTTACCGTGccac                                                                                                                                           >  2:783771/1‑236 (MQ=255)                                                                                        tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGAGAACCGCAAACAGCAACGGAGAAATGAACAGGAAGGCAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTGTA‑GGGA‑TCGGG‑AAACACCAAG‑GCCGCCACCTTTA‑ACCGCTTGTCTCGGGCGGGGG‑AGAAATACAATGA‑CAAAAAAAATGCCTAGGTCA‑CA‑ACAGCTTTTT‑AGGTACCGTGTCAGGCaaaaaat                                                                                                                                  >  2:536223/1‑234 (MQ=255)                                                                                          cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑AGGT‑GCGGG‑GAACAGCAAG‑CCCGCCCCCTTTA‑CCCGCATTTCCCGGGCGGCAG‑GGAAGTAC‑ATCG‑CTAAAAAAAAGCCCACCGCA‑CC‑AAAAAATTTTTAGGTAACGCGCAGGGGAAAAACGCCt                                                                                                                               >  1:678/1‑238 (MQ=255)                                                                                               cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGG‑GATCAGCAAG‑CCCGCCCCTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAAAGCCCACCGCA‑CC‑AAAGAATTTTGAGTTACCGTGCAGGGCAAAACCGCCCTCCcgg                                                                                                                         >  1:111066/1‑239 (MQ=255)                                                                                               cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑GGGT‑TCGGG‑GAACAGCAAGCCCCCCCCCTTTTA‑CCC‑ATTTTTCGGTGCGGGAGCGGAAGACA‑ACCG‑CTAAA‑AAAAGCCCCCCGCA‑CCAAAAGATTTTTTATTTACGGTCGGGGGACAAACGCCCTTCCgg                                                                                                                          >  2:349646/1‑239 (MQ=255)                                                                                               cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCAAGGC‑GGTTAA‑AGTT‑GCGGG‑AATCAACAAG‑CCCCCCCCCTTTT‑CCCCGTTTTTCGGTGGGGGAG‑AGGAGGAA‑ATCG‑CTCAAGAAAATGCCACGGCG‑AC‑AAAGATATTTTAGGTAACGTGCAGGGCCAAAAAGCCCTccccg                                                                                                                         >  2:465469/1‑236 (MQ=255)                                                                                                  tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACCGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTT‑GCGGG‑AAACGCAAAG‑CCCCCCCCTTTT‑‑CCCGGTTTTTCCGTGCGGCCG‑GGGAGTAC‑AACC‑CTAAAAAAATACCCCCCGCC‑CC‑CAAGAATTTTTAGTTACCCTGGAGGGCAAAAACGCCTTCCCCGaac                                                                                                                      >  2:921200/1‑239 (MQ=255)                                                                                                  tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAA‑GGTT‑GGGGG‑AA‑AAACAAA‑CCCCCCCCCTTTA‑CCCCGTTTTTCGGTGGGGGCG‑GGGAGGAA‑AACGCCTAAAAAAAACCCACCGCCC‑CC‑AAAGATTTTTAGGTTACGCGCCGGGGCAAACCCCCCTCCCGGGac                                                                                                                       >  2:583944/1‑238 (MQ=255)                                                                                                     aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCCGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑AGGT‑GCGGG‑GAACAGCAAG‑CCCCCCACTTTTA‑CCCGATTTTTCGGTGCGGCAG‑TGAAGGAC‑AACG‑CCAAAAAAAAGCCCACGGCA‑CC‑CAAGATTTTTTAGTTACCGTGCAGGGCaaaaac                                                                                                                                   >  1:937956/1‑221 (MQ=255)                                                                                                                gccAGCACGGCGTGTCCCGCAAACAGCAAGGGTAAAATGAACAGGAAGGTAAATCCCA‑GCGGTTCGGTAATCCCCCCCAGCAGGGC‑GGTTAA‑GGTT‑GGGGG‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTa                                                                                                                                                    <  1:481297/195‑1 (MQ=255)                                                                                                                     cACGGCGTGTACCGCAAACAGAAACGGTAAAATGAACAGGAAGGTAAATCCCA‑GCGTTTCGGTAATCCCCCCCAGCGTGGC‑GGTTAA‑GGTT‑GGGGA‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGaa                                                                                                                                                                                    <  1:913278/160‑1 (MQ=255)                                                                                                                         gCGTGTCCCAAAAACAGCAAGGGAGAAAAGAACAGGAAGGTAAATCCCG‑GCGGTGCGGTAATCCCCACCAGCAGGGC‑TGTTAA‑GGTG‑GCGGG‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCgg                                                                                                   <  1:1011111/235‑1 (MQ=255)                                                                                                                          ggtgacccaaaaaaAGAAAGGGTAAAAGGAAAA‑GAAGGGAAATTCCAGGCGTTCCGGAATTTCCCCCCAGCAGGGC‑GGTTAA‑GGTG‑GGGGG‑AATCAGAAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGACAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaaca                                                                                                                     <  1:563612/206‑1 (MQ=255)                                                                                                                          cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑AGGT‑TCGGG‑GAAAAACAAA‑CCCGCCCCCTTTA‑ACCGGTTTTTCCGTGGGGCAG‑TGAAGGAA‑ATCG‑CCAAAGAAATACCCCCCGCACCC‑AAAAATTTTTGATTAACCggtcggg                                                                                                                                           >  2:481297/1‑189 (MQ=255)                                                                                                                              tACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑GGGT‑GCGGG‑GAACAGGAAA‑ACCCCCCCTTTTACCCCGATTTT‑CGGGGCGGGAG‑GGAAGGAC‑AACC‑CCAAAAAAAAGCCCCCCCCC‑CC‑CACGAATTTTTGGGTAGCCTGCAGGGGAAAAAAGCCCTCCCGGGACAACCCTTTCTAGGGCTCCGCGCGCaaaa                                                                                          >  2:417287/1‑238 (MQ=255)                                                                                                                              aaCCGCAAACAGGAGCGGAAAAAAGAACAGGAAGGTAACTTCCC‑GATGTACGGTATTCCCCACCGTCTCGGC‑GTTGCG‑GGGG‑GGGAG‑AAGCGTAACG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGc                                                                                                                                                                                                                <  2:633256/125‑1 (MQ=255)                                                                                                                                    aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑AGGT‑GTGGG‑AATCAGCAAG‑CCCCCCACTTTTA‑CCCGATTTTTCGGGGGGGCAG‑GGAAGTAA‑AACG‑CCAAAAAAAATCCCACGGCC‑CC‑AAAGATTTTTTAGTTAACGGGCAGGGGAAAACCCCCCTCCGGGAAAAAAGATTTTAg                                                                                                           >  2:563612/1‑216 (MQ=255)                                                                                                                                        aGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTA‑GGTT‑GCGGG‑GATCCGCAAT‑ACCGCCACTTTTA‑CCCGCTTTTCCCGTGCGGCAG‑TGAAg                                                                                                                                                                                                     >  1:633256/1‑126 (MQ=255)                                                                                                                                           aaggggaaaaGGAACAGGAAGGTAAATCCCA‑GCGTTTCGGTAATCCCCCCCAGCATGGG‑GGTTAA‑GGTG‑GGGGA‑AATCAGCACG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTg                                                                                <  1:717555/229‑1 (MQ=255)                                                                                                                                                     tGAAGAGGAAGGAAATTCCCG‑GCGTTTCGGTAATCCCCCCCAGCATGGG‑GGTTAA‑GTTG‑GCGGG‑AATCACCACG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:665828/239‑1 (MQ=255)                                                                                                                                                     tGAACGGGAAAGAAATTTCCG‑GGTGTTTGGTTATCCCCCCCAGCATGGC‑TGTTAG‑GGTT‑GCGGA‑AACCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:536223/239‑1 (MQ=255)                                                                                                                                                     tGAACGAGAAGGAAAATTCCA‑GGGTTTCGGTAATCCCCACCAGCGTGGG‑GGTAAA‑GGTG‑GAGGA‑AATCACCAAA‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:692211/239‑1 (MQ=255)                                                                                                                                                            gAAGGTAAATCCCA‑GCGGTTCGTTAATCCCCCCCGGCTTGGG‑GGTT‑A‑GGTT‑GGGGGAAATCAGCACC‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                              <  1:300921/237‑1 (MQ=255)                                                                                                                                                             aaGGTAAATTCCA‑GCGTTTCGTTAATTCCCACCATCTTGGC‑GTTTAG‑GGTT‑GCGGG‑AATCAGAAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:125169/231‑1 (MQ=255)                                                                                                                                                                 tAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGG‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑AGAAGAAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCAGCTCCGCGCTCAAACTGAACTCCTGCAGAAGCTGCGCCCAGTACATCTGa                                                       >  1:723256/1‑239 (MQ=255)                                                                                                                                                                    tgtCCG‑GCGGTGCGTTAATCCCCACCAGCTGGGC‑GTTAAG‑GGGT‑GGGGG‑ACTAAAAAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                    <  1:694836/237‑1 (MQ=255)                                                                                                                                                                     tcccA‑GCGTTTCGATACTCCCCCCCAGCAGGGC‑GGTTAG‑GGGG‑GCGGG‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgc                                                   <  2:576201/237‑1 (MQ=255)                                                                                                                                                                           gCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGT‑GCGGG‑AAACAGCCAA‑CCCCCCCCCTTTT‑CCCCA‑TTTCCCGTGCGGGAG‑GGGAGGAC‑ATCC‑CTCAAAAAAAGCCCACCGCC‑CC‑AAAAATTTTTTAGTTACCCTGCAGGGCAAAACCGCCCTCCCGGAACAACGATTTTAGCGGCCCCGCGCCCCAACAGAACACCTCTCGAAGCTGCGCCCCGCAAATCTCAATGCCCCCt                                              >  2:1024254/1‑238 (MQ=255)                                                                                                                                                                                           ccccaGCATGGC‑GGTTAA‑GGTT‑GCGGG‑AATAAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                              <  1:926383/237‑1 (MQ=255)                                                                                                                                                                                           ccccaGCAGGGC‑GGTTAA‑GGTG‑GGGGG‑AACAAGAAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                              <  1:172302/237‑1 (MQ=255)                                                                                                                                                                                           ccccaGAATGGG‑GGTTGA‑TGTC‑GAGGG‑ATACACCAGC‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                              <  1:727156/237‑1 (MQ=255)                                                                                                                                                                                           caccggcggggC‑GGTTAG‑GGTG‑GCGGA‑AATCAGAAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                              <  1:80671/230‑1 (MQ=255)                                                                                                                                                                                                      ggggTTAG‑GGTG‑GGGGA‑AACCACCAAC‑CCCGCCACTTTTA‑CCCGA‑TTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAGATGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                    <  1:1024254/236‑1 (MQ=255)                                                                                                                                                                                                       gggTTAA‑GGTG‑GCGGG‑AACCAGCAAC‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                    <  1:921200/237‑1 (MQ=255)                                                                                                                                                                                                        ggTTAA‑GGTTGGGGGG‑AAACAAAAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                    <  1:349646/238‑1 (MQ=255)                                                                                                                                                                                                         gtgga‑GGTT‑GGGGA‑AAACGACAAC‑CCCGACACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGGCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGGAAAACCGCCTGCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACt                  <  1:77092/234‑1 (MQ=255)                                                                                                                                                                                                                 gtcccGGA‑AATCAACAAC‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCCCGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt            <  2:208560/234‑1 (MQ=255)                                                                                                                                                                                                                        g‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa     <  1:783771/238‑1 (MQ=255)                                                                                                                                                                                                                         aaaTCAGCCAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa     <  1:465469/237‑1 (MQ=255)                                                                                                                                                                                                                           aTCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtg  >  1:231369/1‑239 (MQ=255)                                                                                                                                                                                                                            tCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCa                                                                                                            >  1:578489/1‑132 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCA‑GCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGG‑AATCAGCAAG‑CCCGCCACTTTTA‑CCCGATTTTCCGGTGCGGCAG‑TGAAGTAC‑ATCG‑CTAAAGAAATGCCCACGGCA‑CC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTG  >  NC_000913/3862456‑3862907

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