Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,661	T→G	42.2%	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,661	0	T	G	42.2%	25.1 / 17.6	52	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 (7/21);  new base G (19/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 6.71e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-02

AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCG  >  NC_000913/3862423‑3862896                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                               aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGt                                                                                                                                                                                                                                                                >  1:241543/1‑239 (MQ=255)   cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTA‑AGGGTg                                                                                                                                                                                                                                                              >  1:1150412/1‑239 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATCa                                                                                                                                                                                                                                                     >  1:868914/1‑238 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGAAGCa                                                                                                                                                                                                                                                     >  1:922972/1‑238 (MQ=255)                           catATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGa                                                                                                                                                                                                                                                         <  1:678741/220‑1 (MQ=255)                                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTA‑AAGGTGCGGGAATCAG‑CA‑AACCCGCCCC‑C‑TTTTAccc                                                                                                                                                                                                                           >  2:458316/1‑238 (MQ=255)                                     tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATTAA‑CA‑AGGCCG‑CCC‑C‑CTTTTcccagtt                                                                                                                                                                                                                       >  1:20325/1‑235 (MQ=255)                                                           aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCCAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGGGGTTCGGTAATTCCCACCAGCATGGCGGG‑TTA‑GGGTTGGGGGAA‑CAG‑AC‑AGACCC‑CCA‑C‑TTTTACACCA‑GTTTCCGGGGCgggcgggaagcac                                                                                                                                                                                                >  1:680231/1‑228 (MQ=255)                                                                  cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTA‑AAGGTGCGGGAAACAG‑GC‑AGCCCC‑CCA‑C‑TTTTTACCGA‑TTTTTCCGGGCGGC‑A‑GTGGA‑GGACAT‑ccgtc                                                                                                                                                                                     >  1:692513/1‑235 (MQ=255)                                                                                           ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAAAAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGAATGGC‑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‑GGTGAAATGAAGAGAAAGAAAAATGCCGGCGGTTCGGTAATCCCCCCCAGCATGGG‑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‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑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‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑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‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTTCGGGGATCAG‑CC‑AGCCCC‑CCA‑C‑ATTTTCCCCA‑ATTTCCGGTGCGGG‑C‑GTGTA‑GACCATACCCCAAAAAAAT‑‑CAC‑ACGGA‑C‑CCA‑AAAATTTTTTAGCTATCCTGCACGGCAAAAc                                                                                                                             >  2:605043/1‑238 (MQ=255)                                                                                                                      agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑CTT‑AGGGTGCGGGAATCAG‑CA‑AACCCC‑CCC‑CTTTTTACCCG‑‑TTTTTCGGGGCGGC‑A‑GGGAA‑GTAAAT‑CGCTTAAAAAAT‑GCCC‑ACCGC‑C‑CCC‑AAGATTTTTTAGTTACCCTGCAGGGCcaaaa                                                                                                                             >  1:503473/1‑237 (MQ=255)                                                                                                                        aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTA‑AGGGTGCGGGGAATAG‑CA‑AGCCCC‑CCA‑C‑TTTTTACCGC‑TTTTCCTGGGGGGG‑C‑GCGAA‑GTACAT‑CCCTAAAAAAAT‑ACCC‑ACCGC‑ACCCC‑ACGATTTTTG‑GTTTCCATGCAGGGCGAAAccc                                                                                                                           >  2:110230/1‑237 (MQ=255)                                                                                                                           cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTA‑AGGTT‑CGGGGATCAG‑CA‑AGCCCG‑CCC‑C‑TTTTTCCCCA‑TTTTTCCGGGCGGG‑A‑GTGGA‑GGACAT‑AGCTAAAGAAAA‑TCCC‑ACAGC‑C‑CCA‑AAAATTTTTTAGTTTCCGTGCAGGGCCAAACCG‑ccc                                                                                                                       >  1:1054337/1‑237 (MQ=255)                                                                                                                           cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGa                                                                                                                                                                                                                                                         <  1:349287/124‑1 (MQ=255)                                                                                                                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGTGGGGATCAC‑CA‑AGACCC‑CCA‑C‑TTTTACCCAC‑TTTTCCGGGTCGGCGA‑GGAAG‑GTACAC‑CCCAAAAAAAAT‑‑CAC‑ACGGC‑C‑CCACAAGATTTTTTA‑TTACCGTGCCGGGCCAAAcacccctcccgg                                                                                                                  >  2:344489/1‑228 (MQ=255)                                                                                                                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGAAGCCG‑CC‑AACCCC‑CCA‑C‑CTTTTCCCCG‑TTTTTCGGTGCGGC‑C‑GGGGA‑GGAAAA‑CGCCAAAAAAAT‑ATCC‑CCGGC‑C‑CCA‑AAGCTTTTTTTGTTGCCGTGCCGGGGAAAACCC‑CCTCTCC‑gg                                                                                                                >  2:407310/1‑239 (MQ=255)                                                                                                                                    cGACATAGAGGCCGCCAGCACGG‑GTGTACCAGCAAACAGAAAG‑GGAGAAATGAACAGGAAGATAAATGCCAGCGGTTCGGTAATTCCCCCAAGCAGGGC‑GG‑TTA‑AGGTGGCGGAAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGaaca                                                                                                           <  2:8122/239‑1 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTT‑AAGGTGCGGGGATCAG‑CA‑AGCCCC‑CCA‑C‑CTTTTCCCCA‑TTTTTCCGGGCGGC‑A‑GGGAA‑GTACAT‑CGCTAAAAAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGAATTTTTGGTTTCCGTGGCGGGGAAAAACG‑CCC‑TCC‑GGGgacaa                                                                                                          >  1:854681/1‑239 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTTCGGGGAACAG‑CA‑AGCCCC‑CCC‑C‑TTTTTCCCCA‑ATTTTCGGTGCGGC‑A‑GACAAGGAAAAC‑CCCTAAAAAAAT‑GCCC‑AC‑GC‑G‑CCA‑AAAGTTTTTTTGTTTCCGTGGCGGGGACAACCGCCCT‑TCG‑GGGaaca                                                                                                           >  2:985629/1‑239 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATCAG‑CA‑AGCCCG‑CCC‑CTTTTTACCCAT‑TTTTCCGG‑GCGGC‑A‑GGGAA‑GGACAT‑CGCTAAAGAAAT‑ACCC‑CCCGC‑A‑CCC‑AAGAATTTTTAGTTTACGGGCCAGGCAAAAACC‑CCC‑TCC‑GGGaaaaa                                                                                                          >  1:316317/1‑239 (MQ=255)                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTA‑AGGTTGCGGGAAAAAG‑CA‑AAACCG‑CCA‑C‑CTTTAACCGACTTTTCGGCGGGGGC‑A‑G‑GAA‑GGACAT‑CGCTAAAGAAAA‑GCCC‑CCGGG‑A‑CCC‑AAAAATTTTTAGATAACGGGCAGGGGAAAACAC‑CCC‑TTC‑GGGaaaaa                                                                                                          >  2:742988/1‑239 (MQ=255)                                                                                                                                          tGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGGATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCCGTGCGGC‑A‑GTGGA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCC‑AAAAATTTTTAGTTACCGTGCAGGGCGAAACCG‑CCT‑TCC‑GGGAACAACGAtt                                                                                                     >  1:647463/1‑239 (MQ=255)                                                                                                                                                 gccAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGAATCAA‑CA‑AACCCC‑CCC‑C‑CTTTACCCGA‑TTTTCCGGGGCGGC‑A‑GGGGA‑GGACAT‑CGCTTAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTTAGGTACCGTGCAGGGGAAAACCG‑CCC‑TCC‑GGGGACAAAGATTTCAGccg                                                                                              >  1:352671/1‑237 (MQ=255)                                                                                                                                                 gccAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTA‑AGGATGCGGGGAACAG‑CACAGCCGC‑CCC‑C‑TTTAACCGCG‑TTTTCGGGGGCGCC‑A‑GTGAG‑GACAAT‑CGCAAAAAAAAT‑‑CCC‑CCGGCAA‑CAA‑AAAATTTTTAGGTTACGTGGCAGGGCAAAACC‑‑CCT‑TCCGGGGaacaacatttttgcg                                                                                                >  2:461480/1‑229 (MQ=255)                                                                                                                                                  ccAGCACGGCGTGTAAC‑GCCAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGA‑GGCTTA‑AGGTCCGGGGAGTCAG‑‑A‑AGGCCG‑CCC‑C‑ATTTCACTGA‑GTTTCCGGGGGCGC‑C‑GTGAA‑GGCCCA‑CCCTCAAGACAT‑CCCCAACGCC‑C‑CCA‑AACTTCTTTGCG‑TACCATGTAGGGCAAAAACG‑CCTCTCC‑CGGaccaacatttttactgg                                                                                              >  2:836490/1‑228 (MQ=255)                                                                                                                                                    ggCAGGGCGAGTGCC‑GAAAACACCACA‑GGGAAAATGAACAGAAAGGTAACGTCCAGCGGTTGAATAACTCCCACCAGCGTGGA‑GT‑TTA‑AGGTGGGGGGAACCAG‑AA‑CGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTc                                                                                           <  2:854681/238‑1 (MQ=255)                                                                                                                                                    acggcGGCGTGTACC‑‑TAAACAACAACGGGAGAAATGAAGAAGAAGGAAAATTCGAGCGTTGCGATATTCCCCCCCAGCA‑GGC‑GG‑GTAGGGGGGGCGGAAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTc                                                                                           <  2:352671/235‑1 (MQ=255)                                                                                                                                                               tACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGAACAG‑CA‑AGCCCC‑CCACC‑TTTTACCC‑A‑TTTTTCGGGGCGGG‑AGGGGAA‑GAAAAC‑CCCAAAAAAAAG‑CCCC‑CCGCC‑C‑CCA‑AAGTTTTTTAGGTACCCGGGAGGGGCAAAACC‑‑CCC‑TCC‑GGGAACAAACACTTTCGCGGCCCGGC‑GCTCaaaa                                                                               >  2:38055/1‑238 (MQ=255)                                                                                                                                                               tACC‑GAAAACAGCAAG‑GATGAAATGAGGAGGAAGAAAAATTCCGGCGGTGCGTAAAGCCCCCCCAGCATGGC‑GT‑TGA‑GGGTCGCGGAAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCaaa                                                                                <  2:647463/238‑1 (MQ=255)                                                                                                                                                                    gaaaaCAGAAAC‑GGGAAAAGAAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGaa                                                                           <  1:461480/237‑1 (MQ=255)                                                                                                                                                                     cAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGATGCGGGGATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAG‑GGACAT‑CGCTCAAAGAAT‑GGCC‑ACGGC‑A‑CCC‑AAGAGTTTTGAGATATCGGGCAGGGGAAAAACG‑CCG‑TCC‑GGGGACAAAGATTTCACTGGGTCGGC‑GCCCACACCGaa                                                                           >  2:141589/1‑238 (MQ=255)                                                                                                                                                                      aaaCAGCAGC‑GGTAAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATAAA‑CA‑ACCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAAct                                                                         <  1:613355/239‑1 (MQ=255)                                                                                                                                                                        acaaaaaC‑GGTAAAATGAACGGGAAGGTAAATCCCAGCGGTTCGTTAATCCCCACCAGCAGGGC‑GT‑TTA‑AGTTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAActct                                                                       <  1:458316/234‑1 (MQ=255)                                                                                                                                                                               c‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGGTGCGGGGATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTTAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGAt                                                                >  1:625382/1‑239 (MQ=255)                                                                                                                                                                                    gAAATAACCAGCAAGGTACATCCCAGTGGTGGTATACTCCCCACGAGCAGGGG‑GG‑TTG‑AGGTTGGGGGAAGCAC‑CA‑GTCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGAtgct                                                             <  2:922972/238‑1 (MQ=255)                                                                                                                                                                                     aaaTGAACAGGAAGGTAAATTCCAGCGGTACGGTAATTCCCACCAGCATGGC‑GG‑TCA‑AGGTTGCTGGGATCAG‑CC‑AGCCCC‑CCC‑C‑TTTTTTCCGA‑TTTTCCGGGGCTGC‑C‑GTGCA‑GTGCAC‑CCCTAACGCAAA‑CCCC‑CCAGC‑A‑CCC‑AAGTATTTTGAATTTATGGAAAGGGCCAAACCA‑CCC‑TCC‑CCGGCAAACGATTTCAAGCGCTCGGCGGCTCAAACCGAAccctcctcgaaccgc                                                            >  2:829189/1‑228 (MQ=255)                                                                                                                                                                                        tGAACAGGAAGATAAATTCAAGCGGTTCGATAATCCCCACCAGCATGGC‑GT‑TTA‑AGGTTGCGGAAATCAG‑CA‑GGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  1:141589/239‑1 (MQ=255)                                                                                                                                                                                        tGAACAGAAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCAGGGC‑GG‑TTA‑AGGTGGGGGGAACCAG‑AA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  1:110230/239‑1 (MQ=255)                                                                                                                                                                                             agaaaGGTAAATCCCGGCGTTTGGGTAATCCCCCCCAGCATGGG‑GG‑TTA‑GGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                   <  1:481540/236‑1 (MQ=255)                                                                                                                                                                                             aGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGC‑GG‑TTA‑GGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                   <  1:1017159/239‑1 (MQ=255)                                                                                                                                                                                                    tAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCCA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAA‑GCCC‑CCGGC‑A‑CCA‑AAGATTTTTTAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCCGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGa                                            >  1:440178/1‑239 (MQ=255)                                                                                                                                                                                                                             ccccaGCAGGCC‑GG‑TTA‑AGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  1:601558/237‑1 (MQ=255)                                                                                                                                                                                                                             cacaaacATGCC‑GG‑TTA‑GGTTAGCGGGAACCAA‑CA‑CGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  1:829189/233‑1 (MQ=255)                                                                                                                                                                                                                               ccagatgggC‑GG‑TTA‑AGGTGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                 <  1:877122/232‑1 (MQ=255)                                                                                                                                                                                                                               ccaGCAGGGC‑GG‑TTA‑AGGGGGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                 <  1:407310/239‑1 (MQ=255)                                                                                                                                                                                                                                     tGGC‑GG‑TAA‑ATCTTGCGGAAACCGC‑CA‑ATCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGAt            <  1:239702/238‑1 (MQ=255)                                                                                                                                                                                                                                           g‑TTA‑AGGTGGCGGAAAACAG‑CA‑GGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACt       <  1:742988/238‑1 (MQ=255)                                                                                                                                                                                                                                              tA‑AGGTTGCGGGAATAAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc    <  1:38055/239‑1 (MQ=255)                                                                                                                                                                                                                                                 aGGTGGCGAGAATCAC‑CA‑GCCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCg  <  1:257015/239‑1 (MQ=255)                                                                                                                                                                                                                                                    ttGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCCCAAACTGAACTCTTGCAGCTGCTGCGCCCCGTACATACGAATGCCGCCTTCAACAGCTGCCGTACCAAAg              >  2:767903/1‑224 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                               AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC‑GCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTA‑AGGTTGCGGGAATCAG‑CA‑AGCCCG‑CCA‑C‑TTTTACCCGA‑TTTTCCGGTGCGGC‑A‑GTGAA‑GTACAT‑CGCTAAAGAAAT‑GCCC‑ACGGC‑A‑CCA‑AAGATTTTTGAGTTACCGTGCAGGGCAAAACCG‑CCT‑TCC‑GGGAACAACGATTTCAGCGGCTCGGC‑GCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCG  >  NC_000913/3862423‑3862896

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