Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,662	T→G	21.2%	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	21.2%	84.1 / 3.0	58	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 (30/14);  new base G (4/10);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862425‑3862889                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTg                                                                                                                                                                                                                                               <  1:1213083/239‑1 (MQ=255)      tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑TAGGGTGC‑gggg                                                                                                                                                                                                                                         >  1:207405/1‑238 (MQ=255)      tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGa                                                                                                                                                                                                                                         <  2:684848/239‑1 (MQ=255)      tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑CT‑TAGAGTGC‑gggg                                                                                                                                                                                                                                         >  1:986196/1‑238 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATT‑AAGGTTG‑‑GGGAATCac                                                                                                                                                                                                                                    >  2:248536/1‑238 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGGATCa                                                                                                                                                                                                                                     >  1:185378/1‑238 (MQ=255)                   cAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGGTGCGGGGAATCacaaaccc                                                                                                                                                                                                                              >  1:502906/1‑231 (MQ=255)                          atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AAGGTTGC‑GGGAATCA‑AC‑AAGcccc‑ccccctt                                                                                                                                                                                                                  >  2:1081795/1‑234 (MQ=255)                                 cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AAAGGTGC‑GGGGATCA‑GC‑AAACCCC‑CCAC‑TTTTT‑CCCgg                                                                                                                                                                                                         >  1:1097435/1‑238 (MQ=255)                                     gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGTGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑TAGGTTGC‑GGGAATCA‑GC‑AAACCCG‑CCAC‑tttttacccttttt                                                                                                                                                                                                      >  1:109894/1‑229 (MQ=255)                                     gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑TAAGTGGC‑GGGGATCA‑GC‑AACACCC‑CCCC‑TTTTT‑CCCCAttt                                                                                                                                                                                                      >  1:773953/1‑238 (MQ=255)                                     gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCACATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAACAACGGAGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GA‑AAACTTGC‑AGGGACAA‑AC‑AAGCTCG‑CCCC‑ATTTA‑ACCCAtttt                                                                                                                                                                                                     >  2:1040278/1‑239 (MQ=255)                                            aaCGTCGGGATCCGTTTTTGCAGAAAACCCTAGCCAATCAGCCCCCCCCCCATGTCCCCCACCACACCAAAGAAATACATGCCGGTCGACATGGAGGCCGCCAGCACGGGGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGTTAAATCCCGGCGGTTGGGTAATCCCCCCAAGCTGGGG‑GG‑TT‑AAGGTTGG‑GGAAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTg                                                                                                                                                                                               <  2:835676/238‑1 (MQ=255)                                             cgggcGGGATCCAGTTTCGCGAAAAACCCTACCAAACCAGCCCCCCCCCCATTGTCCCCCCCACACAAAAGAGATACATTACGGTCGACATGGAGGCCGCCGGCACGGCGTGACCCAAAAACAGCAACGGTGAAATGAACGGGAAGGAAAACCCCGGCGGTTCGGTAATCCCCACCAGAAGGGG‑GG‑TT‑AGGGTGGG‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTg                                                                                                                                                                                               <  2:639637/233‑1 (MQ=255)                                                   ggATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑TA‑AAGGTGGG‑GGGAATAACGA‑AAGCGC‑‑CCCC‑TTTTA‑‑CCAATTTTCTGGTgggggcgag                                                                                                                                                                                       >  2:1249681/1‑234 (MQ=255)                                                        cattttgGCGGAAAAACATCGTAAACCACACCCCCCCCAATGTCCCCCACCACACAAAAGAAATACATCGGGGTCGACATTGAGCCCGCCAGCAGGGCGTGTCCCGCAAACAGCAGCGGTAAAAGGAAAGGGAAGGTAAATCCCAGCGGTTCGGTAACCCCCCCCAGCAGGGC‑GG‑TT‑AAGGTTGC‑GGAACTAA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTa                                                                                                                                                                                 <  2:207405/236‑1 (MQ=255)                                                         aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGGTTC‑GGGGATCA‑GC‑AAGCCCC‑CCAC‑TTTTA‑CCCCATTTTTCTGTGCGGC‑AGTGAAGTAc                                                                                                                                                                                >  1:762962/1‑239 (MQ=255)                                                           ttttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGGTGGG‑GGGAAT‑A‑GC‑AAACCCG‑CCCC‑CTTTT‑CCCGATTTTCCGGTgggggcggggaggcg                                                                                                                                                                                 >  1:639637/1‑226 (MQ=255)                                                           tttgggcaaaaaaCCTATACACTCAGCCCCCCGCCCAGGTCCCCCCCCACCCAAAAGAGAAACATTACGCGCGAAATGGAGGCCGCCAGCAGGGGGTGTACCGAAAACAGAAACGGGAAAAGGAACAGGAGGTTAAATCCCGGCGGTGCGGTAATCCCCACCAGCATGGG‑GG‑TT‑AGGGTGGG‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACAt                                                                                                                                                                              <  2:1041887/230‑1 (MQ=255)                                                            tttGGGATAAACCCGGGTCAACCCGACCCCCCCCCTTTTGCCCCCCCCCACAAAGGACATACTTTAGGGTGGACATGGAGGCCGCCAGCCCGGCGGGACCCCAAAACAGACACGGAGAAATGATAAGGAGGGAAAGTCCAAGCGGAGCGTTACTTCCCCCCGGCATGGG‑GT‑TT‑AGGGTGGG‑GGGAAGAA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTAc                                                                                                                                                                                <  2:874743/236‑1 (MQ=255)                                                             tgcgggaaaaaCCTAGTAAACCAGCCCTCCGCCCATGTCGCCCCCCACACAAAAAAGATACATTACGGTGGAAATTGAGGCCGCCGGCCCGGCGTGAACCGAAACCACAAAGGGAGAAAGAAACAGGAAGTAAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GG‑TT‑AAGGTTGG‑GGAAAACA‑GA‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCg                                                                                                                                                                            <  2:1097435/232‑1 (MQ=255)                                                                     aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AAAGGTGT‑GGGGATCA‑GC‑AAAACCG‑CCCC‑TTTTT‑TCCGATTTTTCCGTGCGGC‑GGTGGAGGACATCGCTaaaaaa                                                                                                                                                                    >  1:231935/1‑239 (MQ=255)                                                                             tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑TT‑AAGGTTGG‑GGGAA‑CA‑AC‑AAAGCCC‑CCCC‑TTTTATCCCGTTTTTCCGGTGGGGC‑GGGAAAGAAC‑TCGCCAAAGAAAAGcccc                                                                                                                                                             >  2:316051/1‑237 (MQ=255)                                                                                 tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAAGGTGG‑GGGGATCA‑GC‑AAGCCCC‑CCCC‑CTTTT‑CCCCATTTTCCGGGGCGGG‑AGGGGAGTACAACCCCTAAAAAAAGCCCCCGGc                                                                                                                                                         >  1:136869/1‑238 (MQ=255)                                                                                 tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑CT‑AAGGATTC‑GGGGATCC‑GC‑ACACCCC‑CCCC‑CTTTA‑CCCCAATTTTCCGTGCGGC‑GGTGAATGACATCGCCAACAAAATGctcccgac                                                                                                                                                         >  2:517063/1‑231 (MQ=255)                                                                                  cAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGGTGC‑GGGGATCA‑GC‑AAGCCCC‑CCCC‑TTTTA‑ACCGCTTTTCCGGGGGGGG‑AGTGAAGGACATCGCTAAAAAAATGCCCACGGCa                                                                                                                                                        >  1:968735/1‑238 (MQ=255)                                                                                          cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AAGGGTGC‑GGGGAACA‑GC‑AAACCCG‑CCCC‑TTTTT‑CCCCGTTTTTCGGTGCGGG‑CGGGGAGTACACCGCCAAAAAAAAGCCCCCCGGCCCCAAGAtt                                                                                                                                               >  1:535343/1‑239 (MQ=255)                                                                                                tGTTGCCCACCACACCAAAGAAATACATTACGGACGACATTGAGGCCGCCAGCACGGCGTGTCCCGCAAACAGCAACGGTGAAATAAACAGGAAGGAAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GT‑TT‑AAGGTTGC‑GGGAATCA‑GA‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑g                                                                                                                                        <  1:1037355/239‑1 (MQ=255)                                                                                                ggtcccccACCACCCAAAAGAGATACATTACGGTCGACATGGAGGCCGCCAGCACGGGGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGATAAATTCCAGCGGTTCGGTAACTCCCACCAGCAGGGC‑GG‑TA‑AGGGTTGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTg                                                                                                                                           <  1:1427858/232‑1 (MQ=255)                                                                                                        accacaccaAGGAGCAACATTGCGGTAGGAATGGGGGCCGCCGGCACGGCGGGTCCCGCAAACAGCAACGGTAAAATGAAAAGGAAGGTAAATCCAAGGGGTTCGGTAATCCCCACCAGCAGGGC‑GT‑TT‑AGGGTGGG‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑g                                                                                                                               <  2:1098257/239‑1 (MQ=255)                                                                                                          cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AAAGTTGC‑GGGGAACA‑GC‑AAGCCCC‑CCAC‑TTTTA‑CCCGATTTTTCGGGGGGGC‑GGGGAAGGAAAACGCTAAAAAAATGCCCACCGCACCACAGAATTTTGAGGTTACCGT‑ga                                                                                                                              >  1:1151775/1‑238 (MQ=255)                                                                                                             accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCGCGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCCCCAGCAGGGC‑GG‑TT‑AAGGTGGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCAc                                                                                                                                                            <  1:44161/207‑1 (MQ=255)                                                                                                               caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAAGGTGC‑GGGGATCC‑GC‑AAACCCG‑CCCC‑TTTTT‑CCCGCTTTTTCCGGGCGGC‑AGGGAGATACATCGCCAAAGAAACGGCCCCCGCACCCAAACTATTTTA‑GGTACCCT‑GCA‑gggga                                                                                                                       >  2:1030318/1‑237 (MQ=255)                                                                                                                      aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGGAGTACATCGCTAAAGAAATGCCCACGGCACCCAAGATTTTTGA‑GTTACCGT‑GCA‑CGG‑GAAAACCg                                                                                                                >  2:770245/1‑238 (MQ=255)                                                                                                                       tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGGTGC‑GGGGATCA‑GC‑AAACCCG‑CCCCTTTTTC‑CCCGTTTTTCCGG‑GCGGG‑AGGGGAGGAAACCGCCCAAAAAAAGCCCCCCGGCCCCAAGATTTTTTT‑GTTTCCCTGGCA‑GGG‑CAAAAccc                                                                                                                >  1:696034/1‑237 (MQ=255)                                                                                                                         cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAAGGTGC‑GGGGAACA‑GC‑AAACCCC‑CCCC‑TTTTA‑CCCCATTTTTCGGGGCGGG‑AGTGAAAGACAACGCCAAAAAAATGACCACGGCACCAAAAATTTTTTA‑GTTACCCT‑GCACGGG‑CAAAAcccct                                                                                                              >  2:537825/1‑235 (MQ=255)                                                                                                                              cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑TA‑AAGGTGCC‑GGGAACAA‑CC‑AAGCCC‑‑CCCC‑CTTTA‑TCCCATTTTTCGGGGCGGCGCGGAAAGAACAACCTCAAAAAAATCCCACCGGCCCCAAAAATTTTTTA‑GTTTCCGG‑GAG‑GGG‑ACAAACCCCCTTcccg                                                                                                        >  2:491274/1‑237 (MQ=255)                                                                                                                              cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATG‑AAGGTTG‑‑GGGAATTA‑GG‑AAACCCC‑CCAC‑CTTTT‑CCCGATTTTTCCGGGGGGG‑AGTGGAGTAAAACGCTAAAAAAAAACCCACAGCACCAAAGAATTTTGA‑GATAACGT‑GGC‑AGG‑CAAAAACGCCTTCCggg                                                                                                       >  2:1382028/1‑239 (MQ=255)                                                                                                                                gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGGTGG‑GGGGATCA‑GCAAAGCCCC‑CCCT‑TTTTA‑CCC‑ATTTTTCCGTGCGGC‑CGTGGAGGAAAACGCCAAAAAAAAGCCCCCCGCCCCCAAGAATTTTTA‑GGTACCGG‑GCA‑GGGCCAAAACCCCCTCCCGGGa                                                                                                      >  1:1048976/1‑239 (MQ=255)                                                                                                                                   gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑TAGGGTGC‑GGGAATCA‑AC‑AAAACCGCCCCC‑TTTTC‑CCCCATTTT‑CGGTGCGGG‑CGTGAAGGACAACCCTAAAAAAAAACCCACCGCACCCCAAATTTTTTA‑AGTAACGT‑TGA‑GGG‑CAAAACCCCCTCCCGGGaaaca                                                                                                  >  2:1424731/1‑236 (MQ=255)                                                                                                                                   gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AAGGTTGC‑GGGGAACA‑GC‑AAGACCC‑CCAC‑CTTTT‑CCCGATTTTTCGGTGCGGC‑AGGGAAGTAACACGCTAAAGAAAAGCCCCCCGCACCAAAGATTTTTTA‑GGTAACGG‑GCC‑AGG‑GAAAACCCCCCTCCCGGaacaa                                                                                                  >  2:234796/1‑239 (MQ=255)                                                                                                                                               gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGGTTC‑GGGGAACA‑GC‑AAGCCCC‑CCAC‑ATTTA‑ACCGATTTTTCGGTGGGGC‑AGGGAAAGAAAATGCCAAAAAAAAGCCCCAGGCCCCAAAAATTTTTTG‑GTTCCCGG‑CCG‑GGG‑CAAAACCCCCTTCCGGGAAAAACAAttctaccg                                                                                       >  2:1072488/1‑232 (MQ=255)                                                                                                                                                                 caaaaCAGCAAGGGTGAAATGAACAGGAAGGAAAATCCCGGCGGTTCGGAAATTCCCACCAGCATGGC‑GG‑TT‑AGGGTGGG‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                    <  1:994202/237‑1 (MQ=255)                                                                                                                                                                 aaaaaCAGCAGCGGTGAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGG‑GG‑TT‑AAGGTTGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                    <  1:1438721/237‑1 (MQ=255)                                                                                                                                                                  cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAAGTTTC‑GGGGATTA‑AC‑AAACCCC‑CCCC‑CTTTT‑CCCGGTTTTTCGGTGCGGC‑AGGGAAGTAAATCGCCACAAAAAAACCCCCAGCACCAAACACTTTTTA‑GTTACCCC‑TCA‑GAG‑CAAAACCGCCTGCCCGGGACAACCAATTTAGCGGAGCGGCGCtacacactca                                                                    >  2:1386523/1‑229 (MQ=255)                                                                                                                                                                     aaaGCAACGGTAAAAGGAACATGAAGGTAATTCCCAGCGGTCCGGAATCCCCCCCCAGCGTGGG‑GG‑TT‑AAGGTGGG‑GGAAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAActct                                                                <  2:136869/237‑1 (MQ=255)                                                                                                                                                                          aaCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTTC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCCCCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCCCCAACTGAACTCTTGccg                                                           >  1:309975/1‑237 (MQ=255)                                                                                                                                                                            cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGAATCA‑GC‑AAGACCC‑CCCC‑CTTTA‑CCCGATTTTCCGGTGCGGC‑CGGGAAGTACATCGCTCAAGAAATTCCCACGGCACCCAAGATTTTTTA‑GGTAACCC‑GCA‑GGG‑GAAAAACGCCTTCCGGGAAAAACGATTTCCACGGCCCGGAGCCCAAACAGaaacacttcagat                                                         >  2:446498/1‑228 (MQ=255)                                                                                                                                                                                      aaaagaaaGGAAAATCCCGGCGGTCGGGAAATTCCCCCCAGCATGGC‑GT‑TA‑AGGGGTGC‑GGGAACCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                 <  1:530942/231‑1 (MQ=255)                                                                                                                                                                                      aaaagaaaGGAAAATACCAGCGTTTGGTAAATCCCCACCAGCAGGGG‑GG‑TT‑AGGGTGGG‑GGAAAACA‑GA‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                 <  2:1419186/231‑1 (MQ=255)                                                                                                                                                                                      aaaaGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GA‑TT‑AGGGTGGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                 <  1:770245/234‑1 (MQ=255)                                                                                                                                                                                       agaggatgtaaaaTCCCAGCGGTACGGTACTCCCCACCGGCATGGC‑GG‑TT‑AAGGTTGC‑GGAAAGCA‑GC‑AGGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                              <  1:780352/229‑1 (MQ=255)                                                                                                                                                                                       aCAGGAAGTTAAAGTCCAGCGGTTCGTTAATCCCCCCCAGCATGGC‑GG‑TT‑AAGGGTGC‑GGGAATAA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                              <  1:650117/239‑1 (MQ=255)                                                                                                                                                                                           gAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑TAAGTTGC‑GGGAATCA‑GC‑AAACCCC‑CCCC‑TTTTT‑CCCCATTTTCCGGTGCGGC‑CGTGAAGTACATCGCTAAAAAAATTCCCACGGCACCAAAAATTTTTTA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACCACGATTTCAGCGGCTCCGCGCTCCAACTTAAATCCTGCAGATGCTCTGccccctac                                           >  2:530942/1‑233 (MQ=255)                                                                                                                                                                                                          ggtttggttAATCCCCCCCAGCAGGGC‑GG‑TT‑AGGGTTGC‑GGAAAACA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCt                            <  1:316051/230‑1 (MQ=255)                                                                                                                                                                                                                         caccaGAATGGC‑GG‑TA‑AAGGTGGG‑GGGAATCA‑GA‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc            <  1:9660/239‑1 (MQ=255)                                                                                                                                                                                                                         cacaaGGATGGG‑GG‑TT‑ACGCGTGT‑GGGAATGC‑GA‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTGCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATTCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACATCTGCCGGAc            <  2:107189/239‑1 (MQ=255)                                                                                                                                                                                                                                     ggg‑TT‑AGGGTTGG‑GAGAAGCA‑GC‑AGGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa  <  1:84386/237‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862425‑3862889

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