Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,672	G→A	19.5%	pseudogene (931/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,672	0	G	A	19.5%	62.8 / 2.8	48	pseudogene (931/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 G (21/17);  new base A (3/6);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTG  >  NC_000913/3862450‑3862906                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGAT‑GCGGG‑GATTA‑GCAAAG‑CCCC‑CCAt‑ttt                                                                                                                                                                                                                                     >  2:202776/1‑239 (MQ=255)      gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCCAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAGGGT‑GGGGG‑AATCACGC‑AA‑‑CCCG‑CGCC‑TTTTA‑ccc                                                                                                                                                                                                                               >  2:954372/1‑238 (MQ=255)      gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTT‑GGGGG‑GAACA‑GC‑CAG‑CCCG‑CCCC‑TTTTA‑ccc                                                                                                                                                                                                                               >  2:1138245/1‑238 (MQ=255)       cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑GATTA‑GA‑GAA‑CCCC‑CCCC‑TTTTT‑cccca                                                                                                                                                                                                                             >  2:436527/1‑237 (MQ=255)           gTTGCTGAACATCGGGATACAGTTTTGCGGTAAAACCTAGTCAATCAGACCCCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCCACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGGAATTCCCACCAGCATGGC‑GGGTCAGGGT‑GCGGG‑GATTA‑GC‑AAgaaactcacatttatacgcattt                                                                                                                                                                                                                             >  2:614148/1‑217 (MQ=255)                          gATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GAACA‑GC‑AAC‑ACCGCCCCC‑TTTTA‑CCC‑‑ATTTTTCGGTGC‑gggaggga                                                                                                                                                                                                        >  2:797896/1‑239 (MQ=255)                            tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GATCA‑GC‑AAA‑CCCG‑CCCC‑TTTTT‑CCCC‑CTTTTCCCGGGC‑GGCAGGGaaa                                                                                                                                                                                                      >  2:750574/1‑238 (MQ=255)                              cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATAAGGTT‑GCGGG‑GAACC‑GC‑AAA‑CCCC‑CCCC‑CTTTA‑ACCG‑CGTTTTCTGTGC‑ggtcgtggagga                                                                                                                                                                                                    >  1:140429/1‑229 (MQ=255)                               attttcGCGGAAAAACATAGCAAAACACCCCCCCGCCCGTGTCCCCCACCCCACAAAAGAGAACCATTACGGCCGACATTGAGGCCGCCAGCAC‑GGCGTGTCCCGCAAACAGCAGAGGTAAAATGAACAGGAAGGTAAA‑TTCAAGCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAGGGTT‑GCGGA‑AATAA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTAc                                                                                                                                                                                                   <  1:291425/237‑1 (MQ=255)                                 ttttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GATCA‑GC‑AAA‑CCCG‑CCCC‑TTTTT‑CCCC‑ATTTTTCGGTGC‑GGCAGGGAAGTACa                                                                                                                                                                                                  >  2:998706/1‑238 (MQ=255)                                           aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAGGGT‑GCGGG‑GAACA‑GC‑AAG‑CCCG‑CCCC‑TTTTT‑ACCC‑ATTTTTCGGGGC‑GGCAggggaggaacttggctaagaa                                                                                                                                                                                       >  2:236363/1‑219 (MQ=255)                                                                cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTCAAGGG‑GCGGG‑GAGCA‑GC‑ACA‑CCCC‑CCAC‑TTTTA‑CCCC‑TATTTTCGCGGC‑GGCACGGACGGACAT‑CGCTAAAGAACTACCCA‑CGGGCCCAAAGagt                                                                                                                                                                >  1:498663/1‑237 (MQ=255)                                                                 cccagttttGCCCACCCCCCCAAAGAAAAACTTTAGGGGCGAAGTTGA‑GCCCCCCGCGCGGGCGTGTACCGCAAACAACAAGGGTAAAAGGAACAGGAGTGAAAA‑TTCCGGTGGTTCGGAATTCCCCACGAGTATGGC‑GGTTAAGGTT‑CGGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAt‑tt                                                                                                                                                              <  1:599458/232‑1 (MQ=255)                                                                    caggTTGCCCACCCCACAAAAAAAATACATTACGGCCAACATGAAGCCCCCCAGCAC‑GGCGTGTCCCGCAAACAGAACCGGGGAAAGGAACAGAAAGTTAAA‑TTCCGGCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAAGTTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTg                                                                                                                                                           <  2:449534/236‑1 (MQ=255)                                                                              accacacaaAAGAGATACATAACGGTCGACATTGAGGCCCCCAGCAC‑GGGGTTTCCCGCAAACAGCAACGGTAAAAGGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTg                                                                                                                                                 <  2:335246/239‑1 (MQ=255)                                                                                cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAAGGT‑GCGGG‑GAACA‑GC‑AAGACCCC‑CCCT‑TTTTA‑CCCG‑ATTTT‑CGGTGG‑GGCAGTGAAGGACAA‑CGCTAAAAAAAAGCCCA‑CCGGCCCCAAGAA‑TTTTTAGtttaccgggc                                                                                                                                               >  2:1182898/1‑231 (MQ=255)                                                                                            aTACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AATAA‑GC‑AAA‑CCCC‑CCACTTTTTA‑CCCG‑TTTTTCGGGGCG‑GGCGGTG‑AGAACAA‑ACCTAAAAAAAAGCCCC‑CAGGCCCAAAGAT‑TTTTTAGTGTGCGGGCCAGGC‑AAAACCg                                                                                                                                   >  1:1110160/1‑238 (MQ=255)                                                                                                aTTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑Gc                                                                                                                                                                                                                                                       <  2:926555/129‑1 (MQ=255)                                                                                                      gTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTG‑CGGGG‑AACAA‑CC‑AAG‑CCC‑‑CCCC‑CTTTA‑CCCC‑CTTTTCCGGTGG‑GGCAAGGAAGAACAC‑CCCTAAAAAAAAACCCC‑CGGCACCAAAGATGTTTTGAGTTACC‑TGGCGGGC‑AAAACCCCCCTCCcggca                                                                                                                        >  1:644905/1‑237 (MQ=255)                                                                                                      gTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTTAGGGT‑TCGGG‑AATCA‑GC‑CAG‑CCCC‑CCCC‑CTTTA‑ACCGAATTTTCCGG‑GC‑GGGAGGGAAGGACATCCCCCAAAAAAAACCCCA‑CGCCCCCAAAGAT‑TTTTGAGTCA‑CGTGCCCGGC‑AAAAcaccccctcgggaa                                                                                                                        >  1:1115232/1‑226 (MQ=255)                                                                                                       tCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AATCG‑AC‑AAG‑CCCC‑CCCC‑TTTTT‑CCCG‑TTTTTCTGGGGC‑GGCCGAGAAAGACAA‑AACTTAAAAAATTGCCA‑CAGCACCCAAGAT‑TTTTAAGGTACCGTTCTCTGC‑ATAAACGCCCTCCggggaa                                                                                                                       >  1:954132/1‑236 (MQ=255)                                                                                                              tGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTT‑CCGGG‑AA‑CA‑AC‑AAGCCCCC‑CCCC‑TTTTC‑CCCG‑TTTTTCGGGGGC‑‑GCAGGGGAGGACAT‑ACCTCAACAAATGCCCACCGGACCCAAAAAT‑TTTTTGTTTCCCGCGCGGGGA‑ACAACCCCCTCCCGGGaacacgcat                                                                                                                 >  1:840190/1‑234 (MQ=255)                                                                                                                          cAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGG‑T‑GCGGG‑AATCA‑GC‑AAC‑ACCG‑CCCC‑CTTTT‑CCCC‑ATTTTCTCGTGCGGGCGGGGAAG‑AAAA‑CGCCTAAAAAATGCCCC‑CCGCACCAAAAAA‑TTTTTAGTTACCCCGCCGGGG‑AAAAACACCCTCCGTGgacaacaactttcacggccccg                                                                                                    >  1:870343/1‑223 (MQ=255)                                                                                                                           ccgcgggtgtACCGAAAACAACAACTGTAGAATAAACAGGAATGTAAA‑TCCCAGC‑TTTCCGTAACTCCCCCCAAGGAGGG‑GGCTTAGGTT‑GGGGGCAACAA‑GC‑AAC‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCgg                                                                                                   <  1:1107366/233‑1 (MQ=255)                                                                                                                           cc‑GGCGTGTACCGAAAACAGCAAGGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTGGGTAATCCCCACCAGCATGGC‑GGTTAAGGTT‑GGGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTc                                                                                                     <  2:140429/235‑1 (MQ=255)                                                                                                                            agggggtgacccgaaaacggataCGGTGCATTGACCAGGAAGGTAAC‑TCCCAGCGGTCCGTCATTCCCCCCCGGCTGGGC‑GATAAAGGTT‑GGGGG‑AACCG‑AC‑ACC‑GCCG‑CCAC‑GTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTc                                                                                                     <  1:954372/214‑1 (MQ=255)                                                                                                                              ggCGTGTACCGAAAACAGCAACGGTGAAATAAACAGGAAGGTAAA‑TCCCAGCGGTTCGTTAATCCCCACCAGCAGGGC‑GGTAAGGTTT‑GGGGA‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTg                                                                                                                                                                                                                  <  2:761077/131‑1 (MQ=255)                                                                                                                                        gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGT‑GCGGG‑GATCA‑GC‑AAG‑CCCG‑CCAC‑CTTTACCCCC‑ATTTTCC‑GCGC‑GGCCGTGAAAGAAAT‑ACCCAAAAAAAAGCCCC‑ACGCACCCAAAAA‑TTTTTAGTTACAGGTCGGGGCAAAAACCCCCTTCCGGG‑ACAACGAATTTCacgggccccggcgccaacctg                                                                                        >  1:1167679/1‑218 (MQ=255)                                                                                                                                          aaaCAGCAACGGGGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAGGTT‑TCGGG‑AAATC‑GC‑AAA‑CCCC‑CCAC‑TTTTA‑CCCG‑AGTTTTCCGTGC‑GGCCGGGAAGGACAT‑CGCTAAAAAAAAAACCA‑CGCTCCCCAAAAT‑TTTTTAGTTATCGTGCAGGGC‑AAAAACGCCCTCCGGGAACAAAG‑ATTTTAggcgcgccgggtccaatccacac                                                                                     >  1:851352/1‑216 (MQ=255)                                                                                                                                                       gAAATGAAGAGGAAGGTAAA‑TCCCGGCGTTTCGTAAATTCCCACCAGTATGCG‑GGTTAAGGTT‑GGGGA‑AATCA‑GA‑AAA‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                        <  2:851352/238‑1 (MQ=255)                                                                                                                                                       aaaagaaaaaggaaggaAAA‑TTCCGGGGGGTGTGTAATCCCCCCCCGCGGGGG‑GGTTAGGGGG‑GCGGA‑AACCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                        <  1:797896/229‑1 (MQ=255)                                                                                                                                                        aaaGGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAC‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                             <  2:367797/232‑1 (MQ=255)                                                                                                                                                             aaCAGGAAGGAAAA‑TTCCGGGGTTTCGGTAATTCCCCCAAGCATGGC‑GTTAAAGGTG‑GCGGA‑ACTCA‑GA‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:1182898/237‑1 (MQ=255)                                                                                                                                                               cgaagaagTAAATTTCCGGCGGTTCGAAAACCCCACCCGGCAGGGG‑GGTTAAGTGT‑GGGGT‑AATCA‑AC‑GAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                                <  1:113379/232‑1 (MQ=255)                                                                                                                                                                       tAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAAGGT‑GCGGG‑GATCA‑AC‑AAA‑CCCC‑CCCC‑TTTTA‑CCCG‑ATTTTTCGGTGC‑GGGAGTGAACTACAT‑CGCCCAAGAAAATACCA‑CCGGCCCCAAGAA‑TTTTGAGTTACCGTGCAGGGC‑AAAAAAGCATTCCCGGGACAAAC‑ATTTCAGACGCTCCGTGGTCAAAATGAAATCTTTCACAtggcggcccccggtcacacag                                                       >  1:1024910/1‑220 (MQ=255)                                                                                                                                                                           ttccgtgtggttgtaatatccccACCAGCATGGG‑GTATAAGGTG‑GCGGG‑AAAGA‑AC‑AAC‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                    <  1:781462/219‑1 (MQ=255)                                                                                                                                                                             tcgatgggTACGGTAATCCCCCCCAGCAGGGC‑GGTTAAGGTG‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgcc                                                  <  2:498663/233‑1 (MQ=255)                                                                                                                                                                             tcccggggtcgggaaaTTCCCACCAGCATGGC‑GTTTAGGGGT‑GCGGG‑AATCA‑AC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgcc                                                  <  2:1010759/225‑1 (MQ=255)                                                                                                                                                                                            acccccACCGGCAGGGC‑GTTAAAGGTG‑GCGGA‑AACCA‑AC‑AGG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                                   <  2:209973/236‑1 (MQ=255)                                                                                                                                                                                                 caccaGCA‑GGG‑GGTTTAGGTTGGCGGG‑AACAA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                              <  2:1115232/239‑1 (MQ=255)                                                                                                                                                                                                       cATGGC‑GGTTAAGGTG‑GCGGG‑AATCA‑GC‑ACG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                        <  2:981100/191‑1 (MQ=255)                                                                                                                                                                                                             ggttaaaGGTG‑GCGGG‑AAGCA‑GC‑AGC‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                    <  2:398086/233‑1 (MQ=255)                                                                                                                                                                                                             cgtttAGGGTG‑GGGGA‑AATCA‑GA‑ACG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                    <  2:1124370/235‑1 (MQ=255)                                                                                                                                                                                                                          cGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCTGCGCTCAAACTGAACTCTTTCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg       >  2:145511/1‑239 (MQ=255)                                                                                                                                                                                                                          cGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg       >  2:334974/1‑239 (MQ=255)                                                                                                                                                                                                                             g‑AACCA‑AC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa     <  2:1151556/238‑1 (MQ=255)                                                                                                                                                                                                                                aTCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtg  >  2:797070/1‑239 (MQ=255)                                                                                                                                                                                                                                aTCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtg  >  1:818932/1‑239 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       TATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTG  >  NC_000913/3862450‑3862906

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