Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,739,943	T→C	intergenic (‑133/+65)	lepA ← / ← rseC	GTP binding membrane protein/RseC protein involved in reduction of the SoxR iron‑sulfur cluster

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,739,943	0	T	C	100.0%	216.6 / NA	59	intergenic (‑133/+65)	lepA/rseC	GTP binding membrane protein/RseC protein involved in reduction of the SoxR iron‑sulfur cluster
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/27);  total (32/27)

AACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAT  >  minE/1739876‑1740008                                                                    |                                                                  aaCGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATg                                                                >  1:3159570/1‑71 (MQ=255) aaCGTCTGTATTAATAGAAACGCCGAATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATg                                                                >  1:1314456/1‑71 (MQ=255)   cGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGa                                                              >  1:94273/1‑71 (MQ=255)   cGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGa                                                              >  1:2346149/1‑71 (MQ=255)         tATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATtg                                                           <  1:2486507/68‑1 (MQ=255)         tATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTgg                                                        >  1:1405857/1‑71 (MQ=255)          aTTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTgg                                                        >  1:1177708/1‑70 (MQ=255)          aTTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTgg                                                        >  1:248528/1‑70 (MQ=255)           ttAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTa                                                      <  1:938749/71‑1 (MQ=255)              aTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    <  1:1404854/70‑1 (MQ=255)              aTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    <  1:3187348/70‑1 (MQ=255)              aTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    <  1:1791393/70‑1 (MQ=255)                aGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATtgt                                                          <  1:1425978/62‑1 (MQ=255)                  aaaCGCCGCATTCTACACTACAACATTGAGGCGAGTAAATGTTCATACCGCATGGATTGTgg                                                        >  1:2539519/1‑62 (MQ=255)                  aaaCGCCGCATTCTACACTACAACATTGAGGAGAGGAAATGTTCATACCGCATGGATTGTgg                                                        >  1:187916/1‑62 (MQ=255)                        cgcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTcc                                         <  1:3187000/71‑1 (MQ=255)                         gcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAAc                                             <  1:1869849/66‑1 (MQ=255)                         gcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCt                                        >  1:2115667/1‑71 (MQ=255)                         gcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCt                                        >  1:2536234/1‑71 (MQ=255)                           aTTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGt                                           >  1:3149429/1‑66 (MQ=255)                            ttCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCg                                      >  1:3056112/1‑70 (MQ=255)                            ttCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCg                                      >  1:93204/1‑70 (MQ=255)                              ctacactacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   >  1:559236/1‑71 (MQ=255)                                 cactacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAtt                                  <  1:336597/69‑1 (MQ=255)                                 cactacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   <  1:265086/68‑1 (MQ=255)                                  actacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGt                               >  1:3019288/1‑71 (MQ=255)                                  actacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGt                               >  1:1635028/1‑71 (MQ=255)                                  actacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGt                               <  1:1430383/71‑1 (MQ=255)                                      caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    >  1:1789607/1‑46 (MQ=255)                                      caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    >  1:2795365/1‑46 (MQ=255)                                      caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    >  1:1331977/1‑46 (MQ=255)                                      caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    >  1:1340894/1‑46 (MQ=255)                                      caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    >  1:1309220/1‑46 (MQ=255)                                      caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                    >  1:794817/1‑46 (MQ=255)                                         caTTGAGGCGAGGAAATGTTCATAGCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   <  1:2753120/60‑1 (MQ=255)                                          aTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaaa                       <  1:1938219/71‑1 (MQ=255)                                          aTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaa                        >  1:2779871/1‑70 (MQ=255)                                          aTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaa                        >  1:1166269/1‑70 (MQ=255)                                           ttGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaaaa                      <  1:1428199/71‑1 (MQ=255)                                             gAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAAt                     <  1:636128/70‑1 (MQ=255)                                             gAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATg                    >  1:167109/1‑71 (MQ=255)                                             gAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATg                    >  1:2800428/1‑71 (MQ=255)                                             gAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATg                    >  1:707097/1‑71 (MQ=255)                                              aGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGc                   <  1:2866418/71‑1 (MQ=255)                                                  gAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGGAACGTCCTCGCATTTGTTATGCAAAATGcaac                >  1:1416246/1‑70 (MQ=255)                                                  gAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGcaaca               <  1:1955616/71‑1 (MQ=255)                                                    ggAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCa                         <  1:822476/59‑1 (MQ=255)                                                     gAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGcaaca               <  1:2652879/68‑1 (MQ=255)                                                         tGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGcaac                <  1:2310213/63‑1 (MQ=255)                                                         tGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGcaac                <  1:1854910/63‑1 (MQ=255)                                                             cATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGa     >  1:319531/1‑70 (MQ=255)                                                             cATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGa     >  1:3457/1‑70 (MQ=255)                                                             cATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGa     >  1:1321878/1‑70 (MQ=255)                                                               tACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAt  <  1:2608947/71‑1 (MQ=255)                                                                aCCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   <  1:2804463/37‑1 (MQ=255)                                                                aCCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   <  1:377392/37‑1 (MQ=255)                                                                aCCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   <  1:740559/37‑1 (MQ=255)                                                                aCCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   <  1:220171/37‑1 (MQ=255)                                                                aCCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                   <  1:1454862/37‑1 (MQ=255)                                                                    |                                                                  AACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAT  >  minE/1739876‑1740008

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