Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,677	G→C	26.2%	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,677	0	G	C	26.2%	25.3 / 2.9	37	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/18);  new base C (3/6);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862456‑3862905                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          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)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AATCAG‑C‑AAG‑CCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862456‑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: ‑