Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,714,239	C→A	R41L (CGC→CTC)	nirB ←	nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,714,239	0	C	A	100.0%	126.8 / NA	38	R41L (CGC→CTC)	nirB	nitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/28);  total (10/28)

CCAGCGACAGCTCTTCGGCGGTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGCGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGCCGCATCAGATTTATCAAG  >  minE/2714171‑2714294                                                                     |                                                        ccAGCGACAGCTCTTCGGCGGTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATgagc                                                       <  1:196536/71‑1 (MQ=255)       aCAGCTCTTCGGCGGTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGttct                                                 <  1:2770433/71‑1 (MQ=255)         aGCTCTTCGGCGGTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGttctt                                                >  1:90391/1‑70 (MQ=255)              ttCGGCGGTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGtt                                                   <  1:1544082/62‑1 (MQ=255)                cggcggTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGaa                                         >  1:2826058/1‑70 (MQ=255)                 ggcggTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGaaaa                                       <  1:2712868/71‑1 (MQ=255)                   cggTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGCTCATAAGCCATGAGCCGTTCTTCACAGAAAACg                                     >  1:2371737/1‑71 (MQ=255)                     gtgtGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATgagc                                                       <  1:1019373/51‑1 (MQ=255)                     gtgtGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATgagc                                                       <  1:553342/51‑1 (MQ=255)                     gtgtGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATgagc                                                       <  1:2333182/51‑1 (MQ=255)                     gtgtGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATgagc                                                       <  1:1255845/51‑1 (MQ=255)                     gtgtGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACgg                                    >  1:905096/1‑70 (MQ=255)                     gtgtGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGt                                   <  1:112007/71‑1 (MQ=255)                      tgtgGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTa                                  <  1:97359/71‑1 (MQ=255)                          gtgAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAata                               <  1:2861409/70‑1 (MQ=255)                          gtgAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAata                               <  1:764034/70‑1 (MQ=255)                              gagaAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            >  1:1600775/1‑69 (MQ=255)                                gaAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGt                        <  1:2800770/71‑1 (MQ=255)                                gaAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGT‑ATATCAAAGtt                       <  1:1434007/71‑1 (MQ=255)                                 aaGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGtt                       <  1:2901114/71‑1 (MQ=255)                                 aaGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGtt                       <  1:1208407/71‑1 (MQ=255)                                 aaGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGt                        <  1:141588/70‑1 (MQ=255)                                  aGTAAGACGAGAGGTGTACGCGGTCATAAGCGATgagc                                                       <  1:1656764/38‑1 (MQ=255)                                  aGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTg                      <  1:488243/71‑1 (MQ=255)                                   gTAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTgg                     <  1:2373443/71‑1 (MQ=255)                                    tAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGaaaa                                       <  1:647060/52‑1 (MQ=255)                                    tAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGaaaa                                       <  1:3069442/52‑1 (MQ=255)                                    tAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGaaaa                                       <  1:2840315/52‑1 (MQ=255)                                    tAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGaaaa                                       <  1:2886508/52‑1 (MQ=255)                                    tAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGc                    <  1:619664/71‑1 (MQ=255)                                    tAAGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGc                    <  1:1794643/71‑1 (MQ=255)                                     aaGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            >  1:3181377/1‑62 (MQ=255)                                     aaGACGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            >  1:1127251/1‑62 (MQ=255)                                         cGAGAGGTGTACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGCCGCa                <  1:1029774/70‑1 (MQ=255)                                                  tACGCTGTCATAAGCGATGATCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGCCGCATCAGATTTa       <  1:1035143/70‑1 (MQ=255)                                                  tACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            <  1:1952236/49‑1 (MQ=255)                                                  tACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            <  1:192921/49‑1 (MQ=255)                                                  tACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            <  1:890600/49‑1 (MQ=255)                                                  tACGCGGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            <  1:1378122/49‑1 (MQ=255)                                                    cgcgGTCATAAGCGATGAGCGGTTCTTCACAGAAAACg                                     >  1:2461711/1‑38 (MQ=255)                                                    cgcgGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCa                            <  1:2701743/47‑1 (MQ=255)                                                     gcgGTCATAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGCCGCATCAGATTTATc     >  1:2098310/1‑69 (MQ=255)                                                           aTAAGCGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGCCGCATCAGATTTATCaag  >  1:1832242/1‑66 (MQ=255)                                                                cGATGAGCGGTTCTTCACAGAAAACGGTAATATCAAAGt                        <  1:965736/39‑1 (MQ=255)                                                                     |                                                        CCAGCGACAGCTCTTCGGCGGTGTGGTGAGAGAAGTAAGACGAGAGGTGTACGCGGTCATAAGCGATGCGCGGTTCTTCACAGAAAACGGTAATATCAAAGTTGGCCGCATCAGATTTATCAAG  >  minE/2714171‑2714294

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: ‑