Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,823,561	A→C	intergenic (+306/+652)	nrfG → / ← pheU	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,823,561	0	A	C	100.0%	304.8 / NA	82	intergenic (+306/+652)	nrfG/pheU	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (42/40);  total (42/40)

AAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGC  >  minE/2823494‑2823627                                                                    |                                                                   aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCAc                                                                  >  1:1070739/1‑70 (MQ=255) aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCAc                                                                  >  1:197823/1‑70 (MQ=255) aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCAc                                                                  >  1:68122/1‑70 (MQ=255) aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCAc                                                                  >  1:2914839/1‑70 (MQ=255)  aaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgg                                                                <  1:3180138/71‑1 (MQ=255)    aaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgggg                                                              <  1:2977473/71‑1 (MQ=255)      aaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCAc                                                                  <  1:2227090/65‑1 (MQ=255)      aaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGaa                                                            <  1:1609463/71‑1 (MQ=255)      aaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGaa                                                            <  1:2510090/71‑1 (MQ=255)      aaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGaa                                                            <  1:2529280/71‑1 (MQ=255)        tCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAAc                                                           >  1:1686007/1‑70 (MQ=255)        tCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACg                                                          >  1:1486597/1‑71 (MQ=255)          tACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGgc                                                        >  1:247020/1‑71 (MQ=255)          tACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGgc                                                        >  1:761760/1‑71 (MQ=255)          tACAAAACAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGgc                                                        >  1:411571/1‑71 (MQ=255)            cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACg                                                          <  1:1038682/67‑1 (MQ=255)            cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACg                                                          <  1:2374111/67‑1 (MQ=255)              aaCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTc                                                    >  1:1865893/1‑71 (MQ=255)                 cAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTcc                                                   >  1:1741086/1‑69 (MQ=255)                 cAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCAt                                                 >  1:3117562/1‑71 (MQ=255)                  aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTcc                                                   <  1:389160/68‑1 (MQ=255)                  aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTcc                                                   <  1:518679/68‑1 (MQ=255)                  aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCAtt                                                <  1:2047172/71‑1 (MQ=255)                   aCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTg                                               <  1:2289737/71‑1 (MQ=255)                   aCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTg                                               <  1:297510/71‑1 (MQ=255)                   aCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTg                                               <  1:831/71‑1 (MQ=255)                    cGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTgt                                              <  1:66878/71‑2 (MQ=255)                      cAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCAtt                                                <  1:1495402/67‑1 (MQ=255)                      cAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             >  1:194292/1‑70 (MQ=255)                        acacAATTCATGCCCTGGCAGTATCTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAgg                                            >  1:1294841/1‑69 (MQ=255)                         cacaATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAg                                         >  1:58262/1‑71 (MQ=255)                            aaTTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      >  1:108365/1‑71 (MQ=255)                            aaTTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      <  1:1292833/71‑1 (MQ=255)                            aaTTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      >  1:1074635/1‑71 (MQ=255)                            aaTTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      >  1:1312729/1‑71 (MQ=255)                                  tGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCata                                <  1:219588/71‑1 (MQ=255)                                    cccTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                               >  1:1752741/1‑70 (MQ=255)                                    cccTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                               >  1:1970038/1‑70 (MQ=255)                                     ccTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATT‑AGGAAGTCATTCataa                               >  1:564364/1‑68 (MQ=255)                                      cTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                               >  1:61903/1‑68 (MQ=255)                                       tGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCa                                                  <  1:2238524/48‑1 (MQ=255)                                         gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           <  1:341522/69‑1 (MQ=255)                                         gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           <  1:1368784/69‑1 (MQ=255)                                         gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           <  1:143987/69‑1 (MQ=255)                                          cAGTATGTCACGTTCTCGCGTTTCTGCACGGGGaa                                                            >  1:1502577/1‑35 (MQ=255)                                          cAGTATGTCACGTTCTCGCGTTTCTGCACGGGGaa                                                            >  1:390926/1‑35 (MQ=255)                                          cAGTATGTCACGTTCTCGCGTTTCTGCACGGGGaa                                                            >  1:2649859/1‑35 (MQ=255)                                          cAGTATGTCACGTTCTCGCGTTTCTGCACGGGGaa                                                            >  1:741206/1‑35 (MQ=255)                                           aGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGgcgc                                                      <  1:1727228/40‑1 (MQ=255)                                              aTGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcattca                                  >  1:345137/1‑57 (MQ=255)                                                 tCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTc                                                    <  1:2650176/36‑1 (MQ=255)                                                 tCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             >  1:1915101/1‑43 (MQ=255)                                                 tCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             >  1:3063999/1‑43 (MQ=255)                                                 tCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             >  1:2850734/1‑43 (MQ=255)                                                 tCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             >  1:241762/1‑43 (MQ=255)                                                 tCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             >  1:451345/1‑43 (MQ=255)                                                 tCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCg                 <  1:2484160/71‑1 (MQ=255)                                                 tCACGTTCTAGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             >  1:2400107/1‑43 (MQ=255)                                                   aCGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      <  1:998794/48‑1 (MQ=255)                                                   aCGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaa               <  1:1534777/71‑1 (MQ=255)                                                    cGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaaa              <  1:985402/71‑1 (MQ=255)                                                     gTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcatt                                    >  1:2926202/1‑48 (MQ=255)                                                     gTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcatt                                    >  1:830935/1‑48 (MQ=255)                                                     gTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcatt                                    >  1:1407566/1‑48 (MQ=255)                                                      ttCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAgg                                            <  1:1249023/39‑1 (MQ=255)                                                      ttCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAgg                                            <  1:1621195/39‑1 (MQ=255)                                                      ttCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAgg                                            <  1:1898904/39‑1 (MQ=255)                                                      ttCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGa                <  1:3158657/67‑1 (MQ=255)                                                       tctcGCGTTTCTGCACGGGGAACGGCGCTCCATTGAgg                                            <  1:722396/38‑1 (MQ=255)                                                       tctcGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAgg           >  1:2927743/1‑71 (MQ=255)                                                       tctcGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAgg           >  1:2277146/1‑71 (MQ=255)                                                        ctcGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      >  1:2323747/1‑43 (MQ=255)                                                         tcGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           <  1:1413764/53‑1 (MQ=255)                                                         tcGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           <  1:2692051/53‑1 (MQ=255)                                                          cgcgTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGaa        >  1:1563141/1‑71 (MQ=255)                                                           gcgTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGa         <  1:2331576/69‑1 (MQ=255)                                                           gcgTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAat       >  1:1556047/1‑71 (MQ=255)                                                             gTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAatga     <  1:1837071/71‑1 (MQ=255)                                                                tCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAatgatg   <  1:2029679/70‑1 (MQ=255)                                                                 cTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGc  <  1:1174866/70‑1 (MQ=255)                                                                 cTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGc  >  1:387710/1‑70 (MQ=255)                                                                 cTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCAATCGAAAAGGGAATGATGc  <  1:445153/70‑1 (MQ=255)                                                                    |                                                                   AAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGC  >  minE/2823494‑2823627

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