Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,674:1	+A	6.9%	pseudogene (930/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,674	1	.	A	6.9%	190.4 / 3.9	58	pseudogene (930/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 . (24/30);  new base A (2/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

AGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTG  >  NC_000913/3862444‑3862907                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  aGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGGTAA‑GG‑TTGCTGGGATCAGC‑AAACCC‑g                                                                                                                                                                                                                                           >  1:266808/1‑238 (MQ=255)   catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAAAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGATCGGTAATTCCCACCAGCAT‑GGC‑GGGTAA‑GGCATGCGGGATGCAGC‑AACCTC‑G‑c                                                                                                                                                                                                                                         >  2:141770/1‑238 (MQ=255)           ccGGGTGGTTGAAGAACTGCGGGATCCTTTGTTGCGGAAAACCCTAG‑CAAACCGCCCCCCCCCCATTGTCCCCCCCCACCCC‑AAAGAAATACATTCGGGTCAACATGGGGCCCGCCAGCAGGGCGTGTACCGAAAACAGCAAC‑GGTAAAATGAACAGGAAGAAAAATCCC‑A‑GC‑GTTTCGGTATTCCCCCCCAGCAG‑GGC‑GTTTAA‑GG‑TGGCGGAAATCAGC‑AAGCCC‑G‑CCA‑CTTTTAcc                                                                                                                                                                                                                              <  2:254878/239‑1 (MQ=255)             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGGTAA‑AG‑TTGCGGGGAACAGC‑AAACCC‑G‑CCA‑CTTTTTcccc                                                                                                                                                                                                                            >  1:729306/1‑238 (MQ=255)              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGCGGGTTAA‑GG‑TGTGGGGAAACAGC‑AAGaccccacccttta                                                                                                                                                                                                                                   >  2:440129/1‑225 (MQ=255)              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGGTCGGTAATTCCCACCAGCAT‑GGCGGGTAAA‑GG‑TTTGGGGAA‑CAGC‑GAGCCCGG‑CCA‑Attttcacca                                                                                                                                                                                                                            >  2:928899/1‑234 (MQ=255)                     gCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGGGATACGC‑AAAACC‑C‑CCA‑CCTTTATCCGCATTTTcg                                                                                                                                                                                                                    >  1:234113/1‑238 (MQ=255)                         gacctCGGGATCCCGTTATGCCGTAAACCCCAAGCAAACCGACCCCCCGCCA‑TGCTCCCCCCCACACCAAAAGAAATCCTTTACGGTGGACATGGAGCCCGCCAGCGGGGCGTGTCCCCCAACAGGCACC‑GGTAAAAAGAACAGGAAGAAAAATTCG‑A‑GG‑GGTTGGGTAATCCCCCCCAGCAT‑GGG‑GTTTAG‑GG‑TGGCGGAAACCAGC‑AGGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTg                                                                                                                                                                                                                <  2:1040394/235‑1 (MQ=255)                         aaCATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGCGGGTTAA‑GG‑TTGTGGG‑ATTAGC‑AAACCC‑C‑CCACCTTTTCCCCC‑CTTTTCCgtgt                                                                                                                                                                                                                >  1:1040394/1‑236 (MQ=255)                                 gaccctgTTGTGCGAAAAAGCCAAGACAATCCCCCCCCCGCCCA‑TTTCCCCCACCACACA‑AAAGAAATAATTTCCGGACGACATAGAGCCCGCCAGCACGGCGTCTCCCGCAAACAGCACG‑GGGAAAATGAACGGAAAGGAAAATCCC‑G‑GC‑GTTTCGGTATTCCCCCCCAGCAG‑GGG‑GTTTAG‑GG‑TTGCGGGAACCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCTGC‑AGTg                                                                                                                                                                                                       <  2:413728/232‑1 (MQ=255)                                      aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGGTAAGGG‑TTTGGGGAATAGGA‑AACCCC‑C‑CCC‑CTTTTCCCCC‑GTTTTCC‑GTGGGGC‑AGGGacgga                                                                                                                                                                                                  >  1:575811/1‑234 (MQ=255)                                       gTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGCGGGTGAG‑GG‑TTTGTGGGAGATGC‑AAGCAC‑C‑CCA‑CCTTTACCCA‑AGTTTCCGGggcggcggggaggat                                                                                                                                                                                                   >  2:1065788/1‑229 (MQ=255)                                           tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGCAATTCCCACCAGCAT‑GGC‑GGGTAA‑GG‑TGTCGGGGATCAGC‑AAGATC‑GACCA‑CTTTTTCCCA‑CTTTTCTGGGTGGGC‑AGTggggaaatca                                                                                                                                                                                              >  2:496646/1‑229 (MQ=255)                                             cGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGGTTA‑AG‑TTGCGGGGATCAGC‑AAAGCC‑G‑CCACCTTTTACCCG‑ATTTTCC‑GTGCGGC‑CGTGAAgga                                                                                                                                                                                                  >  1:769573/1‑229 (MQ=255)                                                  aaaCCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGGTAA‑AG‑GTGCGGGGATCAGC‑AAACCC‑G‑CCA‑CCTTTTACCG‑CTTTTCCGGTGCGGC‑AGGGGAGTACAATGCTCAAG‑aa                                                                                                                                                                                    >  1:1047408/1‑239 (MQ=255)                                                              tCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTCAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTTA‑AG‑GTGCGGGGATCAGC‑AAACCC‑G‑CCA‑CTTTTACCCG‑CATTTCCGGGGGCGC‑AGTGAAGGACATAGCTACAG‑CAATACCCCCGGc                                                                                                                                                                         >  1:449493/1‑238 (MQ=255)                                                                       cGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑GTGCGGGGATAAGC‑AAACCC‑G‑CCA‑CCTTTAACCC‑ATTTTTCGGTGCGGC‑AGTGAAGTAAATCGCTAAAG‑AAAAGcc                                                                                                                                                                               >  1:648891/1‑223 (MQ=255)                                                                       cGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGGAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGCGGGTTAA‑GG‑GT‑TGGGGATCAGCAAACCCC‑C‑CCC‑CTTTTCCCCG‑GTTTTCCGGT‑CGGCGAGGGCAGTAAATCGCCAAAG‑AAATGCCCCCGGCAcccagaat                                                                                                                                                                >  2:540223/1‑233 (MQ=255)                                                                                                     aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGGTAA‑GG‑GTGCGGGGATCAGC‑CAAACC‑G‑CCC‑CCTTTACCCG‑CTTTTTCGGTTGGGC‑AGGGAAAGAAACAGCTAAAA‑AAAAATCCCCACGACCACAAAATTTTGAGTTATCGCGCAGGGCGAAAcag                                                                                                                                    >  2:1097253/1‑236 (MQ=255)                                                                                                      tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑GTTGGGGGATCAGC‑CAGCCC‑G‑CCA‑CTTTTACCCC‑CTTTTTCGGTGCGGC‑AGTGAAGTACATCCCCAAAA‑AAAAGGCCACCGGCCCaaaaa                                                                                                                                                                 >  1:1071099/1‑206 (MQ=255)                                                                                                      tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGGTAA‑GG‑GTTCGGGGATCAGCAAAACCC‑C‑CCC‑CTTTTCCCCC‑TTTTTCTGGTGC‑GC‑CGTGAAAGACAACGCCAAAGAAAAATCCCCCGGC‑CCCAAGAATATTGTGTTACCGCGCAcggtacacaacc                                                                                                                                   >  2:963942/1‑229 (MQ=255)                                                                                                              gtaccaCATGGAGACCGGCAGCACGGCTTTCACAGAAAACGACAAG‑GGAAAAAGGAAGAGAAATGTAACTTCC‑A‑GC‑GTTTCAGTACCCCCCCCCGGCAG‑GGC‑GTTTAA‑GG‑GGGCGGAAAACACA‑AACCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGa                                                                                                                          <  2:480623/234‑1 (MQ=255)                                                                                                               ggcGACATTGGAGCCCCCAGGCCGGGGTGTACCCAAAAACGGAACGGGAAAAAGGAAAGGGAAGAAAATTTCC‑G‑GG‑GTTTCGAAAATCCCCCCCAGC‑G‑GGG‑GGTTGA‑GG‑GTTGGGGAACCCCCAAAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGa                                                                                                                          <  2:296738/237‑1 (MQ=255)                                                                                                               gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑GTGCGGGAAACAGC‑AAAACC‑C‑CCC‑CTTTTACCCG‑ATTTTTCCGGGCGGC‑CGTGGAgga                                                                                                                                                                                                  >  1:1136407/1‑165 (MQ=255)                                                                                                                                      ggggtgtCCCGAAACAAGCAAC‑GGTAAAAAGAAGAGGAAGGAAAATCCC‑A‑GC‑GGTTCGGTAATCCCCACCAGCAG‑GGC‑GGTTAA‑GG‑TGGGGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaacaa                                                                                                                      <  1:124145/216‑1 (MQ=255)                                                                                                                                           cgaccaaaaaCAGAAAG‑GGAAAACTGAAAAGGAA‑GTAAACTCCGA‑GC‑GTTTCGAAAATCCCACCCAGCGG‑GGG‑GTTTAG‑GG‑TGGCGGAAACAAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑aaa                                                                                           <  2:726184/232‑1 (MQ=255)                                                                                                                                                          aagggTAAAAGGAACGGGAAGATAATTCCC‑G‑GG‑TGTTCGGTATTCCCCACCAGCTG‑GGG‑GGTTAG‑GG‑TGGGGGAAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAg                                                                             <  1:220937/236‑1 (MQ=255)                                                                                                                                                               tGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGGGATCAGC‑CAGCCC‑C‑CCC‑CTTTTACCCG‑CTTTTTCGGGGCGGC‑AGTGGAGGACATCGCTAAAA‑AAATGCCCACGGCACCCAAGATTTTTTAGTTACCGTGCAGGGCaaaa                                                                                                                                       >  1:1114816/1‑178 (MQ=255)                                                                                                                                                                    tGACCAGAAAGGTAAATCCC‑G‑GG‑GGTTCGGTACTCCCCCCCAGATT‑GGG‑GGTTAG‑GG‑TGGGGGGAATCAGA‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:467996/239‑1 (MQ=255)                                                                                                                                                                    tGAACAGGAAGGTAAATCCA‑G‑GC‑GTTCCGGTAATCCCCACCAGCAG‑GGC‑GGTTAA‑GG‑TGGGGGAAAACAGA‑AGGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:1241820/239‑1 (MQ=255)                                                                                                                                                                      aagaGGAAGATAATTCCC‑G‑GG‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTAAG‑GT‑TTGCGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  1:824830/234‑1 (MQ=255)                                                                                                                                                                           aatggTAAATCCC‑G‑GG‑GGTTCGGTAATCCCCACCAGATG‑GGC‑GTTTAA‑GG‑TGGGGGAAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAg                                                                <  1:43876/232‑1 (MQ=255)                                                                                                                                                                            gaGGGAAAATCC‑A‑GGGGGTTCGGAAATCCCCACCAGCAG‑GGC‑GTTAAA‑GG‑TTGCGGAAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTgcg                                                                     <  1:157002/229‑1 (MQ=255)                                                                                                                                                                             atgataaTTCC‑A‑GG‑GTTTGGGTAACCCCCCCCCGGCA‑GGC‑GGTTAA‑GG‑GTGGGGGAAACAGCAAACCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGcc                                                                   <  2:205807/226‑1 (MQ=255)                                                                                                                                                                                 aaaTTCC‑CGGC‑GGTTGGGTAATCCCCCCCAGCAG‑GGG‑GGTTAG‑GG‑TGGGGGAAACCAGA‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACAt                                                                                                                                                                                               <  1:753474/106‑1 (MQ=255)                                                                                                                                                                                   atccc‑A‑GC‑GTTTCGGTAATCCCCACCAGCAG‑GGG‑GGTTAA‑GG‑TTGGGGAAATCAAG‑AACCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                    <  1:1027523/236‑1 (MQ=255)                                                                                                                                                                                     tCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑GTGCGGGGATCAGC‑AAACCC‑C‑CCACCTTTTACCC‑‑ATTTTTCGGTGCGGC‑AGTGAAGGACAACGCTCAAA‑AAAAGCCCCCCGCCCCAAAGAATTTTTAGTTAACGGGCAGGGGCAAAACGCCTTCCGGGGACAAACATTTCTACGGCTCGGGGGCCC‑AAACTGAACACTTCCCGATGTTGGCCCCCGAACAt                                                           >  2:157002/1‑231 (MQ=255)                                                                                                                                                                                     atcccggC‑GTTTCGGTAATCCCCACCAGCAG‑GGC‑GGTTAG‑GG‑TGGCGGAAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                    <  1:440129/233‑1 (MQ=255)                                                                                                                                                                                        ccggC‑GGTTGGATAATCCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGAAATCACC‑AGACCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccg                                                 <  1:496646/236‑1 (MQ=255)                                                                                                                                                                                                ttggatAATCCCCCCCAGCTG‑GGC‑GTTAAG‑GG‑TGGGGGGAACCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGAtgct                                                                        <  1:456072/204‑1 (MQ=255)                                                                                                                                                                                                    gTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGGGATCAGC‑AAACCC‑C‑CCC‑CCTTTTCCCG‑CTTTTTCCGTGCGGG‑AGTGAAGTACATCGCTAAAG‑AAAAGCCCACGGCACCACAGATTTTTTAGTTACCCTGCAGGGCAAAACCGCCTTCCCGGGAAAACGATTTTAGCGGGTC‑CGCGGTCTAAACGAAAACCTTGCGAAGACTGGGCCCGGTACATCTAATTGCGCCCtttcaaa                                        >  2:348385/1‑233 (MQ=255)                                                                                                                                                                                                        atCCCACCAAGATGGGG‑GGTTAG‑GG‑TTGGGGAAACGAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                                   <  1:963942/238‑1 (MQ=255)                                                                                                                                                                                                            caccgGGAT‑GGC‑TGTTAA‑GG‑TTGCGAGACTCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                              <  1:88713/239‑1 (MQ=255)                                                                                                                                                                                                            caccaGCAT‑GGC‑GGTTAG‑GG‑TTGCGGAAATCGGC‑AACCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                              <  1:892883/239‑1 (MQ=255)                                                                                                                                                                                                                     gggC‑TGTTAG‑GG‑TTGGGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTGCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGAt                       <  1:590362/237‑1 (MQ=255)                                                                                                                                                                                                                        c‑GGTAAG‑GG‑TTGGGGAAAACAGC‑AACCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                    <  1:348385/238‑1 (MQ=255)                                                                                                                                                                                                                         cgtttAG‑GG‑TGGGGGAAACCAGC‑ACGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                    <  1:519811/235‑1 (MQ=255)                                                                                                                                                                                                                           gTTAA‑GG‑TTGGGGAAATCAGA‑AGCCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACt                  <  1:408881/238‑1 (MQ=255)                                                                                                                                                                                                                              taggg‑TTGGGGAAACCAAC‑AGGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc               <  1:1193931/236‑1 (MQ=255)                                                                                                                                                                                                                                    tggcggaaacaaGC‑ACGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt            <  1:368530/226‑1 (MQ=255)                                                                                                                                                                                                                                      gCGGGAATGAGC‑AAGCCC‑G‑CCA‑CTTTGACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGa         <  1:928899/238‑1 (MQ=255)                                                                                                                                                                                                                                      cgggaaaccgac‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGa         <  1:540223/227‑1 (MQ=255)                                                                                                                                                                                                                                       cGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTTCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg       >  1:1243109/1‑239 (MQ=255)                                                                                                                                                                                                                                       cGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                            >  2:908576/1‑218 (MQ=255)                                                                                                                                                                                                                                       cGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAACTTCCGGACCAAAAGTAAACTGTCCGTAGATg       >  1:392126/1‑239 (MQ=255)                                                                                                                                                                                                                                       cGGAAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                            <  1:908576/218‑1 (MQ=255)                                                                                                                                                                                                                                          gAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa     <  1:1097253/238‑1 (MQ=255)                                                                                                                                                                                                                                          aaaTAAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa     <  1:1195148/237‑1 (MQ=255)                                                                                                                                                                                                                                            aTCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtg  >  2:243785/1‑239 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  AGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTG  >  NC_000913/3862444‑3862907

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