Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,676	G→C	32.7%	pseudogene (927/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,676	0	G	C	32.7%	24.5 / 7.9	44	pseudogene (927/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 G (6/20);  new base C (4/9);  total (14/30)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862455‑3862904                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTAAGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTACCAC‑CAGAATGGC‑GG‑GTAAGGTT‑GCGGG‑AGTGAG‑C‑AAgacacgcccattttcacc                                                                                                                                                                                                                          >  2:473457/1‑222 (MQ=255) gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGAAACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACCGAACCGGGGAAATGAACAGGAAGGCAAATTCACGCGGTTT‑GGTAATTACCAC‑CAGGATGGC‑GG‑AGTAGCGT‑TCGGG‑AATAAC‑C‑AAA‑CCCCCC‑Acccctaccac                                                                                                                                                                                                                        >  1:1001272/1‑230 (MQ=255)      gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTAACGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAAGTT‑GCCGG‑GAACAG‑C‑AAG‑CCCGCCAACTTTAA‑CCGGAtt                                                                                                                                                                                                                    >  2:963981/1‑238 (MQ=255)                  cGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑TAAAGGTT‑G‑GGG‑GATTAG‑A‑CAA‑CCCGCC‑CCCTTTACCCCCATTTTCC‑ggggggcag                                                                                                                                                                                                      >  2:372195/1‑232 (MQ=255)                    ggATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTTAAGGT‑GCGGG‑GAACCG‑CAAAC‑CCCCCC‑CCTTTTC‑CCCGATTTTCC‑‑GTGCgggccgg                                                                                                                                                                                                   >  2:553702/1‑234 (MQ=255)                       tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑GTAAGGTG‑TCGGG‑GAACAG‑CAAAG‑CCCCCC‑CCTTTTA‑CCCGTTTTTCC‑GGGGGGGcggggag                                                                                                                                                                                                 >  2:660673/1‑233 (MQ=255)                         cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TGAAAGGT‑GCGGG‑GATCAGCC‑AAC‑CCCCCC‑CCTTTTC‑ACCCATTTCCC‑GTGGGGGCAGGGAAga                                                                                                                                                                                               >  1:74933/1‑238 (MQ=255)                         cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGTAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAAGTT‑TCGGG‑GATTAC‑CAAAC‑CCCCCC‑CCTTTTT‑CCCCTTTTTCG‑GGTgggcatgcagaga                                                                                                                                                                                               >  2:178867/1‑227 (MQ=255)                          aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCACCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑GGAAGGGT‑GGGGG‑GAACAG‑A‑AAA‑CCCCCC‑ACCTTTA‑CCCCCtttttccgcgggggggaagagaaa                                                                                                                                                                                              >  1:352805/1‑219 (MQ=255)                              ttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GGTTTAAGGTG‑GGGGG‑AATAGG‑A‑AAG‑CCCCCC‑CCTTTTT‑CCCCTTTTTCC‑GGGGGGGcggaggaagcaa                                                                                                                                                                                            >  2:1018607/1‑226 (MQ=255)                                                                                        tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑TTAAGGTT‑ACGGG‑AGCACG‑A‑AAC‑CCCGCC‑CCTTTTT‑CCCGCTTTCCC‑GTTTGGGGGGGGAAGAACACCCCCAAAAAAAACCCACCGCCACCAACACATTTTT‑ATTTTCCCGGCAGGGCaaaaaca                                                                                                                                >  1:869020/1‑235 (MQ=255)                                                                                        tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑TAAAGGTT‑GGGGG‑AGCAAG‑C‑AGA‑CCCCCC‑CCTTTTC‑CCCGATTTTCCGGGGGGGGCAGAG‑AGACAAACCCAAAAAAAAAGCCCAC‑GAACCAA‑AAAATTTTTAGTTACAGTGCAGGAGaaaaaaac                                                                                                                               >  1:689107/1‑234 (MQ=255)                                                                                                                                                 aaaagaaaCAGAAAGGAAAATTCCAGC‑GTTT‑GGTAAACCCCCCACGGCAGGGC‑GG‑TTAGGGTT‑GCGGA‑AAACCA‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                      <  2:711228/232‑1 (MQ=255)                                                                                                                                                 aaaaTGAACAGGAAGGTAAATTCCGGGGTTTC‑GGTAATCCCCAC‑CAGCATGGC‑GG‑TTAGGGTG‑GGGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                      <  2:816118/237‑1 (MQ=255)                                                                                                                                                   aaagAACAGGAGGAGAATGTCCCGCGGTTT‑GTTCCTCCCCCC‑CAGCATGGC‑TG‑TTAGGGTG‑GGGGA‑AACCAA‑C‑AAC‑GCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAACGCCCACGGCCCCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                   <  1:660673/236‑1 (MQ=255)                                                                                                                                                   aaTGAAGAGGAAGGTAAATTCCAGCGGTTC‑GGTAATCCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGA‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                   <  2:884440/239‑1 (MQ=255)                                                                                                                                                     tgaagagaaagataatccccgggggTTC‑GGTAATCCCCCA‑CGGATTGGT‑GG‑TTAAGGTG‑GCGGA‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                 <  2:1214562/216‑1 (MQ=255)                                                                                                                                                       aaagggAAGGTAATTCCCGGGGGTTC‑GGAAATTCCCCC‑CAGTGTGGC‑GG‑TT‑AGGTT‑GGGGGAAACCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                 <  2:1023248/233‑1 (MQ=255)                                                                                                                                                        agaagAAGGTGAATGCCAGCGTTTC‑GGAAACTCCACC‑CGTCAGGGC‑GG‑TAAAGGTG‑GCGGC‑AAGCCA‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                              <  2:1281284/235‑1 (MQ=255)                                                                                                                                                        aaaGGGAGGTAAATTGCAGCCGTTC‑GGTAAATCCCCC‑CAAGATGGG‑GG‑TTAAAGTTGGGGGG‑AACCAG‑A‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                 <  2:1164445/235‑1 (MQ=255)                                                                                                                                                             gaagTAAAATCCCGGCGTTCGGGTAATCCCCCC‑CGGCAGGGG‑GT‑TTAGGGTG‑GGGGG‑AAACAC‑A‑AAT‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                            <  2:495890/234‑1 (MQ=255)                                                                                                                                                                             tgtcgggTACTCCCCAA‑CGGCTTGGT‑GG‑TGAAGGTT‑GGGGA‑AATCGG‑A‑AAC‑CCCGCC‑ACTTTTA‑CCCGATTTGCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCtt                                           <  2:954926/233‑1 (MQ=255)                                                                                                                                                                                           ccc‑caGCAGGGG‑GG‑TTAAGGGG‑GCGGG‑AATCAG‑A‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGTCGCCTTCAACAGCTGCCGGAc                            <  2:778274/237‑1 (MQ=255)                                                                                                                                                                                           cac‑ccGCAGGGT‑GG‑TCAGGGTT‑GCGGA‑AATCCG‑A‑AGC‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTAAATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  2:329740/239‑1 (MQ=255)                                                                                                                                                                                            ccgaaggaTGGC‑GT‑TTAAGGTG‑GCGGG‑AATTAG‑C‑AAC‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTGTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  2:786006/232‑1 (MQ=255)                                                                                                                                                                                            cccctgatgggC‑GT‑TTAGGGTG‑GGGGG‑AATAAG‑C‑AAC‑CCCGCC‑CCTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  2:1018749/230‑1 (MQ=255)                                                                                                                                                                                            ccccggctgggC‑GG‑TTAAGGTT‑GCGGA‑AATCAG‑A‑AGG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACAGCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  2:660218/230‑1 (MQ=255)                                                                                                                                                                                             ccccatctGGG‑GG‑GTAAGGGT‑GGGGGAAATCAG‑A‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  2:869020/232‑1 (MQ=255)                                                                                                                                                                                                gcggcgtggggtttAGGGTG‑GCGGA‑AACCAG‑C‑ACG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                          <  2:74933/227‑1 (MQ=255)                                                                                                                                                                                                 ggaTGGC‑TG‑ATAAGGAG‑GCGGG‑AACCGA‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGGGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                          <  2:1279729/234‑1 (MQ=255)                                                                                                                                                                                                       c‑GG‑TTAAGGTT‑GCGAG‑AAACAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                  <  2:437054/238‑1 (MQ=255)                                                                                                                                                                                                        gtg‑TTAGGGTG‑GGGAG‑AATCAA‑C‑AGG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                  <  2:331333/236‑1 (MQ=255)                                                                                                                                                                                                        ggggtgtaaagtgtgcgaaaaccAG‑C‑AAC‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑AGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGAt                     <  2:850970/216‑1 (MQ=255)                                                                                                                                                                                                           ggtgggcggtgcgtacaaccAG‑C‑AAC‑CCCGGC‑CCTTTGA‑CCCGAATTTCC‑GGTGCGGCAGTGAAGTACAGCGGTAAAGAAATGCCCCCCGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCGCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                  <  1:701700/219‑1 (MQ=255)                                                                                                                                                                                                           gggtaggggtgggagaaatgcg‑C‑AAC‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                  <  2:889363/217‑1 (MQ=255)                                                                                                                                                                                                             tAAGGTT‑GCGGT‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc             <  2:1050781/239‑1 (MQ=255)                                                                                                                                                                                                              taGGGT‑GCGGGAAATCAG‑C‑AGG‑CCCGCC‑ACTTTTA‑CCCGATTTGCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc             <  2:640712/238‑1 (MQ=255)                                                                                                                                                                                                               tagtggtgggg‑AATCAA‑A‑AAGCCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GTTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc             <  1:553702/231‑1 (MQ=255)                                                                                                                                                                                                               aGGTG‑GCGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCg           <  2:689107/239‑1 (MQ=255)                                                                                                                                                                                                                      cGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTCCATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAAATTTCCGTaactt     >  2:226305/1‑235 (MQ=255)                                                                                                                                                                                                                      cGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg     >  2:383194/1‑239 (MQ=255)                                                                                                                                                                                                                      cGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACACCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCACAATGAACTCTTGACGATGCCGCGCCCAGTACATCCGCACGCCGCCTTCAACACCCTCGGCACCCAAGATAAAACGTCCGTAAATg     >  1:621296/1‑239 (MQ=255)                                                                                                                                                                                                                           aaccAG‑A‑AAG‑CCCCCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg  <  2:813525/235‑1 (MQ=255)                                                                                                                                                                                                                            aacagac‑AAC‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg  <  2:21322/232‑1 (MQ=255)                                                                                                                                                                                                                            aaaccac‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg  <  2:436397/232‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862455‑3862904

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