Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,502,761	0	G	C	25.9%	124.9 / 14.0	60	A224G (GCG→GGG)	insG	IS4 family putative transposase
Reads supporting (aligned to +/- strand):  ref base G (14/29);  new base C (14/1);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 4.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTAT  >  NC_000913/4502513‑4503010                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                          aCAGATCCCCCATTTCCCCCCCGGGGAAGCGCATGGCGTCCGTCATCGACGGCAGCAGATGGCATACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGGCCCGGCCCCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTT‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGcc                                                                                                                                                                                                                                                                         <  1:71854/251‑1 (MQ=255)               tctcCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑GG‑TTTTCCCCCCAAGCTCCAgcgctt                                                                                                                                                                                                                                                           >  1:370114/1‑246 (MQ=255)                      cGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGt                                                                                                                                                                                                                                                 <  1:295547/251‑1 (MQ=255)                       ggggAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTAATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑AG‑TTTTTCCGCCCCGCTTCAGGCCA‑TTTaaaagc                                                                                                                                                                                                                                                  >  1:420005/1‑247 (MQ=255)                         ggAAGCGCATGGCGTCCGTCATAGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTAACGGTCAGCAGGCGGGCCGTAACTTCCTTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCCGCTTAACCAGATGATCGACCTTACCCAATTTTATG‑ATCTATTCATACTGCGCTCCCTTACTGAGAGGGATATTCCAGTGGCGG‑TG‑GTCTCTCGCCAgacccccagcctttatcaatacc                                                                                                                                                                                                                                                 >  2:105802/1‑229 (MQ=255)                             gctcATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTg                                                                                                                                                                                                                                          <  1:149749/248‑1 (MQ=255)                                             aTCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTTCCCCCCCGCCCCCGGC‑A‑TTTAA‑CAGTCCCCGTGGGGAAA‑AA‑CCTTTAtt                                                                                                                                                                                                                        >  1:360964/1‑250 (MQ=255)                                               cGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCa                                                                                                                                                                                                                      <  1:809845/251‑1 (MQ=255)                                                gACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCCTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑t                                                                                                                                                                                                                    <  1:317461/251‑1 (MQ=255)                                                  cGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑GTCTTCCCCCCCGCTCCCGGC‑A‑TTTTA‑CAATCTCCGTGGGGAAT‑AA‑CCTTTATTCA‑ttt                                                                                                                                                                                                                  >  1:407357/1‑250 (MQ=255)                                                   ggcagcagATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑aa                                                                                                                                                                                                                <  1:662850/249‑1 (MQ=255)                                                     cagcagATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AAc                                                                                                                                                                                                               <  1:693329/250‑1 (MQ=255)                                                        cagATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTTTTCCCCCAGGCCCCCGGC‑C‑TTTTT‑AAGTCCCCGGGAGTAAT‑AA‑CCTTTTTCCC‑TT‑TACGTcc                                                                                                                                                                                                           >  1:137307/1‑250 (MQ=255)                                                               cAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTTTCCCCCCCGGCTCCAGGC‑A‑TTTTA‑CACACCCCGTGAGTAATTAA‑CTTTTTTCCT‑TT‑AActccggggtgt                                                                                                                                                                                                     >  1:115561/1‑241 (MQ=255)                                                                  aCTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATc                                                                                                                                                                                               <  1:584764/251‑1 (MQ=255)                                                                    ttttCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCGGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCTCCAGGCTCCAGGC‑A‑ttt                                                                                                                                                                                                                                                        <  1:701918/199‑1 (MQ=255)                                                                                    aCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGGTG‑TT‑TCCCGCCAGGCTACAGGC‑A‑TTTAA‑CAGTCCCAGAGAGTAAT‑AA‑ACTTTTTTCA‑TT‑AACGTTC‑G‑AGTGGTATCGGCGGGTTGGTCTTT‑aa                                                                                                                                                                             >  1:606238/1‑250 (MQ=255)                                                                                                 gggCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGc                                                                                                                                                              <  1:50694/251‑1 (MQ=255)                                                                                                        cACTTCATTTCCCCGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCATGGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑ttt                                                                                                                                                       <  2:845555/250‑1 (MQ=255)                                                                                                        cACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑ttt                                                                                                                                                       <  1:626640/250‑1 (MQ=255)                                                                                                                        cccGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTTCCCCCCGGCTCCAGGC‑A‑TTTAAACA‑TCCCAGAGGGAAAT‑AA‑CCCTTATCCA‑TT‑AAAGTCA‑G‑AGAGTTTTAGCCGGGTTGTTCTTT‑CA‑GTTGCTCAGCAAGCTCA‑TTTTTGCTGTTTTT‑CATcct                                                                                                                                     >  1:350139/1‑249 (MQ=255)                                                                                                                        cccGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TA‑TTCTCCCCCCCGGCTCCCGGC‑A‑TTTTA‑CAGACCCAGCGAGTAAA‑AAACCTTT‑TCCA‑TTTAAC‑TCA‑G‑AGTGTTTTCTCCCGTTTGTTTTTTTAA‑ATTG‑TCAACAAGCTCA‑TTTTCTCTGTTCTT‑CATCGt                                                                                                                                     >  1:370340/1‑251 (MQ=255)                                                                                                                         ccGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCCCCCCGCTCCCGGCCA‑TTAAA‑CA‑TCCCCGTGAGGAAT‑AA‑CCCTTATCCA‑TT‑AACGTCA‑G‑AGGGTTATCGCCGGGTTGTTCTAT‑AA‑GTTGGTCAGCAAACTCA‑TTTTCGCTGTTCTT‑CCTCGTg                                                                                                                                    >  1:598329/1‑251 (MQ=255)                                                                                                                           gtacccTTTTTTCTGGCCTTCGGGCTGTTTTTCAGCTTCTCCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTGTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCCCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGTCGGTTTGGTCTAT‑AA‑GTGGCGCAGCAAGTTCA‑TTTTCGCTGTTTTT‑CATCGTGcc                                                                                                                                  <  2:786542/246‑1 (MQ=255)                                                                                                                             ccACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCCGGCTCCCGGC‑A‑TTTTA‑CAGTCCCCCTGAGTAAT‑AA‑CCTTTTTCCA‑TT‑TACGTCA‑G‑CGTGTTATCGCCGGGTTGTTTTTT‑TA‑GTTGTTCAGCAAACTCA‑TTTTTGCtttttttcttccccccgc                                                                                                                                 >  1:172078/1‑234 (MQ=255)                                                                                                                              cACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCCCCCTGCCCCCGGC‑A‑TTTCA‑CAGTCCCCGGGAGTAAT‑AA‑CATTTTTCCCCTC‑AAC‑TCA‑T‑AGTGGTTTCGCCGGTGTGTTCTTT‑AT‑GTTTTTCCGCAAGCTCC‑TTTTACCTGttctt‑cctctttccca                                                                                                                                >  2:272615/1‑240 (MQ=255)                                                                                                                                           ccttcGGGCTGGTTTTTATCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AT‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑c                                                                                                                 <  1:741193/247‑1 (MQ=255)                                                                                                                                              gcgtgCTGTTTTTCAGCTTCTCCAGATGATCTCCTTTCCCCAGTTTTCTG‑ATCTCTTCATATTTCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGGCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CAGCGTGCCGAAGGCTGCAGCCGT‑cagc                                                                                                              <  2:870489/247‑1 (MQ=255)                                                                                                                                                gggCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTTTCCCGCCAGGCTCCCGGC‑A‑ATTAA‑AAATCCCCGTGGGTAAT‑AA‑CCTTTTTCCA‑TT‑AACGTCA‑G‑CGGGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAACTCA‑TTTTCGCTGTTCTT‑CCTCGTGCCGAAGGCTGCAGCCGT‑cagc                                                                                                              >  1:26455/1‑249 (MQ=255)                                                                                                                                                    tGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCCCCAAGACCCAGGC‑A‑TTTTA‑CAGTCCCAAGGAGTAAT‑AA‑CCTTTATCCC‑TT‑AACGTCA‑G‑AGAGTTATTCCCGGTTTTTTTTAT‑AC‑ATTGTTCCGCAAGCCCA‑TTTTCCCTGTTCTT‑CTTCGTGcccgatggtgccccccccaaaagatgg                                                                                                         >  2:586888/1‑226 (MQ=255)                                                                                                                                                           cAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCCCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGCGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTc                                                                                                  <  2:172078/250‑1 (MQ=255)                                                                                                                                                                 cACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTcc                                                                                           <  1:865769/250‑1 (MQ=255)                                                                                                                                                                     aGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTTTTCCCCCAGGCTCCCGGC‑A‑TTTTA‑TAGTCCCAGGGAGTAAT‑AA‑ACTTTAACCA‑TTTAAC‑ACT‑GCAGGGTT‑TCGCCGGGTTGTTCTTT‑TA‑GTTGCTCAGCAAGCTCA‑TTTTTTGTGTTTTT‑TTTCGTGCCCAAGGGTGCAGCCGA‑CACCAGATgaggagggaggaggatcgg                                                                                       >  2:795161/1‑233 (MQ=255)                                                                                                                                                                     aGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTTA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATg                                                                                                         <  1:795161/233‑1 (MQ=255)                                                                                                                                                                     aGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTg                                                                                      <  1:163936/251‑1 (MQ=255)                                                                                                                                                                       gtgatCCCCTTT‑CCCAGTTTTCTGGATCTTTTAATTTTGC‑CCCCCCTTCCGTGAGGTTTCATCCAATGGCGG‑TGGTTCCCCCCCAAGGC‑CCAAGC‑A‑TTTAA‑CCGTCCCAAGGAGTTAA‑AA‑CCTTTATCCA‑TT‑AAACTCAAG‑GTTGTTACCGCGGTTTTGTCTTAT‑AA‑‑TTGGTCCGCAAACCCAATTTTTG‑TTTTCCT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCa                                                                                   <  2:566274/250‑1 (MQ=255)                                                                                                                                                                       atgatCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGttctt‑c                                                                                                                                          <  1:281799/199‑1 (MQ=255)                                                                                                                                                                       atgatCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCCCCCGGC‑A‑TTTTA‑AAGTCCCCGAGAGTAAA‑AA‑ACTTTTTCCC‑TT‑AACATCA‑G‑AGTGTTTTCTCCGGGTTTTTTTAT‑AA‑GTTTCTCAGCCAGCTCA‑TTTTTGCTGTTTTT‑TATCGTGCCGGAGGCTGCAGCCGT‑CCGCACAAGGCTGGTccagtccatcttgc                                                                                     >  1:523757/1‑238 (MQ=255)                                                                                                                                                                         gatCGCCTTTACCCAGTTTTCTG‑TTCTCTGCATATGTCGCTCCCTTTCTGAGAGGTATCATCCAGTGGGGT‑TG‑TCCTCCCGCCAGTCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCGGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAg                                                                                  <  2:520445/251‑1 (MQ=255)                                                                                                                                                                                     ccAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTTTCCCGCCAGGCTCCCGGC‑A‑TTTTA‑CAGTCCCCGTGAGGAAT‑AA‑CCCTTTTCCA‑TT‑AACGTTC‑G‑AGTGTTATCGCCGGTTTGTTCTTT‑TA‑ATTTCGCAGCAAGCTCC‑TTTTTGCTGTTCTT‑TCTCGTGCCCGAGGCTGCAGCCGT‑CAACAGAAGGGTGTTCAGGTTCC‑TCTTGCAGACCCATTTTgcc                                                                      >  1:755757/1‑249 (MQ=255)                                                                                                                                                                                      cAGTTTTCTG‑ACCTCTTCATATTGCGCCCCCTTTCTGATAGGTATCACCCAGTGGCGG‑TG‑TTCTCCCTCCAGGCTCCAGGC‑A‑TTTAA‑CAGGCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTG‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCt                                                                    <  2:606238/251‑1 (MQ=255)                                                                                                                                                                                       aGTTTTCTG‑ATCTCTTCATATTTCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGGCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTg                                                                   <  2:242923/251‑1 (MQ=255)                                                                                                                                                                                            ttcttgttgttgcATATTGACCTCCATTTCGGAGAGGTTAAATCCAGTGGTGG‑TG‑TGCTTCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTCTCCA‑TT‑AACGGCA‑G‑AGTGTTATCGTCGGTTGGGTCTAT‑AA‑GGTGCTCAGCAAGCTCA‑GTTTCGCTGTTCTT‑TATCGTGCCTAAGGCTGCAGCCGT‑CAGCATATGGCTGGTCAG‑ATCCACCTGTCAGACC‑ATTTTGTCCTGCgg                                                                <  2:462198/238‑1 (MQ=255)                                                                                                                                                                                             cTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑GTCTCCCCCCAGGCTCCCGGC‑A‑TTTAA‑CAGTCCCAGTGAGGAAT‑AA‑CCTTTATTCC‑TT‑AACGTCA‑G‑AGTGTTATTGCCGGGTTGTTTTAT‑AA‑GGTGCTCCGCACACTCA‑TTTTCGCCGTTCTT‑CATCCTGCCGAAGGCTGCAGCCGT‑CCGCAGAAGGCTGGTCTG‑TTCCCTCTCGCGGACC‑ATTTTTCCCTGCgggg                                                              >  1:494729/1‑249 (MQ=255)                                                                                                                                                                                               g‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGTTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTag                                                            <  1:509993/250‑1 (MQ=255)                                                                                                                                                                                                      ttCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAAGCCCCAGGC‑AATTTAA‑CA‑TCCCAGTGAGGAAT‑AA‑ACTTTATCCC‑TT‑AACGTCC‑G‑AGAGTTATCACCGGTTTGTTCTTT‑AA‑GTTTGTCCGCACGCTCC‑TTTTCGCCGTTCTT‑ACTCGTTCCCAAAGGGGCAGCCGT‑CAGCAAGTGGCTGGTCCG‑GTTCATCTGGCAGGCC‑Ctttttaacttcggggcgagcaccg                                                     >  1:872615/1‑231 (MQ=255)                                                                                                                                                                                                        cATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCgg                                                    <  1:305513/250‑1 (MQ=255)                                                                                                                                                                                                         atatTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCCCCAGGCTCCAGGC‑C‑TTTAA‑AACTCCCCGGGAGTAAT‑AA‑CCCTTTTCCA‑TT‑AACGTCA‑T‑AGGGTTATCACCGGGTTTTTTTAT‑AAAGTT‑CTCAGCAAAGTCA‑ATTTTGCGGTTTTTTCA‑CCTGCCGAAAGATGCAGCCGTACA‑CCGATGGCTGGTCAG‑TTTCATCTGGGAGGCC‑CTTTTGACCTGCgggggaagggcccggg                                                  >  1:43463/1‑238 (MQ=255)                                                                                                                                                                                                          tatTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTTTCCCCCCCGGCTCCAGGC‑A‑TTTTA‑AAGTCCCAGTGAGTAAT‑AA‑CCTTTTTCCA‑TT‑AACGTc                                                                                                                                                                                                            >  1:175059/1‑105 (MQ=255)                                                                                                                                                                                                          tatTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTc                                                                                                                                                                                                            <  2:175059/105‑1 (MQ=255)                                                                                                                                                                                                           atTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTc                                                <  1:322691/251‑1 (MQ=255)                                                                                                                                                                                                                 cTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGGAAT‑AA‑ACTTTTTCCC‑TT‑AACGTCC‑G‑AGGGTTATCGCCGGGTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑GTCCATCTGGCAGACC‑CTTTTGACCTGCGGGGAGAGCGCCGGGTTCCCGGcc                                          >  1:599098/1‑250 (MQ=255)                                                                                                                                                                                                                           gagaGGGATCATCCAGTGGCGG‑TG‑TGCGCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTGCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGc                                <  2:268047/251‑1 (MQ=255)                                                                                                                                                                                                                                    catcCAGTGGCGG‑GG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGGGCCGAAGGCTGCAGCCGT‑CAGCAGATGGTTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAgg                       <  1:505436/251‑1 (MQ=255)                                                                                                                                                                                                                                           tGGCGG‑TG‑TTCTGCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGAAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAAGCTGCACCCGT‑CAGCAGCTGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGTGAGAGCGCCGGGTTCCCGGCAGGTGTCTGGCGGGGGaaaggtggctc                >  2:703373/1‑242 (MQ=255)                                                                                                                                                                                                                                                     ttCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCa                                                                                                                        <  2:830271/142‑1 (MQ=255)                                                                                                                                                                                                                                                     ttCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCa                                                                                                                        >  1:830271/1‑142 (MQ=255)                                                                                                                                                                                                                                                        tcCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGgtg     >  1:204352/1‑251 (MQ=255)                                                                                                                                                                                                                                                        tcCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGgtg     >  1:695554/1‑251 (MQ=255)                                                                                                                                                                                                                                                           cGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCCTGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTAt  >  2:273315/1‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                          ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTG‑ATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGG‑TG‑TTCTCCCGCCAGGCTCCAGGC‑A‑TTTAA‑CAGTCCCAGTGAGTAAT‑AA‑CCTTTATCCA‑TT‑AACGTCA‑G‑AGTGTTATCGCCGGTTTGTTCTAT‑AA‑GTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTT‑CATCGTGCCGAAGGCTGCAGCCGT‑CAGCAGATGGCTGGTCAG‑TTCCATCTGGCAGACC‑ATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTAT  >  NC_000913/4502513‑4503010

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