Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,674	C→A	35.7%	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,674	0	C	A	35.7%	12.5 / 4.8	35	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 (9/14);  new base A (3/9);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-02

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

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