Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,670	T→A	18.1%	pseudogene (934/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,670	0	T	A	18.1%	75.1 / 5.3	59	pseudogene (934/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 (25/18);  new base A (5/4);  total (33/26)
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.49e-01

ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAA  >  NC_000913/3862435‑3862904                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     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)                                                                                                                                                                                                                                      gggC‑GGGAAACA‑GC‑AAGGCCG‑CCAC‑TTTTA‑CCCGGTTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATGTTTGA‑GTTACCGT‑GCA‑GTG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt         <  2:1369119/235‑1 (MQ=255)                                                                                                                                                                                                                                         c‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg    >  2:1442183/1‑239 (MQ=255)                                                                                                                                                                                                                                            ggAATCACGC‑AGGCCCG‑CCAC‑TTGTT‑CCCGATTTTCCGGGGCGGC‑GGTGAAGTACATCGCTGAACAAATGCCCACGGCACCAAAGAGGTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCTGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCTCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg    <  2:1285724/238‑1 (MQ=255)                                                                                                                                                                                                                                             gAATCA‑GC‑AACCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa  <  1:472038/238‑1 (MQ=255)                                                                                                                                                                                                                                             aaaTAA‑GC‑ACCCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa  <  1:234796/237‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGGTTGC‑GGGAATCA‑GC‑AAGCCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGC‑AGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGA‑GTTACCGT‑GCA‑GGG‑CAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAA  >  NC_000913/3862435‑3862904

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