Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,673	C→A	30.8%	pseudogene (930/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,673	0	C	A	30.8%	21.8 / 4.3	43	pseudogene (930/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 C (11/19);  new base A (4/9);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-02

TTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862435‑3862904                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              tttGGGTCAGGATCATATCCGCGTGGTTGCTTAACAACGGTATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AATCAG‑c                                                                                                                                                                                                                                               <  2:1126773/239‑1 (MQ=255) tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑TTAAGTTT‑G‑GGG‑GATCAG‑c                                                                                                                                                                                                                                               >  2:534495/1‑239 (MQ=255)                     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)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              TTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  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 < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑