Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	363,067	T→C	N38S (AAC→AGC)	lacY ←	lactose permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	363,067	0	T	C	100.0%	55.9 / NA	16	N38S (AAC→AGC)	lacY	lactose permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/8);  total (8/8)

CAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGAAAGCAGACCAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGTTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACG  >  NC_000913/362949‑363195                                                                                                                       |                                                                                                                                 cagcagGTATTTGCGCAGCCCGAGTTTGTCAGAAAGCAGACCAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCg                                                                                                             >  2:320738/1‑140 (MQ=255)         aTTTGCGCAGCCCGAGTTTGTCAGAAAGCAGACCAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAAc                                                                                                     <  1:320738/140‑1 (MQ=255)             gcgcAGCCCGAGTTTGTCAGAAAGCAGACCAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGa                                                                                                 <  2:377339/140‑1 (MQ=255)                  gCCCGAGTTTGTCAGAAAGCAGACCAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAg                                                                                            <  2:118244/140‑1 (MQ=255)                   cccGAGTTTGTCAGAAAGCAGACCAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAgg                                                                                           >  1:83521/1‑140 (MQ=255)                                          cAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTaaaag                                                                    <  2:83521/140‑1 (MQ=255)                                                   gTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAAt                                                           >  1:152824/1‑140 (MQ=255)                                                                  agaACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCaa                                            >  1:280593/1‑140 (MQ=255)                                                                    aaCAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCaaaa                                          >  2:189975/1‑140 (MQ=255)                                                                         agaAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTtg                                     >  1:45815/1‑140 (MQ=255)                                                                               aGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGtttt                               <  2:152824/140‑1 (MQ=255)                                                                                    cAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTaaa                          <  1:189975/140‑1 (MQ=255)                                                                                    cAAAAATAATACCCGTATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTaaa                          <  2:179750/140‑1 (MQ=255)                                                                                                   tATCACTTTTGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGAt           >  2:372054/1‑140 (MQ=255)                                                                                                          tttGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTa    >  1:5792/1‑140 (MQ=255)                                                                                                            tGCTGATATGGCTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTAcg  <  1:72355/140‑1 (MQ=255)                                                                                                                       |                                                                                                                                 CAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGAAAGCAGACCAAACAGCGGTTGGAATAATAGCGAGAACAGAGAAATAGCGGCAAAAATAATACCCGTATCACTTTTGCTGATATGGTTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACG  >  NC_000913/362949‑363195

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