Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,990,203	T→C	P188P (CCA→CCG)	hemC ←	hydroxymethylbilane synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,990,203	0	T	C	100.0%	113.6 / NA	33	P188P (CCA→CCG)	hemC	hydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/14);  total (19/14)

GGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAAT  >  NC_000913/3990068‑3990337                                                                                                                                        |                                                                                                                                       ggCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCa                                                                                                                                    >  1:257533/1‑140 (MQ=255)    gcgTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAAcg                                                                                                                                 >  1:280617/1‑140 (MQ=255)      gTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAAcgcg                                                                                                                               <  2:246422/140‑1 (MQ=255)       ttCTGTGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACgcgg                                                                                                                              >  1:362286/1‑140 (MQ=255)             ggTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGaa                                                                                                                        <  2:257533/140‑1 (MQ=255)                  cgcgCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACg                                                                                                                    <  1:148423/139‑1 (MQ=255)                  cgcgCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACg                                                                                                                    >  2:148423/1‑139 (MQ=255)                    cgcAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGa                                                                                                                 >  1:151800/1‑140 (MQ=255)                       agtgcagtTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTc                                                                                                              <  2:36051/140‑1 (MQ=255)                        gtgcagtTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATCCGTGACTcc                                                                                                             <  2:362286/140‑1 (MQ=255)                                 cGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                                                                    >  2:452081/1‑140 (MQ=255)                                   tggtgATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCTGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                                                                                  >  1:352908/1‑140 (MQ=255)                                   tggtgATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                                                                                  >  2:40623/1‑140 (MQ=255)                                          tCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                                                           <  1:155791/140‑1 (MQ=255)                                                cgGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                                                                                     <  2:234023/140‑1 (MQ=255)                                                            cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                                                                         >  2:232588/1‑140 (MQ=255)                                                                   gCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCAt                                                                  <  1:17542/140‑1 (MQ=255)                                                                    cGTGAATCACCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                                                                 >  1:350077/1‑140 (MQ=255)                                                                                aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGtt                                                     <  2:98057/140‑1 (MQ=255)                                                                                 ggcggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTa                                                    >  2:150198/1‑140 (MQ=255)                                                                                    ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATcc                                                 >  2:124968/1‑140 (MQ=255)                                                                                            aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCt                                         <  1:124968/140‑1 (MQ=255)                                                                                             tACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTc                                        <  1:193819/140‑1 (MQ=255)                                                                                                  aCCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGATACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCg                                   >  2:376150/1‑140 (MQ=255)                                                                                                     gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGt                                >  2:457193/1‑140 (MQ=255)                                                                                                                tACCGCCTGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGc                     <  1:325784/140‑1 (MQ=255)                                                                                                                  ccgccgGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCcg                   >  2:276886/1‑140 (MQ=255)                                                                                                                     ccgGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCgcgc                >  2:151086/1‑140 (MQ=255)                                                                                                                     ccgGAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCgcgc                >  2:415518/1‑140 (MQ=255)                                                                                                                       ggAAGAGAAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAg              >  2:410559/1‑140 (MQ=255)                                                                                                                           gagaAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAg          <  2:151800/140‑1 (MQ=255)                                                                                                                            agaAATCTCGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGc         <  1:151086/140‑1 (MQ=255)                                                                                                                                   tcGGGCGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAAt  >  2:146402/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       GGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAAT  >  NC_000913/3990068‑3990337

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑