Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,661	T→G	30.2%	pseudogene (942/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,661	0	T	G	30.2%	31.6 / 2.6	41	pseudogene (942/1617 nt)	glvC	pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (13/14);  new base G (2/12);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 4.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TAAGGTT‑GCGGG‑AATCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCG  >  NC_000913/3862424‑3862895                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   cccAATGGCGATTTGGGCAAGCATCATACCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TAAGGTT‑g                                                                                                                                                                                                                                                  <  2:1197533/239‑1 (MQ=255)      tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGGCCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TAAGGTT‑GCGGG‑a                                                                                                                                                                                                                                            <  2:551593/239‑1 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCCGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑G‑TAAGGGT‑Ggggggacta                                                                                                                                                                                                                                         >  2:758837/1‑233 (MQ=255)            tttGGGTCAGGATCATATCCGCGTGGTTGCTTAACAACGGTATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TAAGGTT‑GCGGG‑AATCAG‑c                                                                                                                                                                                                                                     <  2:1126773/239‑1 (MQ=255)            tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑T‑TAAGTTT‑G‑GGG‑GATCAG‑c                                                                                                                                                                                                                                     >  2:534495/1‑239 (MQ=255)                   cAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGCGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGATTC‑GGTAATTCCCAC‑CCGCATGGC‑GG‑TGAAAGGTT‑GGGGG‑aacaatgacacc                                                                                                                                                                                                                                 >  2:1007622/1‑228 (MQ=255)                                gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTAAGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTACCAC‑CAGAATGGC‑GG‑G‑TAAGGTT‑GCGGG‑AGTGAG‑C‑AAgacacgcccattttcacc                                                                                                                                                                                                                >  2:473457/1‑222 (MQ=255)                                gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGAAACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACCGAACCGGGGAAATGAACAGGAAGGCAAATTCACGCGGTTT‑GGTAATTACCAC‑CAGGATGGC‑GG‑A‑GTAGCGT‑TCGGG‑AATAAC‑C‑AAACCCCCC‑Acccctaccac                                                                                                                                                                                                               >  1:1001272/1‑230 (MQ=255)                                     gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTAACGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TAAAGTT‑GCCGG‑GAACAG‑C‑AAGCCCGCCAACTTTAA‑CCGGAtt                                                                                                                                                                                                           >  2:963981/1‑238 (MQ=255)                                                 cGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑T‑AAAGGTT‑G‑GGG‑GATTAG‑A‑CAACCCGCC‑CCCTTTACCCCCATTTTCC‑ggggggcag                                                                                                                                                                                             >  2:372195/1‑232 (MQ=255)                                                   ggATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TTAAGGT‑GCGGG‑GAACCG‑CAAACCCCCCC‑CCTTTTC‑CCCGATTTTCC‑‑GTGCgggccgg                                                                                                                                                                                          >  2:553702/1‑234 (MQ=255)                                                      tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑G‑TAAGGTG‑TCGGG‑GAACAG‑CAAAGCCCCCC‑CCTTTTA‑CCCGTTTTTCC‑GGGGGGGcggggag                                                                                                                                                                                        >  2:660673/1‑233 (MQ=255)                                                        cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑GAAAGGT‑GCGGG‑GATCAGCC‑AACCCCCCC‑CCTTTTC‑ACCCATTTCCC‑GTGGGGGCAGGGAAga                                                                                                                                                                                      >  1:74933/1‑238 (MQ=255)                                                        cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGTAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TAAAGTT‑TCGGG‑GATTAC‑CAAACCCCCCC‑CCTTTTT‑CCCCTTTTTCG‑GGTgggcatgcagaga                                                                                                                                                                                      >  2:178867/1‑227 (MQ=255)                                                         aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCACCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑G‑GAAGGGT‑GGGGG‑GAACAG‑A‑AAACCCCCC‑ACCTTTA‑CCCCCtttttccgcgggggggaagagaaa                                                                                                                                                                                     >  1:352805/1‑219 (MQ=255)                                                             ttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GGTT‑TAAGGTG‑GGGGG‑AATAGG‑A‑AAGCCCCCC‑CCTTTTT‑CCCCTTTTTCC‑GGGGGGGcggaggaagcaa                                                                                                                                                                                   >  2:1018607/1‑226 (MQ=255)                                                                                                                       tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑T‑TAAGGTT‑ACGGG‑AGCACG‑A‑AACCCCGCC‑CCTTTTT‑CCCGCTTTCCC‑GTTTGGGGGGGGAAGAACACCCCCAAAAAAAACCCACCGCCACCAACACATTTTT‑ATTTTCCCGGCAGGGCaaaaaca                                                                                                                       >  1:869020/1‑235 (MQ=255)                                                                                                                       tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGCGGG‑T‑AAAGGTT‑GGGGG‑AGCAAG‑C‑AGACCCCCC‑CCTTTTC‑CCCGATTTTCCGGGGGGGGCAGAG‑AGACAAACCCAAAAAAAAAGCCCAC‑GAACCAA‑AAAATTTTTAGTTACAGTGCAGGAGaaaaaaac                                                                                                                      >  1:689107/1‑234 (MQ=255)                                                                                                                                                                                aaaagaaaCAGAAAGGAAAATTCCAGC‑GTTT‑GGTAAACCCCCCACGGCAGGGC‑GG‑T‑TAGGGTT‑GCGGA‑AAACCA‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                             <  2:711228/232‑1 (MQ=255)                                                                                                                                                                                aaaaTGAACAGGAAGGTAAATTCCGGGGTTTC‑GGTAATCCCCAC‑CAGCATGGC‑GG‑T‑TAGGGTG‑GGGGG‑AATCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                             <  2:816118/237‑1 (MQ=255)                                                                                                                                                                                  aaagAACAGGAGGAGAATGTCCCGCGGTTT‑GTTCCTCCCCCC‑CAGCATGGC‑TG‑T‑TAGGGTG‑GGGGA‑AACCAA‑C‑AACGCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAACGCCCACGGCCCCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                          <  1:660673/236‑1 (MQ=255)                                                                                                                                                                                  aaTGAAGAGGAAGGTAAATTCCAGCGGTTC‑GGTAATCCCCAC‑CAGCATGGC‑GG‑T‑TAAGGTT‑GCGGA‑AATCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                          <  2:884440/239‑1 (MQ=255)                                                                                                                                                                                    tgaagagaaagataatccccgggggTTC‑GGTAATCCCCCA‑CGGATTGGT‑GG‑T‑TAAGGTG‑GCGGA‑AATCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  2:1214562/216‑1 (MQ=255)                                                                                                                                                                                      aaagggAAGGTAATTCCCGGGGGTTC‑GGAAATTCCCCC‑CAGTGTGGC‑GG‑T‑T‑AGGTT‑GGGGGAAACCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  2:1023248/233‑1 (MQ=255)                                                                                                                                                                                       agaagAAGGTGAATGCCAGCGTTTC‑GGAAACTCCACC‑CGTCAGGGC‑GG‑T‑AAAGGTG‑GCGGC‑AAGCCA‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                     <  2:1281284/235‑1 (MQ=255)                                                                                                                                                                                       aaaGGGAGGTAAATTGCAGCCGTTC‑GGTAAATCCCCC‑CAAGATGGG‑GG‑T‑TAAAGTTGGGGGG‑AACCAG‑A‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                        <  2:1164445/235‑1 (MQ=255)                                                                                                                                                                                            gaagTAAAATCCCGGCGTTCGGGTAATCCCCCC‑CGGCAGGGG‑GT‑T‑TAGGGTG‑GGGGG‑AAACAC‑A‑AATCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                   <  2:495890/234‑1 (MQ=255)                                                                                                                                                                                                            tgtcgggTACTCCCCAA‑CGGCTTGGT‑GG‑T‑GAAGGTT‑GGGGA‑AATCGG‑A‑AACCCCGCC‑ACTTTTA‑CCCGATTTGCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCtt                                  <  2:954926/233‑1 (MQ=255)                                                                                                                                                                                                                          ccc‑caGCAGGGG‑GG‑T‑TAAGGGG‑GCGGG‑AATCAG‑A‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGTCGCCTTCAACAGCTGCCGGAc                   <  2:778274/237‑1 (MQ=255)                                                                                                                                                                                                                          cac‑ccGCAGGGT‑GG‑T‑CAGGGTT‑GCGGA‑AATCCG‑A‑AGCCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTAAATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  2:329740/239‑1 (MQ=255)                                                                                                                                                                                                                           ccgaaggaTGGC‑GT‑T‑TAAGGTG‑GCGGG‑AATTAG‑C‑AACCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTGTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  2:786006/232‑1 (MQ=255)                                                                                                                                                                                                                           cccctgatgggC‑GT‑T‑TAGGGTG‑GGGGG‑AATAAG‑C‑AACCCCGCC‑CCTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  2:1018749/230‑1 (MQ=255)                                                                                                                                                                                                                           ccccggctgggC‑GG‑T‑TAAGGTT‑GCGGA‑AATCAG‑A‑AGGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACAGCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  2:660218/230‑1 (MQ=255)                                                                                                                                                                                                                            ccccatctGGG‑GG‑G‑TAAGGGT‑GGGGGAAATCAG‑A‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                   <  2:869020/232‑1 (MQ=255)                                                                                                                                                                                                                               gcggcgtggggtt‑tAGGGTG‑GCGGA‑AACCAG‑C‑ACGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                 <  2:74933/227‑1 (MQ=255)                                                                                                                                                                                                                                ggaTGGC‑TG‑A‑TAAGGAG‑GCGGG‑AACCGA‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGGGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                 <  2:1279729/234‑1 (MQ=255)                                                                                                                                                                                                                                      c‑GG‑T‑TAAGGTT‑GCGAG‑AAACAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa         <  2:437054/238‑1 (MQ=255)                                                                                                                                                                                                                                       gtg‑T‑TAGGGTG‑GGGAG‑AATCAA‑C‑AGGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa         <  2:331333/236‑1 (MQ=255)                                                                                                                                                                                                                                             tAAGGTT‑GCGGT‑AATCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc    <  2:1050781/239‑1 (MQ=255)                                                                                                                                                                                                                                              taGGGT‑GCGGGAAATCAG‑C‑AGGCCCGCC‑ACTTTTA‑CCCGATTTGCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc    <  2:640712/238‑1 (MQ=255)                                                                                                                                                                                                                                               aGGTG‑GCGGG‑AATCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCg  <  2:689107/239‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCAC‑CAGCATGGC‑GG‑T‑TAAGGTT‑GCGGG‑AATCAG‑C‑AAGCCCGCC‑ACTTTTA‑CCCGATTTTCC‑GGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAA‑AGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCG  >  NC_000913/3862424‑3862895

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