Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,788,645	T→C	Y71H (TAC→CAC)	aroH →	3‑deoxy‑7‑phosphoheptulonate synthase, Trp‑sensitive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,788,645	0	T	C	100.0%	116.2 / NA	33	Y71H (TAC→CAC)	aroH	3‑deoxy‑7‑phosphoheptulonate synthase, Trp‑sensitive
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/12);  total (21/12)

TAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGC  >  NC_000913/1788514‑1788776                                                                                                                                    |                                                                                                                                    tAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCAcc                                                                                                                             >  1:1034366/1‑140 (MQ=255)         ggCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCa                                                                                                                     >  2:423511/1‑140 (MQ=255)         ggCGCCGCCCCCCATGTCCCCGCCTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCa                                                                                                                     <  1:317964/140‑1 (MQ=255)                ccACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTg                                                                                                              >  2:666323/1‑140 (MQ=255)                      aTGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAac                                                                                                        <  1:666323/140‑1 (MQ=255)                       tGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAacc                                                                                                       >  1:528395/1‑140 (MQ=255)                                       cgcAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCt                                                                                       <  2:528395/140‑1 (MQ=255)                                                                   aaGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGaa                                                           >  1:461798/1‑140 (MQ=255)                                                                      aTAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGaaaaa                                                        >  2:625745/1‑140 (MQ=255)                                                                         aGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCa                                                     <  2:354906/140‑1 (MQ=255)                                                                           cGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACg                                                   >  2:809835/1‑140 (MQ=255)                                                                           cGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACg                                                   >  1:838008/1‑140 (MQ=255)                                                                             aCTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGTTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGaa                                                 <  1:423511/140‑1 (MQ=255)                                                                             aCTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGaa                                                 >  2:750720/1‑140 (MQ=255)                                                                             aCTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGaa                                                 <  2:461798/140‑1 (MQ=255)                                                                             aCTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGaa                                                 <  2:305051/140‑1 (MQ=255)                                                                                  tGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACtgtt                                            <  1:809835/140‑1 (MQ=255)                                                                                            tGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGa                                  <  1:726139/140‑1 (MQ=255)                                                                                            tGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGa                                  <  2:244862/140‑1 (MQ=255)                                                                                              gCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGaaa                                >  1:326153/1‑140 (MQ=255)                                                                                                   tGCTCGATCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACt                           >  2:74389/1‑140 (MQ=255)                                                                                                         aTCCACGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAAtctc                     >  2:676951/1‑140 (MQ=255)                                                                                                             aCGATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGAt                 >  2:815431/1‑140 (MQ=255)                                                                                                               gATCTCACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATcc               >  1:189011/1‑140 (MQ=255)                                                                                                                    cACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAtt          >  2:591329/1‑140 (MQ=255)                                                                                                                    cACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAtt          >  2:37017/1‑140 (MQ=255)                                                                                                                    cACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAtt          >  1:949164/1‑140 (MQ=255)                                                                                                                    cACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAtt          >  1:652753/1‑140 (MQ=255)                                                                                                                    cACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAtt          >  1:586808/1‑140 (MQ=255)                                                                                                                    cACCGCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAtt          >  1:583912/1‑140 (MQ=255)                                                                                                                        gCTGCAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTaaa      <  2:586808/140‑1 (MQ=255)                                                                                                                            cAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGc  >  2:1024378/1‑140 (MQ=255)                                                                                                                            cAATGGAGCACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGc  <  1:676951/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    TAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGC  >  NC_000913/1788514‑1788776

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