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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,770,9110TA100.0% 71.6 / ‑6.5 27S73C (AGT→TGT) hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (12/13);  minor base C (1/0);  total (14/13)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCTCGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACTTCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCA  >  CP000730/1770773‑1771029
                                                                                                                                          |                                                                                                                      
gTCTCGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATc                                                                                                                      >  1:84782/1‑141 (MQ=255)
   tcGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGAGTAAATGATTTTAAAACTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACCTCCTa                                                                                                                   >  1:103182/1‑141 (MQ=255)
   tcGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTAAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTa                                                                                                                   >  1:226597/1‑141 (MQ=255)
         aTTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAg                                                                                                             >  1:175005/1‑141 (MQ=255)
            cACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCaa                                                                                                          <  2:233554/141‑1 (MQ=255)
                cGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGtt                                                                                                      <  1:152560/141‑1 (MQ=255)
                cGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGtt                                                                                                      <  2:188841/141‑1 (MQ=255)
                     gCAAAAGTATTTTTAGGCTAGATTGAGTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCaa                                                                                                 >  1:69054/1‑141 (MQ=255)
                       aaaaGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCATCTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACa                                                                                               <  2:17885/141‑1 (MQ=255)
                          aGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTa                                                                                            <  2:128392/141‑1 (MQ=255)
                                 tttGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAACGTCGCTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCtt                                                                                     <  1:54162/141‑1 (MQ=255)
                                 tttGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCtt                                                                                     <  1:52037/141‑1 (MQ=255)
                                         gTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGTCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGatat                                                                             <  1:194076/141‑1 (MQ=255)
                                                          tAAAAATCCTTCTTCAGAAAAATCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTGACATCCTATCACTGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTcac                                                            >  2:196341/1‑141 (MQ=255)
                                                            aacagccTTCTTGATAAAAGTCTTTTGGTATAAATTCAACTCATATGCTAAGTGATTCACAATATTTCTCTTTTTTACATCCTCTCACAGCAATGTGATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAg                                                          <  2:25068/136‑1 (MQ=255)
                                                            aaaaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCACTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTCGATATTTATAAAAGAAATTCACAg                                                          >  2:30430/1‑141 (MQ=255)
                                                                 ccTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCAttt                                                     >  1:144550/1‑141 (MQ=255)
                                                                 ccTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCAttt                                                     >  1:180988/1‑141 (MQ=255)
                                                                                 cGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATaa                                     <  2:175005/141‑1 (MQ=255)
                                                                                         aTAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCa                             >  1:276093/1‑141 (MQ=255)
                                                                                           aaaaTCAAAAAGAATGCAAAAAGAATCACAATATTGCGCTGTATTACAACCTAACACAACAATGTTATCAACACTAAATCCTATTAGATATTTATAAAAGAATTTCACAGCAAAAAAAAAAAAAAAAATAAGCCAGTCATa                           <  2:226597/141‑1 (MQ=255)
                                                                                              aTCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCg                        <  2:155728/141‑1 (MQ=255)
                                                                                                   cTCGAATGCCAAGTAATTCACAATATTGCTCTGTCTTGCATCCTATCACAGCAATGTTATCAACATTAATTCCTCTTAGATAATTATAAAAGAATTTCACAGTATTTTTAGATGAAAGAATAAGACAGTCATAGCGTTGAt                   >  2:213938/1‑141 (MQ=255)
                                                                                                       aaTGCCAAGTGATTCACAATATTTCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTTGTTATTTATAAAAGAATTTCACAGCATTTTTAGATGTAAAAATAAGCCAGTCATAGCGTTGATTTaa               >  1:261084/1‑141 (MQ=255)
                                                                                                              aGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGa        >  1:175459/1‑141 (MQ=255)
                                                                                                                    tcacgatgctGCGCTGTCTTACAGCTTATCACCGCACTGCTATCAGCATTAAATCGTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCa  <  1:273966/132‑1 (MQ=255)
                                                                                                                    tCACAATATTGCGCTGTCTTACGTCATCGTGCAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCCTATTTAGATGAAAAAATAAGCCAGTCATAGCGTGGATTTAACACATGAatatag  >  2:276256/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                      
GTCTCGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACTTCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCA  >  CP000730/1770773‑1771029

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