Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	361,489	G→A	N124N (AAC→AAT)	lacA ←	galactoside O‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	361,489	0	G	A	100.0%	108.8 / NA	30	N124N (AAC→AAT)	lacA	galactoside O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/18);  total (12/18)

ATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTT  >  NC_000913/361244‑361726                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              ataATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCAttttt                                                                                                                                                                                                                                          >  2:395897/1‑251 (MQ=255)                                       ataataGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGaa                                                                                                                                                                                                     >  2:2276/1‑250 (MQ=255)                                                     cccGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGt                                                                                                                                                                                      >  1:36535/1‑251 (MQ=255)                                                         gTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCa                                                                                                                                                                                   >  1:128388/1‑250 (MQ=255)                                                                      agAATACCCCAACAAGGACCGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATc                                                                                                                                                                     <  2:115466/250‑1 (MQ=255)                                                                          aagcccGCCAAGGACCGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAAGTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCGCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCg                                                                                                                                                                 <  1:386461/248‑1 (MQ=255)                                                                                         cGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCgtgt                                                                                                                                                    >  1:139934/1‑249 (MQ=255)                                                                                              gccgccACGACGTTTGGTGGAAGGTCGTTTGTGACGATCCTACCCGCGCCAATAACAGAAGTATCCCCGATGGTTACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCGCGCCATTTGTTCTCAATTCACGGTGTACAGGGTGTCACGTAACGGAGAGAGTAACGTTGGGTGCAAGCAGTGAGTAATAACCGAATGTTACCGTGCAGTCATc                                                                                                                                             <  2:39539/251‑1 (MQ=255)                                                                                               ccgccACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCg                                                                                                                                            <  2:36535/251‑1 (MQ=255)                                                                                                 gccACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGAc                                                                                                                                          <  1:374820/251‑1 (MQ=255)                                                                                                       acgTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGt                                                                                                                                    <  1:317940/251‑1 (MQ=255)                                                                                                        cgTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGtt                                                                                                                                   <  1:411868/251‑1 (MQ=255)                                                                                                           ttGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTaa                                                                                                                                 <  1:167855/250‑1 (MQ=255)                                                                                                                     gTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATGACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAAttt                                                                                                                       <  1:186604/250‑1 (MQ=255)                                                                                                                                  aTACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACAGCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCg                                                                                                         <  2:13297/251‑1 (MQ=255)                                                                                                                                      tacCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCt                                                                                                     >  1:138745/1‑251 (MQ=255)                                                                                                                                       acCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCta                                                                                                    <  1:387454/251‑1 (MQ=255)                                                                                                                                         ccGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTc                                                                                                                                                                                                                                 >  1:105822/1‑124 (MQ=255)                                                                                                                                         ccGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTc                                                                                                                                                                                                                                 <  2:105822/124‑1 (MQ=255)                                                                                                                                                   aaCAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTgg                                                                                        <  2:199330/251‑1 (MQ=255)                                                                                                                                                             cccccGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAaga                                                                              <  2:11277/250‑1 (MQ=255)                                                                                                                                                                               ggATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTc                                                            <  2:93162/251‑1 (MQ=255)                                                                                                                                                                                     ataaCCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCCAccc                                                       >  1:101035/1‑250 (MQ=255)                                                                                                                                                                                           acaTGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGtt                                                <  2:16432/251‑1 (MQ=255)                                                                                                                                                                                                  ttCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGc                                                                                                          >  2:311173/1‑186 (MQ=255)                                                                                                                                                                                                  ttCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGc                                                                                                          <  1:311173/186‑1 (MQ=255)                                                                                                                                                                                                                      ccAATCGTTATCGGAAAAGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTACATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTAACCCGGCGTTTTTCCCCACCCTGGCAAAActttcttt                     >  2:254976/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTCACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCGGCGTTTTTCCCTACCCGGGCAAACATTTCTTTTATCAGGCTTTCTCtt     >  1:202675/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaGAGTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCtt     <  2:139934/251‑1 (MQ=255)                                                                                                                                                                                                                                         agTACATCTCGCCATTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCAGGCGTTTTTCCCTACCGGGGGAAACATTTTTTTTATCAGGCTTTCTCttttt  >  2:191876/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              ATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTT  >  NC_000913/361244‑361726

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