Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,660	T→G	41.3%	pseudogene (943/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,660	0	T	G	41.3%	29.6 / 18.6	59	pseudogene (943/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 (9/23);  new base G (20/3);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 2.36e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCG  >  NC_000913/3862422‑3862895                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                 aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGt                                                                                                                                                                                                                                                                  >  1:241543/1‑239 (MQ=255)   cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑GTA‑AGGGTg                                                                                                                                                                                                                                                                >  1:1150412/1‑239 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATCa                                                                                                                                                                                                                                                       >  1:868914/1‑238 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGAAGCa                                                                                                                                                                                                                                                       >  1:922972/1‑238 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAAAAGGAAGGGAAATTCCAGCGGTTCGGTAA‑TTACCACCAGCATGGC‑TG‑TGT‑AGGtccggcttttttt                                                                                                                                                                                                                                                      >  2:657012/1‑227 (MQ=255)                        catcctATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCCAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCCCCAGCATGGC‑GG‑GTT‑AGGTTACGGGGAtgcgcacaccccac                                                                                                                                                                                                                                            >  2:790407/1‑225 (MQ=255)                           catATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGa                                                                                                                                                                                                                                                           <  1:678741/220‑1 (MQ=255)                                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑GTA‑AAGGTGCGGGAATCAG‑CA‑AACCCGCCCC‑C‑TTTTAccc                                                                                                                                                                                                                             >  2:458316/1‑238 (MQ=255)                                   cgTGGATTCTGACCATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGCGTACC‑GCAAACAGCAAC‑GGTGAAATGAACCGGAAGGGAAATTCCAGCGGTTCGGTAC‑TTCCCACCAGCCTGGC‑GG‑TTA‑AGGTTGCGGGAGGCAG‑CA‑AGACCG‑ccccctgttcccca                                                                                                                                                                                                                             >  2:257015/1‑227 (MQ=255)                                     tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATTAA‑CA‑AGGCCG‑CCC‑C‑CTTTTcccagtt                                                                                                                                                                                                                         >  1:20325/1‑235 (MQ=255)                                                           aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCCAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGGGGTTCGGTAA‑TTCCCACCAGCATGGCGGG‑TTA‑GGGTTGGGGGAA‑CAG‑AC‑AGACCC‑CCA‑C‑TTTTACACCA‑GTTTCCGGGGCgggcgggaagcac                                                                                                                                                                                                  >  1:680231/1‑228 (MQ=255)                                                                  cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑GTA‑AAGGTGCGGGAAACAG‑GC‑AGCCCC‑CCA‑C‑TTTTTACCGA‑TTTTTCCGGGCGGC‑A‑GTGGA‑GGACAT‑ccgtc                                                                                                                                                                                       >  1:692513/1‑235 (MQ=255)                                                                            taataacacacgccTCCGCCCACTTAGCCAACCACACAAAAGAACTACAGTACGGTCGCAAGTGAGGCCGCCGGCACGGCGTGTCCC‑GAAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATCCCAGCGTTTCGGTAC‑TCCCCTCCAGCTGGGG‑GG‑TTA‑GGGTTGGGGAAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑Gcc                                                                                                                                                                            <  2:1058402/227‑1 (MQ=255)                                                                                           ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAAAAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGAATGGC‑GG‑GCA‑AAGGTGCGGGGATCAA‑CA‑ACCCCG‑CCC‑C‑CTTTACCCCC‑GTTTCCCGGGTGGC‑A‑AAGGA‑GTAACC‑CCCTAAAGAAATACCCC‑CCGGC‑C‑CCa‑aaaa                                                                                                                                                          >  2:613355/1‑237 (MQ=255)                                                                                                     tGCCCACCACACCAAAGATATACATGACGGTCGACACGGACGCCGCCACAGCGGCGTTTACC‑GCAAACAGCAAC‑GGTGAAATGAAGAGAAAGAAAAATGCCGGCGGTTCGGTAA‑TCCCCCCCAGCATGGG‑GG‑TTA‑GGGTTGCGGAACTCAG‑AA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTa                                                                                                                                               <  2:1132705/239‑1 (MQ=255)                                                                                                      gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGAACAG‑CA‑AGCCCG‑CCA‑C‑CTTTTCCCGA‑TTTTTCGGTGCGGC‑C‑GTGAA‑GGACAT‑CGCTAAAGAAAA‑GACC‑ACCGC‑A‑CCC‑AAGAATTTTTAGTTa                                                                                                                                               >  1:1132705/1‑238 (MQ=255)                                                                                                               accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AAGATTGGGGAATTAGAAA‑AGCCGC‑CCC‑C‑TTTAACCCGA‑TTTTCCGG‑GGGGA‑C‑GTGAA‑GAACAC‑ACCCAAAAAAAT‑ACCC‑ACGGC‑C‑CCC‑AAAAGTTTTTAGGTatccctccggg                                                                                                                                     >  2:1140427/1‑229 (MQ=255)                                                                                                                      agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTTCGGGGATCAG‑CC‑AGCCCC‑CCA‑C‑ATTTTCCCCA‑ATTTCCGGTGCGGG‑C‑GTGTA‑GACCATACCCCAAAAAAAT‑‑CAC‑ACGGA‑C‑CCA‑AAAATTTTTTAGCTATCCTGC‑ACGGCAAAAc                                                                                                                              >  2:605043/1‑238 (MQ=255)                                                                                                                      agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑CTT‑AGGGTGCGGGAATCAG‑CA‑AACCCC‑CCC‑CTTTTTACCCG‑‑TTTTTCGGGGCGGC‑A‑GGGAA‑GTAAAT‑CGCTTAAAAAAT‑GCCC‑ACCGC‑C‑CCC‑AAGATTTTTTAGTTACCCTGC‑AGGGCcaaaa                                                                                                                              >  1:503473/1‑237 (MQ=255)                                                                                                                        aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑GTA‑AGGGTGCGGGGAATAG‑CA‑AGCCCC‑CCA‑C‑TTTTTACCGC‑TTTTCCTGGGGGGG‑C‑GCGAA‑GTACAT‑CCCTAAAAAAAT‑ACCC‑ACCGC‑ACCCC‑ACGATTTTTG‑GTTTCCATGC‑AGGGCGAAAccc                                                                                                                            >  2:110230/1‑237 (MQ=255)                                                                                                                           cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GGTTTA‑AGGTT‑CGGGGATCAG‑CA‑AGCCCG‑CCC‑C‑TTTTTCCCCA‑TTTTTCCGGGCGGG‑A‑GTGGA‑GGACAT‑AGCTAAAGAAAA‑TCCC‑ACAGC‑C‑CCA‑AAAATTTTTTAGTTTCCGTGC‑AGGGCCAAACCG‑ccc                                                                                                                        >  1:1054337/1‑237 (MQ=255)                                                                                                                           cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGa                                                                                                                                                                                                                                                           <  1:349287/124‑1 (MQ=255)                                                                                                                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGTGGGGATCAC‑CA‑AGACCC‑CCA‑C‑TTTTACCCAC‑TTTTCCGGGTCGGCGA‑GGAAG‑GTACAC‑CCCAAAAAAAAT‑‑CAC‑ACGGC‑C‑CCACAAGATTTTTTA‑TTACCGTGC‑CGGGCCAAAcacccctcccgg                                                                                                                   >  2:344489/1‑228 (MQ=255)                                                                                                                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGAAGCCG‑CC‑AACCCC‑CCA‑C‑CTTTTCCCCG‑TTTTTCGGTGCGGC‑C‑GGGGA‑GGAAAA‑CGCCAAAAAAAT‑ATCC‑CCGGC‑C‑CCA‑AAGCTTTTTTTGTTGCCGTGC‑CGGGGAAAACCC‑CCTCTCC‑gg                                                                                                                 >  2:407310/1‑239 (MQ=255)                                                                                                                                    cGACATAGAGGCCGCCAGCACGG‑GTGTACCAGCAAACAGAAAG‑GGAGAAATGAACAGGAAGATAAATGCCAGCGGTTCGGTAA‑TTCCCCCAAGCAGGGC‑GG‑TTA‑AGGTGGCGGAAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGaaca                                                                                                            <  2:8122/239‑1 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTT‑AAGGTGCGGGGATCAG‑CA‑AGCCCC‑CCA‑C‑CTTTTCCCCA‑TTTTTCCGGGCGGC‑A‑GGGAA‑GTACAT‑CGCTAAAAAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGAATTTTTGGTTTCCGTGG‑CGGGGAAAAACG‑CCC‑TCC‑GGGgacaa                                                                                                           >  1:854681/1‑239 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTTCGGGGAACAG‑CA‑AGCCCC‑CCC‑C‑TTTTTCCCCA‑ATTTTCGGTGCGGC‑A‑GACAAGGAAAAC‑CCCTAAAAAAAT‑GCCC‑AC‑GC‑G‑CCA‑AAAGTTTTTTTGTTTCCGTGG‑CGGGGACAACCGCCCT‑TCG‑GGGaaca                                                                                                            >  2:985629/1‑239 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATCAG‑CA‑AGCCCG‑CCC‑CTTTTTACCCAT‑TTTTCCGG‑GCGGC‑A‑GGGAA‑GGACAT‑CGCTAAAGAAAT‑ACCC‑CCCGC‑A‑CCC‑AAGAATTTTTAGTTTACGGGC‑CAGGCAAAAACC‑CCC‑TCC‑GGGaaaaa                                                                                                           >  1:316317/1‑239 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑GTA‑AGGTTGCGGGAAAAAG‑CA‑AAACCG‑CCA‑C‑CTTTAACCGACTTTTCGGCGGGGGC‑A‑G‑GAA‑GGACAT‑CGCTAAAGAAAA‑GCCC‑CCGGG‑A‑CCC‑AAAAATTTTTAGATAACGGGC‑AGGGGAAAACAC‑CCC‑TTC‑GGGaaaaa                                                                                                           >  2:742988/1‑239 (MQ=255)                                                                                                                                          tGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGGATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCCGTGCGGC‑A‑GTGGA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCC‑AAAAATTTTTAGTTACCGTGC‑AGGGCGAAACCG‑CCT‑TCC‑GGGAACAACGAtt                                                                                                      >  1:647463/1‑239 (MQ=255)                                                                                                                                                 gccAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGAATCAA‑CA‑AACCCC‑CCC‑C‑CTTTACCCGA‑TTTTCCGGGGCGGC‑A‑GGGGA‑GGACAT‑CGCTTAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTTAGGTACCGTGC‑AGGGGAAAACCG‑CCC‑TCC‑GGGGACAAAGATTTCAGccg                                                                                               >  1:352671/1‑237 (MQ=255)                                                                                                                                                 gccAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑GTA‑AGGATGCGGGGAACAG‑CACAGCCGC‑CCC‑C‑TTTAACCGCG‑TTTTCGGGGGCGCC‑A‑GTGAG‑GACAAT‑CGCAAAAAAAAT‑‑CCC‑CCGGCAA‑CAA‑AAAATTTTTAGGTTACGTGGC‑AGGGCAAAACC‑‑CCT‑TCCGGGGaacaacatttttgcg                                                                                                 >  2:461480/1‑229 (MQ=255)                                                                                                                                                  ccAGCACGGCGTGTAAC‑GCCAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGA‑GGCTTA‑AGGTCCGGGGAGTCAG‑‑A‑AGGCCG‑CCC‑C‑ATTTCACTGA‑GTTTCCGGGGGCGC‑C‑GTGAA‑GGCCCA‑CCCTCAAGACAT‑CCCCAACGCC‑C‑CCA‑AACTTCTTTGCG‑TACCATGT‑AGGGCAAAAACG‑CCTCTCC‑CGGaccaacatttttactgg                                                                                               >  2:836490/1‑228 (MQ=255)                                                                                                                                                    ggCAGGGCGAGTGCC‑GAAAACACCACA‑GGGAAAATGAACAGAAAGGTAACGTCCAGCGGTTGAATAA‑CTCCCACCAGCGTGGA‑GT‑TTA‑AGGTGGGGGGAACCAG‑AA‑CGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTc                                                                                            <  2:854681/238‑1 (MQ=255)                                                                                                                                                    acggcGGCGTGTACC‑‑TAAACAACAACGGGAGAAATGAAGAAGAAGGAAAATTCGAGCGTTGCGATAT‑TCCCCCCCAGCA‑GGC‑GG‑GTAGGGGGGGCGGAAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTc                                                                                            <  2:352671/235‑1 (MQ=255)                                                                                                                                                               tACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGAACAG‑CA‑AGCCCC‑CCACC‑TTTTACCC‑A‑TTTTTCGGGGCGGG‑AGGGGAA‑GAAAAC‑CCCAAAAAAAAG‑CCCC‑CCGCC‑C‑CCA‑AAGTTTTTTAGGTACCCGGGA‑GGGGCAAAACC‑‑CCC‑TCC‑GGGAACAAACACTTTCGCGGCCCGGC‑G‑CTCaaaa                                                                               >  2:38055/1‑238 (MQ=255)                                                                                                                                                               tACC‑GAAAACAGCAAG‑GATGAAATGAGGAGGAAGAAAAATTCCGGCGGTGCGTAAA‑GCCCCCCCAGCATGGC‑GT‑TGA‑GGGTCGCGGAAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCaaa                                                                                <  2:647463/238‑1 (MQ=255)                                                                                                                                                                    gaaaaCAGAAAC‑GGGAAAAGAAACAGGAAGGTAAATCCCAGCGGTTCGGTAA‑TCCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGaa                                                                           <  1:461480/237‑1 (MQ=255)                                                                                                                                                                     cAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGATGCGGGGATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAG‑GGACAT‑CGCTCAAAGAAT‑GGCC‑ACGGC‑A‑CCC‑AAGAGTTTTGAGATATCGGGC‑AGGGGAAAAACG‑CCG‑TCC‑GGGGACAAAGATTTCACTGGGTCGGC‑G‑CCCACACCGaa                                                                           >  2:141589/1‑238 (MQ=255)                                                                                                                                                                      aaaCAGCAGC‑GGTAAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TCCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATAAA‑CA‑ACCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAAct                                                                         <  1:613355/239‑1 (MQ=255)                                                                                                                                                                        acaaaaaC‑GGTAAAATGAACGGGAAGGTAAATCCCAGCGGTTCGTTAA‑TCCCCACCAGCAGGGC‑GT‑TTA‑AGTTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAActct                                                                       <  1:458316/234‑1 (MQ=255)                                                                                                                                                                          aGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑GTA‑AGGGTGCGGGAAACAG‑CA‑AGCCCG‑CCA‑C‑CTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GGACAA‑CGCTAAAAAAAT‑GCCC‑ACGGC‑C‑CCA‑AAGAATTTTTTGTTAACGGGC‑CGGGCAAAACCG‑CCC‑TCC‑GGGGAAAACGATTTCAGCGGCTCgcgacccaaaacttacactt                                                                        >  2:183424/1‑219 (MQ=255)                                                                                                                                                                               c‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTTAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGAt                                                                >  1:625382/1‑239 (MQ=255)                                                                                                                                                                                    gAAATAACCAGCAAGGTACATCCCAGTGGTGGTATAC‑TCCCCACGAGCAGGGG‑GG‑TTG‑AGGTTGGGGGAAGCAC‑CA‑GTCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGAtgct                                                             <  2:922972/238‑1 (MQ=255)                                                                                                                                                                                     aaaTGAACAGGAAGGTAAATTCCAGCGGTACGGTAA‑TTCCCACCAGCATGGC‑GG‑TCA‑AGGTTGCTGGGATCAG‑CC‑AGCCCC‑CCC‑C‑TTTTTTCCGA‑TTTTCCGGGGCTGC‑C‑GTGCA‑GTGCAC‑CCCTAACGCAAA‑CCCC‑CCAGC‑A‑CCC‑AAGTATTTTGAATTTATGGAA‑AGGGCCAAACCA‑CCC‑TCC‑CCGGCAAACGATTTCAAGCGCTCGGCGG‑CTCAAACCGAAccctcctcgaaccgc                                                            >  2:829189/1‑228 (MQ=255)                                                                                                                                                                                        tGAACAGGAAGATAAATTCAAGCGGTTCGATAA‑TCCCCACCAGCATGGC‑GT‑TTA‑AGGTTGCGGAAATCAG‑CA‑GGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  1:141589/239‑1 (MQ=255)                                                                                                                                                                                        tGAACAGAAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCAGGGC‑GG‑TTA‑AGGTGGGGGGAACCAG‑AA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  1:110230/239‑1 (MQ=255)                                                                                                                                                                                           aacagaagtgcaaatcaggCGGTTCGGGAACTCCCCCCCAGCAGGGC‑GG‑TTG‑GGGGTGCGGAAACCAG‑CA‑CGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  2:20325/219‑1 (MQ=255)                                                                                                                                                                                             agaaaGGTAAATCCCGGCGTTTGGGTAA‑TCCCCCCCAGCATGGG‑GG‑TTA‑GGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                   <  1:481540/236‑1 (MQ=255)                                                                                                                                                                                             aGGAAGGTAAATCCCAGCGGTTCGGTAA‑TCCCCCCCAGCAGGGC‑GG‑TTA‑GGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                   <  1:1017159/239‑1 (MQ=255)                                                                                                                                                                                                    tAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCCA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAA‑GCCC‑CCGGC‑A‑CCA‑AAGATTTTTTAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCCGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGa                                            >  1:440178/1‑239 (MQ=255)                                                                                                                                                                                                           cAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AAGATGCGGGAAACAG‑CC‑AGCCCG‑CCA‑C‑TTTTTCCCGA‑TTTTCCGGTGCGGC‑A‑GGGAA‑GGACAA‑CGCTAAAAAAAT‑ACCC‑CCCGC‑C‑CCC‑AAGATTTTTTAGTTACCGTGCGAGGGAAAAACAC‑CCT‑CCC‑GGGAA‑AACGATTTCAACGGCTCGGG‑GCCTCAAAATAAAGCCTT‑CAGAtgcagccgcccgacaaaccgtatcccc                                      >  2:255437/1‑217 (MQ=255)                                                                                                                                                                                                                          tgtaaaagaggaggggggg‑TTA‑AGGTTGCGGAAATCAA‑AA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                        <  1:767903/222‑1 (MQ=255)                                                                                                                                                                                                                              ccccaGCAGGCC‑GG‑TTA‑AGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  1:601558/237‑1 (MQ=255)                                                                                                                                                                                                                              cacaaacATGCC‑GG‑TTA‑GGTTAGCGGGAACCAA‑CA‑CGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  1:829189/233‑1 (MQ=255)                                                                                                                                                                                                                                ccagatgggC‑GG‑TTA‑AGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                 <  1:877122/232‑1 (MQ=255)                                                                                                                                                                                                                                ccaGCAGGGC‑GG‑TTA‑AGGGGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                 <  1:407310/239‑1 (MQ=255)                                                                                                                                                                                                                                      tGGC‑GG‑TAA‑ATCTTGCGGAAACCGC‑CA‑ATCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGAt            <  1:239702/238‑1 (MQ=255)                                                                                                                                                                                                                                            g‑TTA‑AGGTGGCGGAAAACAG‑CA‑GGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACt       <  1:742988/238‑1 (MQ=255)                                                                                                                                                                                                                                               tA‑AGGTTGCGGGAATAAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc    <  1:38055/239‑1 (MQ=255)                                                                                                                                                                                                                                                  aGGTGGCGAGAATCAC‑CA‑GCCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCg  <  1:257015/239‑1 (MQ=255)                                                                                                                                                                                                                                                     ttGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CCCAAACTGAACTCTTGCAGCTGCTGCGCCCCGTACATACGAATGCCGCCTTCAACAGCTGCCGTACCAAAg              >  2:767903/1‑224 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                 AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAA‑TTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGC‑AGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑G‑CTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCG  >  NC_000913/3862422‑3862895

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