Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	460,284	G→T	A725S (GCA→TCA)	lon →	DNA‑binding ATP‑dependent protease La

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	460,284	0	G	T	100.0%	216.3 / NA	66	A725S (GCA→TCA)	lon	DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (36/30);  total (36/30)

GCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGC  >  W3110S.gb/460238‑460337                                               |                                                      gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCgcggc                                         >  1:1027507/1‑61 (MQ=255) gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCgcggc                                         >  1:890347/1‑61 (MQ=255) gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCgcggc                                         >  1:4643996/1‑61 (MQ=255) gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCgcggc                                         >  1:3632868/1‑61 (MQ=255) gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCgcggc                                         >  1:2795027/1‑61 (MQ=255) gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCgcggc                                         >  1:2698846/1‑61 (MQ=255) gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCgcggc                                         >  1:2207825/1‑61 (MQ=255)  cGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGcggcgg                                       >  1:1490749/1‑62 (MQ=255)  cGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGcggcgg                                       >  1:416434/1‑62 (MQ=255)  cGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGcggcgg                                       >  1:1164588/1‑62 (MQ=255)  cGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGcggcgg                                       >  1:3062173/1‑62 (MQ=255)         aGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTaa                                 >  1:4080191/1‑61 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGa                                     <  1:320114/50‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGa                                     <  1:4415410/50‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGa                                     <  1:2845680/50‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGa                                     <  1:1625723/50‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGa                                     <  1:2897547/50‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGa                                     <  1:2600151/50‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGa                                     <  1:2922485/50‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTaaaa                               <  1:734085/56‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTaaaa                               <  1:4207502/56‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGt                           >  1:3801109/1‑60 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTg                          >  1:1640028/1‑61 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTg                          >  1:4326677/1‑61 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGc                         <  1:1876736/62‑1 (MQ=255)                gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGc                         <  1:3815648/62‑1 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:3438573/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:3171222/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:3612780/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:1193507/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:1745741/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:1651361/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:4277088/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:816935/1‑62 (MQ=255)                    aTCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                     >  1:1025454/1‑62 (MQ=255)                       ggtggtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTcc                  >  1:4172905/1‑62 (MQ=255)                       ggtggtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTcc                  >  1:1205092/1‑62 (MQ=255)                       ggtggtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTcc                  >  1:2349042/1‑62 (MQ=255)                       ggtggtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTcc                  >  1:4015726/1‑62 (MQ=255)                           gtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:4266515/62‑1 (MQ=255)                           gtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:838488/62‑1 (MQ=255)                           gtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:1391177/62‑1 (MQ=255)                           gtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:1517904/62‑1 (MQ=255)                           gtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:2156000/62‑1 (MQ=255)                           gtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:2751991/62‑1 (MQ=255)                            tttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGtt               >  1:127942/1‑60 (MQ=255)                            tttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGtt               >  1:4616350/1‑60 (MQ=255)                            tttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              >  1:3022566/1‑61 (MQ=255)                             ttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCg             >  1:3123759/1‑61 (MQ=255)                             ttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCg             >  1:1464452/1‑61 (MQ=255)                             ttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCg             >  1:1113244/1‑61 (MQ=255)                             ttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCg             >  1:4297335/1‑61 (MQ=255)                             ttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCg             >  1:2793294/1‑61 (MQ=255)                                  agaaaAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAAt        <  1:1208428/61‑1 (MQ=255)                                  agaaaAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAAt        <  1:1526858/61‑1 (MQ=255)                                  agaaaAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAAt        <  1:3144974/61‑1 (MQ=255)                                  agaaaAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAAt        <  1:1503961/61‑1 (MQ=255)                                  agaaaAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAAt        <  1:3202677/61‑1 (MQ=255)                                       aaCTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAAcgc  <  1:1550601/62‑1 (MQ=255)                                           ccTGTCAGCGCATCGCGGCGGGATTAAAACAGTGc                         <  1:3722618/35‑1 (MQ=255)                                           ccTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:2721876/46‑1 (MQ=255)                                           ccTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:3728423/46‑1 (MQ=255)                                           ccTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:3752417/46‑1 (MQ=255)                                           ccTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:432529/46‑1 (MQ=255)                                           ccTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:1543167/46‑1 (MQ=255)                                           ccTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTc              <  1:3134859/46‑1 (MQ=255)                                               |                                                      GCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGC  >  W3110S.gb/460238‑460337

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