Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,660	T→G	27.2%	pseudogene (943/1617 nt)	glvC ←	pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,660	0	T	G	27.2%	58.7 / 9.5	65	pseudogene (943/1617 nt)	glvC	pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base G (15/2);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 3.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAAC  >  NC_000913/3862422‑3862889                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑TA‑GGt                                                                                                                                                                                                                                                       >  1:490334/1‑239 (MQ=255)   cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TG‑AA‑GGGTg                                                                                                                                                                                                                                                     >  1:579238/1‑239 (MQ=255)     cAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCg                                                                                                                                                                                                                                                   <  1:250532/239‑1 (MQ=255)        tGGGGATTTCGGTCAGCATCATATCCGCGTTGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCGAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAAGT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑a                                                                                                                                                                                                                                               <  2:980669/239‑1 (MQ=255)        tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGCGATACATTACGGTCAACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGACCAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGA‑GG‑GT‑AA‑AGGTGCgggg                                                                                                                                                                                                                                                >  2:1130138/1‑238 (MQ=255)        tGGCGATTTGGGTCAGCATCATACCCGCTTGGTTGCTGACCATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCGTGTTGCCCACCACACCAAAGAGATACATTACGGACGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑a                                                                                                                                                                                                                                               <  2:717549/239‑1 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑GGTTGCGGG‑GATCa                                                                                                                                                                                                                                           >  1:757423/1‑238 (MQ=255)             aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGAAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATT‑AA‑GG‑TGCggggaaca                                                                                                                                                                                                                                            >  1:639742/1‑233 (MQ=255)                  ggTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGAAACTT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAA‑GGTTGC‑GG‑AATCAgacag                                                                                                                                                                                                                                      >  1:764632/1‑234 (MQ=255)                      aGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAACTT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑CA‑GGTTGCGGG‑GATcgcaaaaaccc                                                                                                                                                                                                                                  >  2:343002/1‑228 (MQ=255)                        catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑TA‑GGGTGCGGG‑GATCAAC‑AAGCCCG‑cc                                                                                                                                                                                                                              >  1:402749/1‑238 (MQ=255)                             tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑AGTTTCGGG‑AATCAGC‑CAACCCC‑CCC‑C‑cttt                                                                                                                                                                                                                      >  1:209819/1‑239 (MQ=255)                                  gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑GGGTGCGGG‑AAACAAC‑AAACCCC‑CCA‑C‑TTTTT‑cccc                                                                                                                                                                                                                >  1:260670/1‑238 (MQ=255)                                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCCAACAGCAACGGCGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGAAATTCCCACCAGCATGGC‑GG‑GT‑AA‑GGGTGCCGG‑GATCCAC‑AAGCCCG‑CCA‑C‑ATTTA‑cccc                                                                                                                                                                                                                >  2:1061460/1‑237 (MQ=255)                                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑TT‑AA‑GG‑TGCGGG‑GAACAAG‑AAGACCG‑CCC‑C‑tttttaccca                                                                                                                                                                                                                >  1:1064912/1‑233 (MQ=255)                                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑GA‑GAGTTCTGG‑AAACAGC‑aaaactcccccactcacccca                                                                                                                                                                                                                   >  2:596831/1‑220 (MQ=255)                                       gatGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAACTT‑CCAGCGGTGCGGACATTCACACCAGCATGGCGGG‑AA‑AA‑AGTTGCGGG‑AAACAGC‑CAGCCCG‑CAC‑C‑ttttcaccgtgtt                                                                                                                                                                                                             >  2:768914/3‑230 (MQ=255)                                          gCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑CT‑AAGGGTTGCGGG‑AAT‑ACC‑AAAACCC‑CCC‑C‑CTTTA‑TCCGCTTTTCCg                                                                                                                                                                                                        >  1:430932/1‑239 (MQ=255)                                            tGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGGTGGGGG‑GAACTGCAAAACCCC‑CCC‑T‑TTTTA‑CCGGATTTCCCgt                                                                                                                                                                                                       >  1:1206568/1‑238 (MQ=255)                                                  tCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑GA‑AA‑GGTGCCGGG‑CAACACC‑AACCCCC‑CCC‑C‑TTTAC‑CCCGATTTTCCgtgccgcc                                                                                                                                                                                                 >  2:442019/1‑232 (MQ=255)                                                   cGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAA‑GGTTCGGGG‑ATAAAAC‑AAACCGC‑CCC‑C‑TTTTA‑CCCATTTTTCTGGGGGGGcg                                                                                                                                                                                                >  1:142706/1‑238 (MQ=255)                                                    gcaaagaaGTTTTGCGATAAAGCGTAGCGAATCACGCCCACGCCGGTGTTCCCCCCCACAAAGAAGAGAAAAAGCACGCGCGAGATTGAGGCCACGAGAAGGGCGTGTACCGAAAACGAAAGCGGTAAAATGAGGAGGAGGATAAATT‑CCATGGGTGCGGAATCCCCCACCAGCAGGGT‑GG‑TT‑AA‑GGTGGGGGA‑ATTCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCgg                                                                                                                                                                                                  <  2:222549/228‑1 (MQ=255)                                                        tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGGTGGGGG‑GATCAAC‑AAGACCG‑CCA‑C‑TTTTT‑CCCCATTTTTCGGTGCggggaggaag                                                                                                                                                                                          >  1:712708/1‑231 (MQ=255)                                                               ttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑TA‑GGGTTAGGG‑GATAAAG‑AAAACCG‑CCA‑A‑CTTTT‑CCCCATTTTCCGGGGCGGGCGGGAAGTaaacc                                                                                                                                                                                    >  1:595969/1‑234 (MQ=255)                                                                tgcgcaaaaaCCTAGTCAATCAACCCCCCCCCCATGTCCCCCACCACACCAAAGAGAAACATTACGGACGACAATGGGGCCGCCACCACGGCGTGCACCGAAAGCAGCAACGGTAAAATGAAGAGGAAGGTAATTC‑CCAGCGGTTGGGTACCCCCCACCAGCAGGGG‑TG‑TT‑GA‑GGTGGGGGA‑AATAACC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGc                                                                                                                                                                                 <  2:402749/233‑1 (MQ=255)                                                                  cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGGGGG‑GATCAAC‑AAAACCC‑CCA‑C‑TTTTT‑ACCGGGTTTCCCGTGCGGGAGGGAAGTA‑CACCGccc                                                                                                                                                                               >  1:926009/1‑237 (MQ=255)                                                                  cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGGAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGGTGCGGG‑GATTAGC‑AAGCCCC‑CCA‑C‑TTTTA‑TCCCGTTTTCCGTTGCGGCatgggaagacaatccca                                                                                                                                                                                >  1:347422/1‑223 (MQ=255)                                                                               tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑AGTTGCGGG‑GATCAGC‑AAGCCCC‑CCACC‑TTTTC‑CCCGTTTTTCC‑GTGCGGCCGTGAAGGA‑AAACGC‑TC‑AAGAA‑AAGCCCa                                                                                                                                                                >  1:785806/1‑238 (MQ=255)                                                                                        ccgcggcgcgtGTTGCCCCCCACACAGAGGAGATAGATTGCGTTGCATGTTGGGGCCCCCATCGCGGGCTATGCCGCAAACAGAAGGGGGGAAGTGATCAGGAAGTGAAGTC‑CCAGCGGCGCGGTATGTCCCCCGAGCAGGGG‑TT‑TT‑AA‑GGTGGAGGA‑AATCAGC‑ATCCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCC‑CGGCG‑CCaa                                                                                                                                                     <  2:112315/228‑1 (MQ=255)                                                                                            cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑AGTTGCGGG‑AATCAAC‑AAGCCCG‑CCA‑CATTTTA‑CCCGTTTTTCCGG‑GCGGCCGGGAAAGA‑CATAGC‑TA‑AAAAA‑AAACCCACCGGCA‑CAAAAGA‑t                                                                                                                                                >  1:547921/1‑239 (MQ=255)                                                                                            cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑AGATGGGGG‑GATCAGC‑AAACCCC‑CCC‑C‑CTTTA‑CCCCGTTTTTCGGTGCGGGAGTGGAGGA‑AAACGCGTA‑AGGAA‑AT‑CCCA‑Cgggccccaagcat                                                                                                                                                 >  2:291853/1‑228 (MQ=255)                                                                                            cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑AGATGCGGG‑GATCAGC‑AAACCCC‑CCC‑C‑TTTTT‑CCCGATTTTTCCTGGGGGCAGGGAAGGA‑CATCGC‑CT‑AAGAA‑AAACCCC‑CCGCC‑CCCAAGaat                                                                                                                                                >  1:1054497/1‑237 (MQ=255)                                                                                                         caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑AGGTTCGGG‑AATCAGC‑AAACCCA‑CCA‑C‑TTTTA‑CCCGATTTTTCGGTTGGGCAGTGGAGTA‑CAACAC‑TA‑AAAAA‑AAGCCCA‑CGGCT‑CCCAAGA‑CTTTTTAGGTaaccg                                                                                                                                  >  1:1003611/1‑235 (MQ=255)                                                                                                                    agagATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑AGATGCGGG‑GAACCGC‑AAACACC‑CCACC‑TTTTA‑CCC‑ACTTTTCGGTGCGGGACGGAAAGA‑CAACCA‑CA‑AAGAA‑AACCCCA‑CAGCA‑CCCAAGACTTTTTGGGTACCTATC‑CAGGGCaaa                                                                                                                       >  2:1127342/1‑239 (MQ=255)                                                                                                                        aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑GGGTTCGGG‑AGACAGC‑AAACCCC‑CCC‑C‑CTTTA‑CCCCCTTTTCCGCTGGGGCAGGGCAGTA‑CAACCC‑CC‑AAGAACAACCCCC‑CGCCA‑CC‑AAGA‑CTTTTTTGGTACCATGCCGGGGCAAAAcc                                                                                                                    >  2:178141/1‑238 (MQ=255)                                                                                                                         tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GT‑AA‑AGGTTCGGG‑GAGCAGC‑AAGCCCGCCCC‑C‑TTTTA‑CCC‑AGTTTTTGGTGCGGCAGGGAAAGA‑AAACGC‑CC‑AAACA‑ATGCCCC‑CGGCC‑ACAAAAA‑TTTTTTAGTTAccagccaagagcaaaccgc                                                                                                                   >  2:265192/1‑221 (MQ=255)                                                                                                                              tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑AGTTGCGGG‑GATCAAC‑AAACCCC‑CCC‑C‑TTTTA‑CCCGATTTTTCCGTGCGGGAGGGAAGAA‑CATCGC‑TA‑AAAAA‑ATGCCCC‑CGGCA‑CCAAAGA‑TTTTTGTGTTACCGCG‑CAGGGCaaaaaacccctccc                                                                                                            >  1:1025253/1‑229 (MQ=255)                                                                                                                              tACGGTCGACATTGAGGCCGCCAGCACGGCGTGGACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAAATCCCACCAGCATGGC‑GG‑AT‑AA‑GGTTGGGGG‑GATGCGC‑CAGCCCGCCCA‑C‑TTTAA‑CCC‑ATTTTACGGGGGGGGAGGGGAGGA‑CACTAC‑AACAAAAA‑AGCCCCA‑CGCCA‑CC‑AAGA‑GTTTTGAAAGACCGTG‑CGCGGAAAAAcacccctcac                                                                                                            >  2:404026/1‑230 (MQ=255)                                                                                                                                cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAAGT‑CCAGCGGTTCGGTAATTCCCACCAGAATGGC‑GGATA‑AA‑GGTTCCGGA‑AAT‑ACC‑AACCGCC‑CCA‑C‑CTTTACCCCGCTTTTCTGGT‑AGGCGGTTTAATAGAATCGC‑AA‑AAGAA‑ATACAC‑‑CGGCACCCACGGA‑TTTTTGAGTTATC‑TT‑CCGGGCaaaaacacccctccgg                                                                                                          >  2:503344/1‑227 (MQ=255)                                                                                                                                   tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCAACAGCATGGCGGG‑TA‑AA‑GGTTGGGGG‑AAT‑AAC‑AAAACCC‑CCC‑C‑TTTTA‑CCCCATTTTTCGGGGCGGGGGGGGAGCA‑AATCGC‑TA‑AAAAA‑ATATCCA‑CCGCC‑CACCAGA‑GTTTTGTGGTACCGCG‑CAGGGCaaaaacacactcccgggac                                                                                                       >  2:560639/1‑224 (MQ=255)                                                                                                                                      aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑GATCAGC‑AAGCCCC‑CCC‑C‑TTTTA‑CCCGATTTTTCGGTGCGGCAGGGGAGTA‑CATCGC‑TA‑AAGAA‑AAGCCCC‑CGGCA‑CCCAAGA‑TTTTTTGGTTAACGTG‑CAGGGGAAAACCGCCCTT‑CCGGaaaaa                                                                                                    >  1:767567/1‑236 (MQ=255)                                                                                                                                               ccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑TA‑AA‑GGTTG‑GGG‑AATCAGC‑AAACCCG‑CCC‑C‑CTTTA‑ACCGGTTTTTCGGGGCGGCAGGGAAGTA‑CATCGC‑TC‑AAAAA‑AAGCCCA‑CGGTC‑CCCAAGG‑TTTTTTGGGTAACGGG‑CAAGGCAAAAACC‑CCTTCCCGGGaaaaagcttttcgc                                                                                          >  2:1194356/1‑230 (MQ=255)                                                                                                                                                                   aaaaaCCGCAGCGGTAAAATGAAAGGAAAGTTAAATC‑CCAGCGGTTCGTTAATCCCCCCCAGCAGGGC‑GG‑TT‑AG‑GGTGGGGGA‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                     <  1:47914/237‑1 (MQ=255)                                                                                                                                                                    cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGGTGCGGG‑GATCAGC‑AAGCCCG‑CCA‑C‑CTTTA‑CACGCTTTTCAGGGGCGGCAGTGAAGTA‑CATCGC‑TA‑AAAAA‑AAGCCCA‑CGGCA‑CCcaaaa                                                                                                                                                  >  1:1041919/1‑163 (MQ=255)                                                                                                                                                                    cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑AGGTGCGGG‑AATCAGC‑AAACCCC‑CCA‑C‑CTTTA‑ACCGATTTTTCGGTGCGGCAGTGAAGGA‑CACCGC‑CA‑AAAAA‑AATCCCC‑CGCCC‑CCCAAGA‑TTTTTTTGTTACCGTG‑CAGGGCAAAAACC‑CCCT‑CCGGGAACAAAGATTTTAGCGGCTCGGCGCGCAAAccgcac                                                                    >  2:834236/1‑234 (MQ=255)                                                                                                                                                                     aaaCAGCAACGGTGAAGTGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATCCCCACCGGCATGGC‑GG‑TT‑AA‑GGGTGCGGA‑AATCAGC‑AGCCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAct                                                                   <  1:1127342/239‑1 (MQ=255)                                                                                                                                                                       aCAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑GATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGGTTTTTCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCa                                                                                                                                                      >  1:183287/1‑158 (MQ=255)                                                                                                                                                                       aCAGCAACGGTAAAAGAAAAAGGAAGGAAAATC‑CCAGGGGTCCGGTAATCCCCCCCAGCTGGGC‑GG‑TT‑AG‑GGTGGCGGG‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAActct                                                                 <  1:291853/239‑1 (MQ=255)                                                                                                                                                                                  gAAAGGAACGGGAAGAAAAATT‑CCAGGGGTCCGGTTATTCCCACCAGCATGGC‑GG‑TA‑AA‑GGTGGGGGG‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                       <  1:1106069/238‑1 (MQ=255)                                                                                                                                                                                   aaaTAAACAGGAAGTAAAATT‑CCAGCGGTTCGGAAATCCCCACCAGCATGGC‑GG‑TA‑AG‑GGTTGCGGA‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                       <  1:8900/237‑1 (MQ=255)                                                                                                                                                                                     aataaCAGGAAAGTAAATTCCCGGCGGTCGGGTAATCCCAACCAGCAGGGC‑GG‑TA‑AG‑GGTTGCGGA‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                    <  1:660197/236‑1 (MQ=255)                                                                                                                                                                                      tGAACGGGAGGTAAAATC‑CCAGGGTGTGGGTAATTCCCCCCAGC‑TGGC‑GGTTT‑AG‑GGTTGCGGA‑AACCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                  <  1:768931/239‑1 (MQ=255)                                                                                                                                                                                      tGAACAGGAATGAAAATC‑CCAGCGTTTCGATAATCCCCCCCAGCATGGC‑GT‑GT‑AA‑GGTTGGGGG‑AATCAAC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                  <  1:1194356/239‑1 (MQ=255)                                                                                                                                                                                        aaCAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑AATCAGC‑AAGACCC‑CCA‑C‑TTTTA‑CCCGATTTTTCGGTGCGGCAGGGGAGTA‑CATCGC‑CA‑AAAAA‑AATCCCC‑CCGCA‑CCCAAAA‑TTTTTGAGGTACCGTG‑CAGGGCAAAAAAC‑CCTT‑CCGGGAACAACCATTTTAGCGGCCCGGGGCTCAAACGGAACTCtttccgcagcggggccc                                                 >  2:665205/1‑223 (MQ=255)                                                                                                                                                                                           aGGAAGTAAATTT‑CCAGCGTTGCGGTAATCCCCACCAGCTGGGC‑GT‑TT‑AA‑GGTTGTGGA‑ACTCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                             <  1:106471/239‑1 (MQ=255)                                                                                                                                                                                                      atcccAGCGGTGCGTTACTCCCCCCCAGAAGGGC‑GG‑TT‑AG‑GGTTGCGGA‑AATCAGC‑ACGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                   <  1:740045/236‑1 (MQ=255)                                                                                                                                                                                                         ccagtggtttgtaacatCCCACCAGCAGGGC‑GG‑TC‑AG‑TGTGGGGGG‑AACAGGA‑AACCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccg                                <  1:768914/223‑1 (MQ=255)                                                                                                                                                                                                                    ggacaccaccccacggcgggggg‑TT‑AA‑GGGTGGGGGTAACCGGC‑AACCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACa                       <  2:260670/218‑1 (MQ=255)                                                                                                                                                                                                                            cccccgCATGGC‑GG‑TT‑AA‑GGTGGCGGG‑AACCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc             <  1:665205/234‑1 (MQ=255)                                                                                                                                                                                                                              cgagatgggC‑GG‑TT‑AA‑GGTTGCGGG‑AAGCAGC‑AGCCACG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa           <  1:657388/232‑1 (MQ=255)                                                                                                                                                                                                                              ccaGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa           <  1:721084/239‑1 (MQ=255)                                                                                                                                                                                                                                ccaggttgctg‑TT‑TA‑GGCTGGGGGAAATCAGA‑ACGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAATTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa           <  1:1182449/229‑1 (MQ=255)                                                                                                                                                                                                                                     ggC‑GG‑TT‑AA‑GGTGGCGGG‑AATAAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATa     <  1:442019/238‑1 (MQ=255)                                                                                                                                                                                                                                       c‑GG‑TT‑AA‑GGTTGGGGG‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa   <  1:842811/238‑1 (MQ=255)                                                                                                                                                                                                                                        ggg‑TT‑AA‑GGTGGCGGA‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa   <  2:1025253/237‑1 (MQ=255)                                                                                                                                                                                                                                         gg‑TT‑AA‑GGTTGCGGA‑AATCAGA‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAAc  <  1:834236/238‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TT‑AA‑GGTTGCGGG‑AATCAGC‑AAGCCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTA‑CATCGC‑TA‑AAGAA‑ATGCCCA‑CGGCA‑CCAAAGA‑TTTTTGAGTTACCGTG‑CAGGGCAAAACCG‑CCTT‑CCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAAC  >  NC_000913/3862422‑3862889

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