Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	70,659	G→A	W91* (TGG→TGA)	araC →	DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	70,659	0	G	A	100.0%	77.8 / NA	22	W91* (TGG→TGA)	araC	DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/13);  total (9/13)

GATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCA  >  NC_000913/70421‑70905                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       gATACTCGTTTAACGCCCTTCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCGCTATT‑CGCGTTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGTAGAGATGCATCACTACGTTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTgg                                                                                                                                                                                                                                            <  1:17119/251‑1 (MQ=255) gATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTAACCATT‑CGCGGTAAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTATGCCGACCGGGTGATATTTGGCGGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTgg                                                                                                                                                                                                                                            <  1:10090/251‑1 (MQ=255)               gCCCATCTGGTGGCGGGTTTAACGCCGCTTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGACTGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGGGAAAAATCAGGGTCGAGAATTTGTCTGCCGCCCGGGGGGTATTTTGCTGTTCCCGCCCGGAGAGATTCATCACTACTGGCGTCATCCGGAGGCTCGGGAATGATATCACCAGTGGGTTTACTTTCGTcc                                                                                                                                                                                                                              >  2:199633/1‑251 (MQ=255)                cccATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCcg                                                                                                                                                                                                                             >  2:83076/1‑251 (MQ=255)                   aTCTGGTGGCGGGGTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGGGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCcgcgc                                                                                                                                                                                                                          >  2:165081/1‑251 (MQ=255)                     ttggtggCGGGTTTAACGCCGATTGAGGCCAACGGTTAGCTCGATTTTTTTATCGACCGACCGCTGGGAATTAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCcgctcg                                                                                                                                                                                                                         <  2:206620/249‑4 (MQ=255)                                                      ggTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGAGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCAGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTccat                                                                                                                                                                                       >  2:115312/1‑248 (MQ=255)                                                                    ttttATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCCCCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACg                                                                                                                                                                          <  1:83076/250‑1 (MQ=255)                                                                        aTCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGttt                                                                                                                                                                     <  1:238834/251‑1 (MQ=255)                                                                          cacccgaccgCTGGTAATGAAAGTTTATATGCTCAATCTC‑CCATTTCGCGGTCAGTGGGTGGTGAAAAATCGGGGACGAGAATTTGTTTGCCGCCCGGGTAATATTTTGCTTTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTTGCCGTCAATATTTGCCTGTACGGGtttct                                                                                                                                                                   <  1:51916/248‑1 (MQ=255)                                                                            accgaccgCTGGAAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGTAAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttcttt                                                                                                                                                                 <  1:207742/251‑1 (MQ=255)                                                                             ccgaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctttc                                                                                                                                                                >  1:169240/1‑251 (MQ=255)                                                                                                         ctcAATCTCACCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                                                                                                                    >  2:195952/1‑251 (MQ=255)                                                                                                                      atgcgcgGTCAGGGGGGGGTGAAAATCCGGGGAGGAGAATTTGTCTGCCGACCGGGTGATATTTTGCGGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                                                                                                                        <  2:219478/248‑1 (MQ=255)                                                                                                                                               aTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCa                                                                                                >  1:25442/1‑250 (MQ=255)                                                                                                                                                    ggACGAGAATTTGTCTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCCCCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGa                                                                                          <  1:144663/251‑1 (MQ=255)                                                                                                                                                               tgtgtgccaacggggagATATTTTGGTGTCCCCTCCAGGAGAGATTCATCACGACGGTCGTCATCCGGAGGCTCGCGATTGATATCACCAGTGGGTTTACTTTCGGCCGCGCGCCTACTGGCATGAATGGCTTAGCTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTGTGGGCAAATCATTAACGCCCTGCAAGTGGAAGGGCGCTAt                                                                                <  2:136412/234‑1 (MQ=255)                                                                                                                                                                        aCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAgctgc                                                                      >  2:237335/1‑251 (MQ=255)                                                                                                                                                                         ccGGGTGATATTTGGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACAGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAgctgct                                                                     <  2:164906/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAAGCCGGGCACGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGTAATTGGTACTGAGGCGCATGGAAGCGATTAACGAgg             >  2:7727/1‑250 (MQ=255)                                                                                                                                                                                                                                            aaTGATATCACCAGTGGGTTTACTTTCGTCCGCTCGCCTTCTGTCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTGCTGCGGCGCATGGAAGCGATTAACGAGTCGCtccatcca  <  1:120124/251‑1 (MQ=255)                                                                                                                                                                                                                                            aaTGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCtccatcca  <  1:113507/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       GATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATT‑CGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCA  >  NC_000913/70421‑70905

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