Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,662	T→G	25.6%	pseudogene (942/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,662	0	T	G	25.6%	59.0 / 2.3	52	pseudogene (942/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 (21/16);  new base G (2/13);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 5.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTC  >  NC_000913/3862425‑3862894                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑AG‑ttt                                                                                                                                                                                                                                                  >  1:726184/1‑238 (MQ=255)  ccAATGGCGATTTGGGTCAGCATCATAGCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGc                                                                                                                                                                                                                                                 <  1:774191/239‑1 (MQ=255)                    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)                                                                                                                                                                                          aCAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTg                                                                                                                                                                                                                                                  >  2:910479/1‑57 (MQ=255)                                                                                                                                                                                          aCAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTg                                                                                                                                                                                                                                                  <  1:910479/57‑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)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑TGTTGCCCACCACACC‑AAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAC‑GGTGAAATGAACAGGAAGGTAAATTCC‑A‑GC‑GGTTCGGTAATTCCCACCAGCAT‑GGC‑GGTTAA‑GG‑TTGCGGGAATCAGC‑AAGCCC‑G‑CCA‑CTTTTACCCG‑ATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAG‑AAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTC‑GGCGCTC‑AAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTC  >  NC_000913/3862425‑3862894

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