Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,052,367	T→C	L108S (TTG→TCG)	hemN →	coproporphyrinogen III dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,052,367	0	T	C	100.0%	172.3 / NA	49	L108S (TTG→TCG)	hemN	coproporphyrinogen III dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/21);  total (28/21)

GATCAGTATCTGGACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATTGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGCGCGAAAACT  >  NC_000913/4052291‑4052444                                                                             |                                                                              gATCAGTATCTGGACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCAc                                                                           >  1:816131/1‑81 (MQ=255)       tATCTGGACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGggc                                                                     >  2:816149/1‑81 (MQ=255)       tATCTGGACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGggc                                                                     <  2:816089/81‑1 (MQ=255)        aTCTGGACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGggcg                                                                    >  1:816171/1‑81 (MQ=255)         tCTGGACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGggcg                                                                    <  1:816090/80‑1 (MQ=255)            ggACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGgcggcg                                                                 >  1:816139/1‑80 (MQ=255)            ggACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGgcggcg                                                                 >  1:816138/1‑80 (MQ=255)             gACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGggcgg                                                                   <  2:816091/77‑1 (MQ=255)             gACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGcggcgg                                                                >  2:816178/1‑80 (MQ=255)              aCGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGaa                                                              >  1:816156/1‑81 (MQ=255)                gcgcTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACg                                                            <  2:816092/81‑1 (MQ=255)                   cTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCg                                                         <  2:816093/81‑1 (MQ=255)                    tGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGa                                                        >  1:816163/1‑81 (MQ=255)                    tGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGa                                                        >  1:816183/1‑81 (MQ=255)                      gAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGAc                                                       <  2:816094/80‑1 (MQ=255)                         cAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTAt                                                   >  1:816160/1‑81 (MQ=255)                              aaTCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGaa                                              >  1:816169/1‑81 (MQ=255)                               aTCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGaa                                              >  2:816170/1‑80 (MQ=255)                               aTCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAAt                                             <  1:816095/81‑1 (MQ=255)                                 cGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATa                                            >  1:816182/1‑80 (MQ=255)                                 cGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATa                                            >  2:816189/1‑80 (MQ=255)                                     cATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAgc                                        <  1:816096/80‑1 (MQ=255)                                        cGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCa                                     >  2:816141/1‑80 (MQ=255)                                         gTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCaaa                                   >  2:816179/1‑81 (MQ=255)                                         gTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCaaa                                   >  1:816177/1‑81 (MQ=255)                                         gTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCaaa                                   <  1:816097/81‑1 (MQ=255)                                           gCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAAt                                  >  2:816157/1‑80 (MQ=255)                                            cACCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCa                                <  2:816098/81‑1 (MQ=255)                                             aCCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAg                               >  1:816167/1‑81 (MQ=255)                                             aCCGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAg                               >  1:816159/1‑81 (MQ=255)                                              ccGCTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAgc                              >  1:816161/1‑81 (MQ=255)                                                 cTGTTCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAgccgc                           <  1:816099/81‑1 (MQ=255)                                                    ttCGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAgccgcc                          <  1:816100/79‑1 (MQ=255)                                                      cGCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCtgat                      >  1:816165/1‑81 (MQ=255)                                                       gCCGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCtgatg                     <  1:816101/81‑1 (MQ=255)                                                        ccGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCtgatg                     >  2:816164/1‑80 (MQ=255)                                                        ccGGGCGTCACGTCAGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCtgatg                     >  1:816173/1‑80 (MQ=255)                                                              gtcacgtcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAgctgc              <  2:816102/81‑1 (MQ=255)                                                                cacgtcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGctgctg            >  2:816162/1‑81 (MQ=255)                                                                 acgtcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCtgctgc           <  2:816103/81‑1 (MQ=255)                                                                   gtcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTgcgc         <  2:816105/81‑1 (MQ=255)                                                                   gtcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTgcgc         <  1:816104/81‑1 (MQ=255)                                                                    tcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTgcgc         >  2:816129/1‑80 (MQ=255)                                                                    tcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGcgcg        <  1:816106/81‑1 (MQ=255)                                                                    tcaGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGcgcg        <  1:816107/81‑1 (MQ=255)                                                                      aGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGCGCGaa      <  2:816108/81‑1 (MQ=255)                                                                      aGCCAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGCGCGa       <  1:816109/80‑1 (MQ=255)                                                                         cAATCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGCGCGAAAAc   >  1:816174/1‑81 (MQ=255)                                                                          aaTCGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGCGCGAAAACt  >  1:816180/1‑81 (MQ=255)                                                                             |                                                                              GATCAGTATCTGGACGCGCTGGAGCAAGAAATCGTCCATCGTGCACCGCTGTTCGCCGGGCGTCACGTCAGCCAATTGCACTGGGGCGGCGGAACGCCGACGTATCTGAATAAAGCGCAAATCAGCCGCCTGATGAAGCTGCTGCGCGAAAACT  >  NC_000913/4052291‑4052444

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