Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,051,170	+C	intergenic (+172/+279)	hemN → / ← glnG	O2‑independent coproporphyrinogen III oxidase/GlnG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,051,170	1	.	C	100.0%	181.7 / NA	78	intergenic (+172/+279)	hemN/glnG	O2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (51/27);  total (51/27)

GAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC  >  NC_000913/4051095‑4051239                                                                             |                                                                      gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                       <  2:1655209/76‑2 (MQ=255) gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                       <  1:2178364/76‑2 (MQ=255)   aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgctg                                                                     <  2:2082580/76‑4 (MQ=21)       attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                                 <  2:2723264/76‑1 (MQ=21)       attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                                 <  2:1749076/76‑1 (MQ=18)         tggttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                               <  2:1253639/73‑1 (MQ=18)          gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                              >  1:2275095/1‑76 (MQ=21)           cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                             <  1:1701747/76‑1 (MQ=21)            ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGGCt                                                            >  1:2443486/1‑76 (MQ=18)            ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCCGCGCTCGCCGCGTCt                                                            >  2:410167/1‑76 (MQ=17)             tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                           >  1:3966757/1‑76 (MQ=18)             tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                           >  1:3874497/1‑76 (MQ=18)              ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                          >  1:3672564/1‑76 (MQ=18)               tGTAGGTCGGCTAATGCGTTTACTCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                         <  1:1146967/76‑1 (MQ=17)               tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                         <  1:1169290/76‑1 (MQ=18)                 tAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                       <  2:2961998/76‑1 (MQ=18)                  aGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                      >  1:1998249/1‑76 (MQ=18)                  aGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                      >  1:3455993/1‑76 (MQ=18)                  aGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                      >  2:2907271/1‑76 (MQ=17)                   ggCGGGATGGGGCGTTGACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                                                     <  2:257473/76‑1 (MQ=14)                   ggCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                                                     <  1:4817138/76‑1 (MQ=17)                   ggCCGAATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGTGGCGCTGGTCGCGTCTTATCAgg                                                     <  1:4821296/76‑1 (MQ=17)                    gCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                    <  2:3917072/76‑1 (MQ=17)                      cGGATAAGGCGTTTACGCCGCATCA‑GTCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                  <  2:1888872/76‑1 (MQ=17)                        gATAGGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                <  2:426568/76‑1 (MQ=14)                           aaGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                                             <  2:3029404/76‑1 (MQ=17)                            aGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                            >  2:2724954/1‑76 (MQ=18)                                gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                        >  1:2860202/1‑76 (MQ=18)                                gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                        >  2:3002451/1‑76 (MQ=18)                                gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                        >  2:2968455/1‑76 (MQ=18)                                gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGCCTTCTCCGCGCCACGACtttgc                                        >  2:1746167/1‑74 (MQ=9)                                 tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGCCTTtctc                                       >  2:444669/1‑76 (MQ=17)                                 tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                       <  2:2496165/76‑1 (MQ=18)                                 tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                       <  2:3650588/76‑1 (MQ=18)                                 tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCCACGACTTtctc                                       >  1:4753504/1‑76 (MQ=18)                                 tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGAGCTGGTCGCGTCTTATCAGGCCTACGACCTtctc                                       >  2:566171/1‑76 (MQ=18)                                 tttACGCCGCATCAGG‑CAGCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                       <  2:2760819/76‑1 (MQ=18)                                  ttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCCGGGCTACGACTTTCTcc                                      >  2:878129/1‑76 (MQ=18)                                  ttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                      >  1:1396638/1‑76 (MQ=21)                                   tACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                     >  2:1792341/1‑76 (MQ=21)                                   tACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                     >  2:1757836/1‑76 (MQ=21)                                    aCGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                                    >  2:2821285/1‑76 (MQ=21)                                    aCGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                                    <  2:1594107/76‑1 (MQ=18)                                     cgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                   <  1:673766/76‑1 (MQ=21)                                     cgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                   >  2:3980126/1‑76 (MQ=21)                                      gccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                  >  1:64683/1‑76 (MQ=21)                                      gccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                  >  1:4782650/1‑76 (MQ=21)                                      gccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                  >  2:76887/1‑76 (MQ=21)                                        cgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                                <  1:3355365/76‑1 (MQ=21)                                         gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                               >  1:647711/1‑76 (MQ=25)                                          cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCCTcccg                              >  2:284885/1‑74 (MQ=21)                                          cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                              >  2:277183/1‑76 (MQ=25)                                           aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             >  2:1436396/1‑76 (MQ=25)                                           aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             >  2:3289364/1‑76 (MQ=25)                                           aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             >  1:1465514/1‑76 (MQ=25)                                           aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             >  2:3559681/1‑76 (MQ=25)                                                 gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCCTaaa                        >  2:247354/2‑76 (MQ=25)                                                 gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa                        >  2:2394769/2‑76 (MQ=255)                                                  cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                       <  1:3008734/76‑1 (MQ=255)                                                   aTCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                      >  2:3478174/1‑76 (MQ=255)                                                    tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATCaacac                     >  2:4691471/1‑75 (MQ=25)                                                    tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                     >  1:1967827/1‑76 (MQ=255)                                                    tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                     >  1:765809/1‑76 (MQ=255)                                                    tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                     >  2:2644707/1‑76 (MQ=255)                                                     cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCCTCCGGCATAcacaac                    >  1:3038878/1‑76 (MQ=25)                                                     cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGGCTTACCACGCCTACGACTTTCTCCGCATCCGGCATAaacaac                    >  2:824130/1‑76 (MQ=37)                                                      tGAGCCTGATGCCTGCTTCTGCGCTGTTCTCTTCTTATCCTGTCTACGACTTTCTCCGCATCCGGCATCaacacag                   >  1:4330589/1‑73 (MQ=255)                                                      tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                   >  2:2134315/1‑76 (MQ=255)                                                      tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                   >  1:4314360/1‑76 (MQ=255)                                                       gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa                  >  2:2343043/1‑76 (MQ=255)                                                        aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAg                 >  2:4480889/1‑76 (MQ=255)                                                          cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACCACGAGcc               >  2:74493/1‑75 (MQ=25)                                                          cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCt               >  1:1981523/1‑76 (MQ=255)                                                             gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa            <  2:2292938/76‑1 (MQ=255)                                                                gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt         <  1:3010272/76‑1 (MQ=255)                                                                   tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaaa      >  1:2449714/1‑76 (MQ=255)                                                                     gtGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa    <  2:919520/75‑1 (MQ=255)                                                                     aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATACACAACGAGCTTCAGCTAAACa    >  2:612348/1‑76 (MQ=38)                                                                     aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa    <  1:3519892/76‑1 (MQ=255)                                                                     aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa    <  1:26047/76‑1 (MQ=255)                                                                      tGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAg   <  1:137341/76‑1 (MQ=255)                                                                       gcggcgcTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAANAACGAGCTTCAGCTAAACAGc  <  1:2763564/76‑1 (MQ=255)                                                                       gcggcgcTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc  >  1:4204304/1‑76 (MQ=255)                                                                             |                                                                      GAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC  >  NC_000913/4051095‑4051239

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