Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,658	A→G	18.6%	pseudogene (946/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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,658	0	A	G	18.6%	112.4 / 6.0	65	pseudogene (946/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 A (36/14);  new base G (2/11);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 3.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTC  >  NC_000913/3862423‑3862894                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGATGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑Gt                                                                                                                                                                                                                                                         >  1:484273/1‑239 (MQ=255)        tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TTAA‑GATGC‑gggg                                                                                                                                                                                                                                                 >  2:853302/1‑238 (MQ=255)        tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAG‑GGTGC‑GGGa                                                                                                                                                                                                                                                 >  1:1099932/1‑239 (MQ=255)            gATTTGGGGTAGCATCATATCCGCGTGGTTGCTGAACAACGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCGCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GTTGC‑GGGAATC‑a                                                                                                                                                                                                                                            <  2:480124/239‑1 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAACGAGATACATTACGGTCGACATTGCGACCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGAAAATTCCAGCGCTCCGGTAACTCCCACCAGCATGGC‑GG‑T‑TAAGTGTTGC‑gggcaccg                                                                                                                                                                                                                                             >  2:831559/1‑234 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGC‑GGGGATC‑a                                                                                                                                                                                                                                            >  1:10248/1‑238 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGC‑GGGGATC‑a                                                                                                                                                                                                                                            >  1:584744/1‑238 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAG‑GGTTC‑GGGCATC‑a                                                                                                                                                                                                                                            >  2:1099592/1‑238 (MQ=255)              tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TCAA‑GGTGC‑GGGAAAC‑Agg                                                                                                                                                                                                                                          >  1:671438/1‑238 (MQ=255)              tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GTTGC‑GGGAATC‑AGc                                                                                                                                                                                                                                          <  1:893081/239‑1 (MQ=255)              tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAG‑GGTGC‑ggggaacac                                                                                                                                                                                                                                            >  1:1168003/1‑232 (MQ=255)                      aGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGAGTACCGCAAACAGCAACGGGGAAATGAAAAGGAAGGCAAATTCCAGAGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAA‑GGTGG‑TGGGATT‑AGC‑AAGtc‑cc                                                                                                                                                                                                                                 >  2:29728/1‑237 (MQ=255)                             tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGGTCGGGAATTCCCACCAGCATGGC‑GG‑T‑GTGG‑GTTGC‑GGGCATC‑AAC‑CAGcc‑ccgacccttt                                                                                                                                                                                                                         >  2:929592/1‑230 (MQ=255)                                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAG‑GTTTC‑GGGAATC‑AGA‑AAGCC‑CGCCC‑CCTTT‑A‑ccc                                                                                                                                                                                                                  >  1:143025/1‑237 (MQ=255)                                   cgTGGTTGCTGAACATCGGGATACAGTTTTGCGGCAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGAGTACCGCAAACAGCAACGGTGAAATGAAACGGAAGGTAAATTCCAGCGGTGCGGTAATTCCAACAAGCATGGC‑GG‑T‑TGAA‑GTTGT‑GGGAATA‑AGC‑AAGCCACGTCA‑CTTTC‑A‑ccgg                                                                                                                                                                                                                 >  2:326370/1‑237 (MQ=255)                                       gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGC‑GGGGGTC‑AGC‑AAACC‑CCCCC‑CCTTT‑c‑cccccgttt                                                                                                                                                                                                            >  2:672052/1‑233 (MQ=255)                                          gCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GG‑T‑CACT‑AGTGA‑GGGAGCC‑CGT‑CACCC‑CCCCC‑CTCTT‑A‑CCCGATTTTct                                                                                                                                                                                                          >  2:779348/1‑237 (MQ=255)                                             gAACAGCGGGACCCAGTTTTGCGGTAAAACCTAGTCAATCAGCCCCCCGCCCATGTTGCCCACCACACCAAAGAGATACATTAGGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGAAAACAGCAACGGTGAAATGAACAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGG‑GG‑T‑TAAG‑GTTGC‑GGGAATC‑AGA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTg                                                                                                                                                                                                      <  1:1024951/239‑1 (MQ=255)                                                          cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGC‑GGGGATC‑AGC‑AAGAC‑CCCCC‑CTTTT‑T‑CCCGATTTTTCCGTGGGGCCGGGAAGTa                                                                                                                                                                                         >  1:414044/1‑239 (MQ=255)                                                               ttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TTAG‑GGTGG‑GGGGATC‑AGC‑AAGCC‑CCCCA‑CTTTT‑A‑ACCGATTTTCCGGTGGGGCGGGGAAGGACacc                                                                                                                                                                                     >  1:1183473/1‑236 (MQ=255)                                                                  cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGC‑GGGGAAC‑AGC‑AAACC‑CCCCC‑ATTTT‑A‑CCCGATTTTTCGGGGCGGGAGGGAAGGACATCG‑Ctc                                                                                                                                                                                >  1:725604/1‑238 (MQ=255)                                                                       aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGC‑GGGAATC‑AGC‑CAGCC‑CGCCA‑CATTT‑A‑CCCGGTTTTTCGGGGCGGCAGTGAAAGAAAACG‑CC‑AAAGaa                                                                                                                                                                          >  1:1179920/1‑239 (MQ=255)                                                                                   tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGG‑GGGAATC‑AACAAACCC‑CCCCC‑CTTTT‑A‑CCC‑ATTTTTCGGGGCGGCAGGGAAGGACAACC‑CT‑CAAAAAAAACCCACC‑G‑c                                                                                                                                                             >  1:804903/1‑238 (MQ=255)                                                                                   tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAA‑GTTGC‑GGGGATC‑AGC‑AAAAC‑CGCCC‑CTTTT‑A‑CCCGATTTTTCGGTGGGGCAGGGAAGGAAATCG‑CT‑AAAAAAAAGCCCACC‑G‑c                                                                                                                                                             >  1:1100845/1‑238 (MQ=255)                                                                                    cAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGAT‑TAAG‑G‑TGC‑GGGGATC‑AGC‑AAGGC‑CGCCA‑CCTTT‑A‑TCCGATTTTTCGGTGCGGCCGGGGAGGACAATC‑CT‑AAAGAAATGCCCACG‑G‑Ca                                                                                                                                                            >  1:730278/1‑238 (MQ=255)                                                                                            cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGG‑GGGGAAC‑AAC‑AAAAC‑CCCCC‑CTTTTAA‑CCCCTTTTTCCGG‑GCGGCAGGGAAGGACATCC‑CC‑AAAAAAAACCCCCCGCG‑CACCAAAgtt                                                                                                                                                    >  2:578228/1‑237 (MQ=255)                                                                                                     tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TAGG‑GTTGG‑GGGAA‑C‑AGC‑AAGGC‑CCCCC‑CCTTT‑T‑CCCGGTTTTCCCGTGCGGGCGTGAAATACCTCC‑CT‑AAAGAAAAGCCCCCA‑G‑CACCCAAGATTTTTTAGtt                                                                                                                                           >  1:992772/1‑238 (MQ=255)                                                                                                     tGCCCACAACACAAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTAAAATGAACAGGAAGGAAAATTCCAGCGGTTCGGTAATCCCCACCAGCAGGGC‑GG‑T‑TAAG‑GTTGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑a                                                                                                                                         <  1:672052/239‑1 (MQ=255)                                                                                                      gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TAAG‑GTT‑C‑GGGAATC‑AGC‑AAACC‑CCCCC‑CTTTT‑A‑CCCCATTTTTCCGGGCGGCCGTGAAATAAAACG‑CG‑AAAGAAAAACCCCCC‑G‑CCCCAAAGATTTTTGAGttt                                                                                                                                          >  1:1133954/1‑237 (MQ=255)                                                                                                      gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGCGGGGAATA‑GCC‑AACCC‑CCCCC‑CTTTT‑C‑CCCG‑TTTTTCGGTGGGGCGGGGAAGGACATCG‑CT‑AAAAAAAAGCCCACC‑G‑CCCCCAAAATTTTTTGGTT‑a                                                                                                                                         >  2:494058/1‑238 (MQ=255)                                                                                                                      agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TTAG‑GGTGC‑GGGGATC‑AGC‑AAGCC‑CCCCC‑CCTTT‑A‑CCCCATTTTCCGGTGCGGCAGGGTAGAAAAACG‑CT‑AAAAAAAAGCCCCCC‑G‑CACCCAAGAGTTTTTAGTT‑TCCGTGCCGGGG‑aaaaa                                                                                                                        >  2:1118046/1‑237 (MQ=255)                                                                                                                        aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TTAA‑GTTGC‑GGGGATC‑AGC‑AAAGC‑CCCCA‑CCTTT‑T‑CCCCATTTTTCCGGGCGGCAGGTGAGTAAAACGCCC‑AAAAAAAACCCCACG‑C‑ACCCAAAAATTTTTTGGTA‑CCCTTGCAGGGA‑AAAAcc                                                                                                                       >  1:692032/1‑238 (MQ=255)                                                                                                                         tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTGC‑GGGGATC‑AGC‑AAACC‑CGCCA‑CTTTT‑T‑CCCGCTTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAAAAATGCCCCCG‑G‑CACCAAAGATTTTTTAGTT‑ACCGTGCAGGGC‑AAAACCG‑c                                                                                                                    >  1:1083222/1‑238 (MQ=255)                                                                                                                              tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑T‑TAAG‑G‑TGT‑GGGGACC‑AGC‑CAAAC‑CGCCA‑CCTTT‑A‑TCCGATTTTTCCGGGCGGCAGGGAAGTAAAACC‑CG‑CCAGAAAAACCCACG‑GCCCCCAAGGATTTTTGGGTT‑‑TCGGGCGGGGC‑AAAAcaccctctc                                                                                                                >  2:1164301/1‑230 (MQ=255)                                                                                                                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAG‑GATGC‑GGGAATC‑AGC‑CAAAC‑CGCCC‑CTTTT‑ACCCCGATTTT‑CGGGGGGGCCGTGGAGGACAACG‑CT‑CAAAAAAAACCCACG‑G‑CCACAAAGGATTTTTAGTT‑TACGCGCCGGGC‑AAAAAAG‑CCCTCCcgg                                                                                                            >  2:768516/1‑239 (MQ=255)                                                                                                                                 ggTAGACATGGAGGCCGCCAGCACGGCGTGTCCCGAAAACAGCAACGGTAAAATGAACGGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCGGCGTGGG‑GG‑T‑TAGG‑GTGGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCggg                                                                                                            <  2:1100845/238‑1 (MQ=255)                                                                                                                                  gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAA‑GTTGC‑GGGGAAC‑AGC‑CAACC‑CCCCA‑CCTTT‑T‑CCCCCTTTTTCCGTGCGGCAGGGGAGAACACAC‑CT‑CAAAAAAAGCCCCCG‑G‑CACCAAAACATTTTGAGTT‑AACGTGCAGGGC‑AAAAACC‑CCTTCCcggt                                                                                                           >  2:424778/1‑237 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAG‑GGTTC‑GGGAATC‑AGC‑AAACC‑CCCCC‑CTTTT‑T‑CCCGATTTTTCGGTGCGGCAGGGAAATACATCG‑CT‑AAAGAAAAGCCCACG‑G‑CCCCAAAGATTTTTTAGTT‑ACCGTGCAGGGG‑AAAAACC‑CCCTCCGGGaaaaa                                                                                                       >  1:65011/1‑239 (MQ=255)                                                                                                                                                 gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TTAG‑GGTGG‑GGGAAGC‑AAC‑CAGCC‑CCCCACCTTTT‑A‑CCC‑ATTTTTCGGGGCGGCGGGGAAGTACACCG‑CC‑AAAGAAATTACCCCG‑G‑GCCCAAAGAATTTTTTTGT‑ATCGGGCCGGGCGAAAACGA‑CCTTCCGGG‑ACAAAGAttttcaagg                                                                                           >  2:15665/1‑233 (MQ=255)                                                                                                                                                 gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAA‑GGTTC‑GGGGATCAAGA‑AAGCC‑C‑CCC‑CTTTT‑T‑ACCGATTTTTCGGTGCGGCGGGGAAGTAAATAG‑CT‑TAAAAAAATCCCCCG‑G‑CCCCCCAGAATTTTTAGTTAACCGGCCCGGGA‑AAAACGC‑CCTCCCGGGaaaaaggtctctacg                                                                                             >  2:573412/1‑228 (MQ=255)                                                                                                                                                 gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAA‑GGTGC‑GGGGAAC‑AAG‑AAGGC‑CCCCA‑CCTTT‑T‑CCCGACTTTTCTGGGTGGCAGTGGAGGACAACGACT‑AAA‑AAAAACCCCCC‑G‑CCCCCAAGAAGTTTTTGTT‑ACCGCGTCGGGG‑ACAAACG‑CCCTCCCGGAACAACAACTTCAGCg                                                                                            >  2:21287/1‑238 (MQ=255)                                                                                                                                                           cgggtcccGAAAACAGAAACGGTGAAATGAACAGGAGGGTAAATCCCAGCGGTGCGGAAATCCCCACCAGCAGGGC‑GG‑T‑TAGG‑GTGGG‑GGGAATA‑AGC‑AGGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCt                                                                                  <  1:1168160/232‑1 (MQ=255)                                                                                                                                                           cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAG‑GGTGG‑GGGAAAC‑ACG‑AAACC‑CCCCC‑CTTTT‑T‑CCCCGTTTTTCCGTGCGGCAGTGAAGTAACACG‑CTAAAAAAAAGGCCCCCG‑C‑CCCCAAAG‑TTTTTTAGTT‑ACCGGGCAGGGG‑AAAAACGCCCTTCCGGG‑ACAACGATTTTCGCGCGTCGGGGCt                                                                                  >  2:640789/1‑238 (MQ=255)                                                                                                                                                                 cccaaaaCAGCAACGGTAAAAGGAACAGAAAGGAAAATTCCAGCGGTCCGTTAATCCCCACCAGCATGGC‑GG‑T‑TAGG‑GTGGC‑GGAAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                           <  1:599170/235‑1 (MQ=255)                                                                                                                                                                 ccGCAAACACACAGGGTAAAATAAACAGAAAGGAAAAGTCGAGCGGTTCGGTAACTCCCCCCGGCCGGGC‑GG‑T‑AAGG‑GGTGG‑GGAAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                           <  2:1168003/239‑1 (MQ=255)                                                                                                                                                                           cAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GGTTC‑GGGAATC‑AAA‑GAGCC‑CGGCA‑CTTTT‑A‑CCCGGTTTTCCGGTGCGGCAGGGGAGGACAACC‑CT‑CAAGAAATGCCCCCG‑G‑CCCCCAAGAATTTTGTGGT‑ACCGTGCCGGGC‑CAAAACG‑CCCTCCCGGGACAAAGATTTTAGCGGGTCGGCGCGCACACtagaccacctg                                                                  >  2:615608/1‑228 (MQ=255)                                                                                                                                                                            aggggTAAAATAAAGAGAAAGGAAAATTCCAGCGGTCGGTTAATCCCCACCAGCATGGC‑GG‑T‑TAAG‑GTGGG‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                <  1:853302/236‑1 (MQ=255)                                                                                                                                                                                gTGAAATGAACAGAAAGGAAAATTCCAGCGTTTCGGTAATCCCCACCAGCAGGGC‑GG‑T‑TAAG‑GTTGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                             <  1:1118046/238‑1 (MQ=255)                                                                                                                                                                                  aaaaTGAACAGGAAGGAAAATTCCGGCGTGCCGGAAATCCCCACCAGCATGGT‑GT‑T‑TAGG‑GTGGC‑GGGAACC‑AAC‑AACCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                           <  2:804903/237‑1 (MQ=255)                                                                                                                                                                                    aaagAACAGGAAGGTAAATCCCGGCGGTCCGATAATCCCCACCAACAGGGC‑TG‑T‑TAAG‑GTTGG‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                        <  1:15665/236‑1 (MQ=255)                                                                                                                                                                                    aaTGAAGAGAAGGGTAAATCCCGGGGGTTGGGTAATCCCCCCCGGCAT‑GC‑GG‑TGTAAG‑TGGGC‑GGAAACC‑AGA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                        <  1:745465/239‑1 (MQ=255)                                                                                                                                                                                    aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCCCCAGCAGGGC‑GG‑T‑TAAG‑GTGGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                        <  1:471051/239‑1 (MQ=255)                                                                                                                                                                                      taaaCAGGAAGGAAAATCCCAGCGGTTCGATAACCCCCACCAGCATGGC‑GT‑T‑TAAG‑GTGGG‑GGGAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                      <  1:1099592/237‑1 (MQ=255)                                                                                                                                                                                      tGAACAGCAAGGTAAATACCAGCGGTTCGTTAATTCCCACCAGCATGGC‑GG‑T‑TTAG‑GTTAC‑GGGAATC‑AAC‑AAGCC‑CGCCA‑CCTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTGCATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                      <  1:876063/239‑1 (MQ=255)                                                                                                                                                                                      tGAACAGAAAGGAAAATCCCAGCGTGGCGTTAATCCCCCCCAGCTTGTG‑GT‑T‑TAGG‑GTGGC‑GGGAATC‑AGC‑AAGGC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                      <  1:929592/239‑1 (MQ=255)                                                                                                                                                                                      agAACAGTAAGTTAAACCCCAGTGGTGTGATATCCCCCCCCAGCATGGC‑GT‑T‑TAGG‑GGTGC‑GGAAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                      <  2:725604/238‑1 (MQ=255)                                                                                                                                                                                           aGGAAGTTAAATTCCGGCGGTTTGGTAATCCCCACCAGGATGGC‑GT‑T‑TGTG‑GTTGC‑GGAAATC‑ATC‑ATGAC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGTGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                 <  1:144653/239‑1 (MQ=255)                                                                                                                                                                                            gaaaGGTAAATCCCGGCGGTCCGTAAACCCCCCCCAGCGTGGG‑GT‑T‑TAGG‑GTGGC‑GGAAACC‑AAA‑AACCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCGCG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                 <  2:16444/236‑1 (MQ=255)                                                                                                                                                                                                       tCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGGGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCGCG‑G‑CACCAAAGATTTTTTAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGGACAAAGATTTTAGCGGCTCGGGGCCCAAAATGAACTCTTGCAGATGCTGCGCCCAGGACATCTGAATgc                                      >  1:16444/1‑238 (MQ=255)                                                                                                                                                                                                                      atccccACCAGCATGGC‑GG‑T‑TAAG‑GTTGC‑GGGAACA‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                      <  1:615608/236‑1 (MQ=255)                                                                                                                                                                                                                             cgcgaaTGGC‑GT‑T‑AAGG‑GTCGC‑GGAAATA‑AGA‑ACGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa               <  1:327359/234‑1 (MQ=255)                                                                                                                                                                                                                             ccacaagggC‑GG‑T‑AAGG‑GGTGC‑GGAAATC‑AGA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa               <  1:549684/232‑1 (MQ=255)                                                                                                                                                                                                                             ccaGCTGGGC‑GG‑T‑TAGG‑GTGGC‑GGGAACC‑AGC‑ACGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa               <  2:992772/239‑1 (MQ=255)                                                                                                                                                                                                                                      c‑GG‑T‑TAAG‑GTTGG‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa       <  1:1216924/238‑1 (MQ=255)                                                                                                                                                                                                                                             tAAG‑GTTGG‑GGAATTC‑ACA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCATTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc  <  1:520416/239‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAG‑GTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTC  >  NC_000913/3862423‑3862894

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