Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,660	T→G	34.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	34.3%	31.4 / 12.7	50	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 (10/21);  new base G (14/3);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 2.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862424‑3862888                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑g                                                                                                                                                                                                                                                  <  2:262285/239‑1 (MQ=255) cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAC‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAGGTT‑Gc                                                                                                                                                                                                                                                 >  2:885447/1‑239 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTTCCCCCAGCATGGC‑GGTTAAGGGT‑GCGGG‑GAtgat                                                                                                                                                                                                                                       >  2:17323/1‑236 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑GAtac                                                                                                                                                                                                                                        >  2:1025217/1‑236 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGA‑ACATTACGGTCGACATTGAGGCCGCCCGCAC‑GGCGTGTACCGCAAACAGCAACCGTGAAATGAACAGGAAGGTAAA‑TTCCACCGGTTCGGTCATTCCCCCCAGCATTGC‑GGTTAAGGTT‑Gctgtgattcg                                                                                                                                                                                                                                        >  1:1112886/1‑229 (MQ=255)                           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)                                                                                                                                                                   cAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCggg                                                                                                                                                                                                                                              <  2:650812/82‑1 (MQ=255)                                                                                                                                                                   cAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGT‑GCggg                                                                                                                                                                                                                                              >  1:650812/1‑82 (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)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTA‑CCCG‑ATTTTCCGGTGC‑GGCAGTGAAGTACAT‑CGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGC‑AAAACCGCCTTCCGGGAACAACG‑ATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862424‑3862888

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