Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,686	T→C	14.8%	pseudogene (918/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,686	0	T	C	14.8%	156.9 / 4.3	67	pseudogene (918/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 (24/31);  new base C (9/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 4.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTT  >  NC_000913/3862457‑3862920                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    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)                                                                                                                                                                                                                          aaccacaaaccc‑cGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACAGCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg                 <  1:661418/228‑1 (MQ=255)                                                                                                                                                                                                                           c‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCACTACATCTGAATGCCGCATTCAACAGCTGCCAGTCCAAAGATAAACTGGCCGTCGATGaatcggt             >  2:689575/1‑234 (MQ=255)                                                                                                                                                                                                                                   gCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTaaa        >  1:1081743/1‑238 (MQ=255)                                                                                                                                                                                                                                        gCCC‑CTTTT‑T‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGt   <  2:897927/239‑1 (MQ=255)                                                                                                                                                                                                                                         ccA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg                 <  2:1002676/224‑1 (MQ=255)                                                                                                                                                                                                                                         ccA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg                 >  1:1002676/1‑224 (MQ=255)                                                                                                                                                                                                                                         ccA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGtt  >  1:277966/1‑239 (MQ=255)                                                                                                                                                                                                                                          cA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGtt  <  2:414044/238‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    CGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑T‑TAAGGTTGC‑GGGAATC‑AGC‑AAGCC‑CGCCA‑CTTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCG‑CT‑AAAGAAATGCCCACG‑G‑CACCAAAGATTTTTGAGTT‑ACCGTGCAGGGC‑AAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTT  >  NC_000913/3862457‑3862920

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