Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,670	T→C	18.6%	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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,670	0	T	C	18.6%	80.5 / 7.4	66	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 (25/22);  new base C (2/9);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGG  >  NC_000913/3862434‑3862908                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            gATTTGGGGTAGCATCATATCCGCGTGGTTGCTGAACAACGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCGCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAATC‑a                                                                                                                                                                                                                                                          <  2:480124/239‑1 (MQ=255)  aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGC‑GGGGATC‑a                                                                                                                                                                                                                                                          >  1:584744/1‑238 (MQ=255)  aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGC‑GGGGATC‑a                                                                                                                                                                                                                                                          >  1:10248/1‑238 (MQ=255)  aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAGGGTTC‑GGGCATC‑a                                                                                                                                                                                                                                                          >  2:1099592/1‑238 (MQ=255)   tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TCAAGGTGC‑GGGAAAC‑Agg                                                                                                                                                                                                                                                        >  1:671438/1‑238 (MQ=255)   tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAATC‑AGc                                                                                                                                                                                                                                                        <  1:893081/239‑1 (MQ=255)           aGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGAGTACCGCAAACAGCAACGGGGAAATGAAAAGGAAGGCAAATTCCAGAGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAAGGTGG‑TGGGATT‑AGC‑AAGtc‑cc                                                                                                                                                                                                                                               >  2:29728/1‑237 (MQ=255)                  tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGGTCGGGAATTCCCACCAGCATGGC‑GG‑T‑GTGGGTTGC‑GGGCATC‑AAC‑CAGcc‑ccgacccttt                                                                                                                                                                                                                                       >  2:929592/1‑230 (MQ=255)                        cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAGGTTTC‑GGGAATC‑AGA‑AAGCC‑CGCCC‑CCTTT‑A‑ccc                                                                                                                                                                                                                                >  1:143025/1‑237 (MQ=255)                        cgTGGTTGCTGAACATCGGGATACAGTTTTGCGGCAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGAGTACCGCAAACAGCAACGGTGAAATGAAACGGAAGGTAAATTCCAGCGGTGCGGTAATTCCAACAAGCATGGC‑GG‑T‑TGAAGTTGT‑GGGAATA‑AGC‑AAGCCACGTCA‑CTTTC‑A‑ccgg                                                                                                                                                                                                                               >  2:326370/1‑237 (MQ=255)                            gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGC‑GGGGGTC‑AGC‑AAACC‑CCCCC‑CCTTT‑c‑cccccgttt                                                                                                                                                                                                                          >  2:672052/1‑233 (MQ=255)                               gCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GG‑T‑CACTAGTGA‑GGGAGCC‑CGT‑CACCC‑CCCCC‑CTCTT‑A‑CCCGATTTTct                                                                                                                                                                                                                        >  2:779348/1‑237 (MQ=255)                                  gAACAGCGGGACCCAGTTTTGCGGTAAAACCTAGTCAATCAGCCCCCCGCCCATGTTGCCCACCACACCAAAGAGATACATTAGGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGAAAACAGCAACGGTGAAATGAACAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGG‑GG‑T‑TAAGGTTGC‑GGGAATC‑AGA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTg                                                                                                                                                                                                                    <  1:1024951/239‑1 (MQ=255)                                               cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGC‑GGGGATC‑AGC‑AAGAC‑CCCCC‑CTTTT‑T‑CCCGATTTTTCCGTGGGGCCGGGAAGTa                                                                                                                                                                                                       >  1:414044/1‑239 (MQ=255)                                                    ttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TTAGGGTGG‑GGGGATC‑AGC‑AAGCC‑CCCCA‑CTTTT‑A‑ACCGATTTTCCGGTGGGGCGGGGAAGGACacc                                                                                                                                                                                                   >  1:1183473/1‑236 (MQ=255)                                                       cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGC‑GGGGAAC‑AGC‑AAACC‑CCCCC‑ATTTT‑A‑CCCGATTTTTCGGGGCGGGAGGGAAGGACATCG‑Ctc                                                                                                                                                                                              >  1:725604/1‑238 (MQ=255)                                                            aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGC‑GGGAATC‑AGC‑CAGCC‑CGCCA‑CATTT‑A‑CCCGGTTTTTCGGGGCGGCAGTGAAAGAAAACG‑CC‑AAAGaa                                                                                                                                                                                        >  1:1179920/1‑239 (MQ=255)                                                                        tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGG‑GGGAATC‑AACAAACCC‑CCCCC‑CTTTT‑A‑CCC‑ATTTTTCGGGGCGGCAGGGAAGGACAACC‑CT‑CAAAAAAAACCCACC‑G‑c                                                                                                                                                                           >  1:804903/1‑238 (MQ=255)                                                                        tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAAGTTGC‑GGGGATC‑AGC‑AAAAC‑CGCCC‑CTTTT‑A‑CCCGATTTTTCGGTGGGGCAGGGAAGGAAATCG‑CT‑AAAAAAAAGCCCACC‑G‑c                                                                                                                                                                           >  1:1100845/1‑238 (MQ=255)                                                                         cAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGAT‑TAAGG‑TGC‑GGGGATC‑AGC‑AAGGC‑CGCCA‑CCTTT‑A‑TCCGATTTTTCGGTGCGGCCGGGGAGGACAATC‑CT‑AAAGAAATGCCCACG‑G‑Ca                                                                                                                                                                          >  1:730278/1‑238 (MQ=255)                                                                                 cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGG‑GGGGAAC‑AAC‑AAAAC‑CCCCC‑CTTTTAA‑CCCCTTTTTCCGG‑GCGGCAGGGAAGGACATCC‑CC‑AAAAAAAACCCCCCGCG‑CACCAAAgtt                                                                                                                                                                  >  2:578228/1‑237 (MQ=255)                                                                                          tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TAGGGTTGG‑GGGAA‑C‑AGC‑AAGGC‑CCCCC‑CCTTT‑T‑CCCGGTTTTCCCGTGCGGGCGTGAAATACCTCC‑CT‑AAAGAAAAGCCCCCA‑G‑CACCCAAGATTTTTTAGtt                                                                                                                                                         >  1:992772/1‑238 (MQ=255)                                                                                          tGCCCACAACACAAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTAAAATGAACAGGAAGGAAAATTCCAGCGGTTCGGTAATCCCCACCAGCAGGGC‑GG‑T‑TAAGGTTGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑a                                                                                                                                                       <  1:672052/239‑1 (MQ=255)                                                                                           gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TAAGGTT‑C‑GGGAATC‑AGC‑AAACC‑CCCCC‑CTTTT‑A‑CCCCATTTTTCCGGGCGGCCGTGAAATAAAACG‑CG‑AAAGAAAAACCCCCC‑G‑CCCCAAAGATTTTTGAGttt                                                                                                                                                        >  1:1133954/1‑237 (MQ=255)                                                                                           gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGCGGGGAATA‑GCC‑AACCC‑CCCCC‑CTTTT‑C‑CCCG‑TTTTTCGGTGGGGCGGGGAAGGACATCG‑CT‑AAAAAAAAGCCCACC‑G‑CCCCCAAAATTTTTTGGTT‑a                                                                                                                                                       >  2:494058/1‑238 (MQ=255)                                                                                                           agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TTAGGGTGC‑GGGGATC‑AGC‑AAGCC‑CCCCC‑CCTTT‑A‑CCCCATTTTCCGGTGCGGCAGGGTAGAAAAACG‑CT‑AAAAAAAAGCCCCCC‑G‑CACCCAAGAGTTTTTAGTT‑TCCGTGCCGGGG‑aaaaa                                                                                                                                      >  2:1118046/1‑237 (MQ=255)                                                                                                             aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TTAAGTTGC‑GGGGATC‑AGC‑AAAGC‑CCCCA‑CCTTT‑T‑CCCCATTTTTCCGGGCGGCAGGTGAGTAAAACGCCC‑AAAAAAAACCCCACG‑C‑ACCCAAAAATTTTTTGGTA‑CCCTTGCAGGGA‑AAAAcc                                                                                                                                     >  1:692032/1‑238 (MQ=255)                                                                                                              tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTGC‑GGGGATC‑AGC‑AAACC‑CGCCA‑CTTTT‑T‑CCCGCTTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAAAAATGCCCCCG‑G‑CACCAAAGATTTTTTAGTT‑ACCGTGCAGGGC‑AAAACCG‑c                                                                                                                                  >  1:1083222/1‑238 (MQ=255)                                                                                                                   tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑T‑TAAGG‑TGT‑GGGGACC‑AGC‑CAAAC‑CGCCA‑CCTTT‑A‑TCCGATTTTTCCGGGCGGCAGGGAAGTAAAACC‑CG‑CCAGAAAAACCCACG‑GCCCCCAAGGATTTTTGGGTT‑‑TCGGGCGGGGC‑AAAAcaccctctc                                                                                                                              >  2:1164301/1‑230 (MQ=255)                                                                                                                     cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAGGATGC‑GGGAATC‑AGC‑CAAAC‑CGCCC‑CTTTT‑ACCCCGATTTT‑CGGGGGGGCCGTGGAGGACAACG‑CT‑CAAAAAAAACCCACG‑G‑CCACAAAGGATTTTTAGTT‑TACGCGCCGGGC‑AAAAAAG‑CCCTCCcgg                                                                                                                          >  2:768516/1‑239 (MQ=255)                                                                                                                      ggTAGACATGGAGGCCGCCAGCACGGCGTGTCCCGAAAACAGCAACGGTAAAATGAACGGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCGGCGTGGG‑GG‑T‑TAGGGTGGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCggg                                                                                                                          <  2:1100845/238‑1 (MQ=255)                                                                                                                       gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAAGTTGC‑GGGGAAC‑AGC‑CAACC‑CCCCA‑CCTTT‑T‑CCCCCTTTTTCCGTGCGGCAGGGGAGAACACAC‑CT‑CAAAAAAAGCCCCCG‑G‑CACCAAAACATTTTGAGTT‑AACGTGCAGGGC‑AAAAACC‑CCTTCCcggt                                                                                                                         >  2:424778/1‑237 (MQ=255)                                                                                                                          gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAGGGTTC‑GGGAATC‑AGC‑AAACC‑CCCCC‑CTTTT‑T‑CCCGATTTTTCGGTGCGGCAGGGAAATACATCG‑CT‑AAAGAAAAGCCCACG‑G‑CCCCAAAGATTTTTTAGTT‑ACCGTGCAGGGG‑AAAAACC‑CCCTCCGGGaaaaa                                                                                                                     >  1:65011/1‑239 (MQ=255)                                                                                                                                      gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TTAGGGTGG‑GGGAAGC‑AAC‑CAGCC‑CCCCACCTTTT‑A‑CCC‑ATTTTTCGGGGCGGCGGGGAAGTACACCG‑CC‑AAAGAAATTACCCCG‑G‑GCCCAAAGAATTTTTTTGT‑ATCGGGCCGGGCGAAAACGA‑CCTTCCGGG‑ACAAAGAttttcaagg                                                                                                         >  2:15665/1‑233 (MQ=255)                                                                                                                                      gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAAGGTTC‑GGGGATCAAGA‑AAGCC‑C‑CCC‑CTTTT‑T‑ACCGATTTTTCGGTGCGGCGGGGAAGTAAATAG‑CT‑TAAAAAAATCCCCCG‑G‑CCCCCCAGAATTTTTAGTTAACCGGCCCGGGA‑AAAACGC‑CCTCCCGGGaaaaaggtctctacg                                                                                                           >  2:573412/1‑228 (MQ=255)                                                                                                                                      gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAAGGTGC‑GGGGAAC‑AAG‑AAGGC‑CCCCA‑CCTTT‑T‑CCCGACTTTTCTGGGTGGCAGTGGAGGACAACGACT‑AAA‑AAAAACCCCCC‑G‑CCCCCAAGAAGTTTTTGTT‑ACCGCGTCGGGG‑ACAAACG‑CCCTCCCGGAACAACAACTTCAGCg                                                                                                          >  2:21287/1‑238 (MQ=255)                                                                                                                                                cgggtcccGAAAACAGAAACGGTGAAATGAACAGGAGGGTAAATCCCAGCGGTGCGGAAATCCCCACCAGCAGGGC‑GG‑T‑TAGGGTGGG‑GGGAATA‑AGC‑AGGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCt                                                                                                <  1:1168160/232‑1 (MQ=255)                                                                                                                                                cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑G‑TAAGGGTGG‑GGGAAAC‑ACG‑AAACC‑CCCCC‑CTTTT‑T‑CCCCGTTTTTCCGTGCGGCAGTGAAGTAACACG‑CTAAAAAAAAGGCCCCCG‑C‑CCCCAAAG‑TTTTTTAGTT‑ACCGGGCAGGGG‑AAAAACGCCCTTCCGGG‑ACAACGATTTTCGCGCGTCGGGGCt                                                                                                >  2:640789/1‑238 (MQ=255)                                                                                                                                                      cccaaaaCAGCAACGGTAAAAGGAACAGAAAGGAAAATTCCAGCGGTCCGTTAATCCCCACCAGCATGGC‑GG‑T‑TAGGGTGGC‑GGAAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                                         <  1:599170/235‑1 (MQ=255)                                                                                                                                                      ccGCAAACACACAGGGTAAAATAAACAGAAAGGAAAAGTCGAGCGGTTCGGTAACTCCCCCCGGCCGGGC‑GG‑T‑AAGGGGTGG‑GGAAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                                         <  2:1168003/239‑1 (MQ=255)                                                                                                                                                                cAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGGTTC‑GGGAATC‑AAA‑GAGCC‑CGGCA‑CTTTT‑A‑CCCGGTTTTCCGGTGCGGCAGGGGAGGACAACC‑CT‑CAAGAAATGCCCCCG‑G‑CCCCCAAGAATTTTGTGGT‑ACCGTGCCGGGC‑CAAAACG‑CCCTCCCGGGACAAAGATTTTAGCGGGTCGGCGCGCACACtagaccacctg                                                                                >  2:615608/1‑228 (MQ=255)                                                                                                                                                                 aggggTAAAATAAAGAGAAAGGAAAATTCCAGCGGTCGGTTAATCCCCACCAGCATGGC‑GG‑T‑TAAGGTGGG‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                              <  1:853302/236‑1 (MQ=255)                                                                                                                                                                     gTGAAATGAACAGAAAGGAAAATTCCAGCGTTTCGGTAATCCCCACCAGCAGGGC‑GG‑T‑TAAGGTTGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                           <  1:1118046/238‑1 (MQ=255)                                                                                                                                                                       aaaaTGAACAGGAAGGAAAATTCCGGCGTGCCGGAAATCCCCACCAGCATGGT‑GT‑T‑TAGGGTGGC‑GGGAACC‑AAC‑AACCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                         <  2:804903/237‑1 (MQ=255)                                                                                                                                                                         aaagAACAGGAAGGTAAATCCCGGCGGTCCGATAATCCCCACCAACAGGGC‑TG‑T‑TAAGGTTGG‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                      <  1:15665/236‑1 (MQ=255)                                                                                                                                                                         aaTGAAGAGAAGGGTAAATCCCGGGGGTTGGGTAATCCCCCCCGGCAT‑GC‑GG‑TGTAAGTGGGC‑GGAAACC‑AGA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                      <  1:745465/239‑1 (MQ=255)                                                                                                                                                                         aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCCCCAGCAGGGC‑GG‑T‑TAAGGTGGC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                      <  1:471051/239‑1 (MQ=255)                                                                                                                                                                           taaaCAGGAAGGAAAATCCCAGCGGTTCGATAACCCCCACCAGCATGGC‑GT‑T‑TAAGGTGGG‑GGGAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  1:1099592/237‑1 (MQ=255)                                                                                                                                                                           tGAACAGCAAGGTAAATACCAGCGGTTCGTTAATTCCCACCAGCATGGC‑GG‑T‑TTAGGTTAC‑GGGAATC‑AAC‑AAGCC‑CGCCA‑CCTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTGCATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  1:876063/239‑1 (MQ=255)                                                                                                                                                                           tGAACAGAAAGGAAAATCCCAGCGTGGCGTTAATCCCCCCCAGCTTGTG‑GT‑T‑TAGGGTGGC‑GGGAATC‑AGC‑AAGGC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  1:929592/239‑1 (MQ=255)                                                                                                                                                                           agAACAGTAAGTTAAACCCCAGTGGTGTGATATCCCCCCCCAGCATGGC‑GT‑T‑TAGGGGTGC‑GGAAACC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  2:725604/238‑1 (MQ=255)                                                                                                                                                                                aGGAAGTTAAATTCCGGCGGTTTGGTAATCCCCACCAGGATGGC‑GT‑T‑TGTGGTTGC‑GGAAATC‑ATC‑ATGAC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGTGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                               <  1:144653/239‑1 (MQ=255)                                                                                                                                                                                 gaaaGGTAAATCCCGGCGGTCCGTAAACCCCCCCCAGCGTGGG‑GT‑T‑TAGGGTGGC‑GGAAACC‑AAA‑AACCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCGCG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                               <  2:16444/236‑1 (MQ=255)                                                                                                                                                                                            tCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGGGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCGCG‑G‑CACCAAAGATTTTTTAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGGACAAAGATTTTAGCGGCTCGGGGCCCAAAATGAACTCTTGCAGATGCTGCGCCCAGGACATCTGAATgc                                                    >  1:16444/1‑238 (MQ=255)                                                                                                                                                                                                           atccccACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAACA‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                                    <  1:615608/236‑1 (MQ=255)                                                                                                                                                                                                                  cgcgaaTGGC‑GT‑T‑AAGGGTCGC‑GGAAATA‑AGA‑ACGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                             <  1:327359/234‑1 (MQ=255)                                                                                                                                                                                                                  ccacaagggC‑GG‑T‑AAGGGGTGC‑GGAAATC‑AGA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                             <  1:549684/232‑1 (MQ=255)                                                                                                                                                                                                                  ccaGCTGGGC‑GG‑T‑TAGGGTGGC‑GGGAACC‑AGC‑ACGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                             <  2:992772/239‑1 (MQ=255)                                                                                                                                                                                                                           c‑GG‑T‑TAAGGTTGG‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                     <  1:1216924/238‑1 (MQ=255)                                                                                                                                                                                                                              cgttatggtgggC‑GGTAATC‑AGC‑CACCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                     <  2:1179920/227‑1 (MQ=255)                                                                                                                                                                                                                                gtaaggtgggC‑GGAAATA‑ACC‑AACCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACt                   <  1:687554/229‑1 (MQ=255)                                                                                                                                                                                                                                  tAAGGTTGG‑GGAATTC‑ACA‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCATTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc                <  1:520416/239‑1 (MQ=255)                                                                                                                                                                                                                                   agggTTGC‑GGAAATC‑AGC‑AGGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCAGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTcc               <  1:743909/237‑1 (MQ=255)                                                                                                                                                                                                                                   agggTTGC‑GGAAATC‑AGC‑AACCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTcc               <  1:640789/237‑1 (MQ=255)                                                                                                                                                                                                                                      gTTCC‑GGAAATC‑AGC‑AAGCC‑CGCCA‑CCTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt             <  1:664855/238‑1 (MQ=255)                                                                                                                                                                                                                                          c‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGGACATATGAATGCCGCCTTCAACCGCTGCCGGACCAAAGATAACCTGTCCGTAgtcg        >  2:735657/1‑236 (MQ=255)                                                                                                                                                                                                                                                 tC‑AGC‑AAGAC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGTGCTCAAACTGAACTCTTGCAGATGATGCGCCCAGTACATTTGAAAGCCCCCTTCCACAGCTGCCGGACCAAAGCTACATTGGCCGTAGagtaagggg  >  2:775094/1‑231 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            GATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGG  >  NC_000913/3862434‑3862908

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: ‑