Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,770,911	T→A	S73C (AGT→TGT)	hemD ←	uroporphyrinogen‑III synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,770,911	0	T	A	100.0%	91.4 / NA	30	S73C (AGT→TGT)	hemD	uroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (11/19);  total (11/19)

CGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACTTCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGC  >  CP000730/1770777‑1771041                                                                                                                                       |                                                                                                                                   cGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATCCACAATATTGCGCTGTCTTACATCCTAt                                                                                                                              <  1:183613/141‑1 (MQ=255)       ttCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGc                                                                                                                        <  1:357958/141‑1 (MQ=255)            tCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGt                                                                                                                   <  1:40159/141‑1 (MQ=255)                  gCAAAATTATTTTTTGGTTAGTTTGATTAAATGTTTTTAAAAATTCTTCTTGAGAATAGTCGTTTGGCATAAAATCAACTCGTATGCCAAGTGATTCACATTATTGCGTTGTCTTACATCCTATCATAGCAATGTTATCaa                                                                                                             <  2:236595/141‑1 (MQ=255)                     aaaGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACAt                                                                                                          <  1:232583/141‑1 (MQ=255)                      aaGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACAtt                                                                                                         <  2:17339/141‑1 (MQ=255)                             ttttttttttttttATTAAATGATCTTCAAAACCCTTCTTGATATAAGTCGCTTGGCATCAAACCACCTCCAATCCCACGTAATTCACAATATTTTTTTTTCTTACATCCTATCACATCAATGTTATCAACATTAATTCCt                                                                                                  >  2:154485/11‑141 (MQ=255)                             ttttGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCATAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCt                                                                                                  <  2:406397/141‑1 (MQ=255)                               ttGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCttt                                                                                                <  2:171197/141‑1 (MQ=255)                               ttGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCttt                                                                                                <  2:53635/141‑1 (MQ=255)                                            ttAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATaa                                                                                   <  2:128273/141‑1 (MQ=255)                                                       tAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTcac                                                                        <  2:178130/141‑1 (MQ=255)                                                           aaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCa                                                                    >  1:270422/1‑141 (MQ=255)                                                                 tcttGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTa                                                              >  2:186437/1‑141 (MQ=255)                                                                             tCGTTTGGCATAAAATCAACTCGACCGCCAAGTGAGTCACACTATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCACCATTAATTCCTTTTAGACATCTATAAAAAAATTTAACACCACTTTTACATGAAAACATa                                                  >  2:377736/1‑141 (MQ=255)                                                                             tCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATa                                                  >  1:276377/1‑141 (MQ=255)                                                                              cGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATaa                                                 >  1:285510/1‑141 (MQ=255)                                                                                     cATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTc                                          <  2:162295/141‑1 (MQ=255)                                                                                     aaTAAAATCAACTAGAAAGAAAAGTAAAACACAATATTGCGTTGTCTTAAATCCTATAACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGAATAAAAAAAAAAAAAAATAAGCCAGTc                                          <  2:162331/140‑1 (MQ=255)                                                                                      aTAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCa                                         <  1:238994/141‑1 (MQ=255)                                                                                                cTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGAt                               >  1:114429/1‑141 (MQ=255)                                                                                                  cGAATGCCAAGTGATTCACATTATTGCGCTGTCTTAAATTCTAAAACAGAAATGTTATAAATCTTAATTCCTTTTAGATATTTATAAAAGAATTTCAGAGCATTATTAAAAAAAAAAATAAGCCAGTCATAGCGTTGAttt                             <  2:398567/141‑1 (MQ=255)                                                                                                           aGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGa                    >  2:264448/1‑141 (MQ=255)                                                                                                                ttCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATc               <  2:161258/141‑1 (MQ=255)                                                                                                                 tCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCa              >  1:14783/1‑141 (MQ=255)                                                                                                                  cacaATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCaa             >  2:3327/1‑141 (MQ=255)                                                                                                                   acaATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCaaa            >  1:133888/1‑141 (MQ=255)                                                                                                                           gcgcTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTg    <  1:415492/141‑1 (MQ=255)                                                                                                                           gcgcTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTg    <  1:278298/141‑1 (MQ=255)                                                                                                                             gcTGTCTTACATCCTATCACAGAAATGTTATCAACATTAATTCCTGTTAGATATTTATAAAAGAATTTCACAGCAATAAAAGATAAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGc  <  2:65764/141‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   CGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACTTCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGC  >  CP000730/1770777‑1771041

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