Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,769,572	A→G	V295V (GTT→GTC)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,769,572	0	A	G	100.0%	199.1 / NA	55	V295V (GTT→GTC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/26);  total (29/26)

ATCTCATAAAGTAAATTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATAACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTATATCTAAATAACTTAGA  >  CP000730/1769428‑1769707                                                                                                                                                 |                                                                                                                                        atCTCATAAAGAAAATTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACg                                                                                                                                      >  2:267954/1‑148 (MQ=255)   ctcATAAAGTAAATTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATgacgac                                                                                                                                    >  2:412514/1‑148 (MQ=255)     cATAAAGTAAATTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGacac                                                                                                                                  >  2:323781/1‑148 (MQ=255)     cATAAAGTAAATTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGacac                                                                                                                                  >  2:55065/1‑148 (MQ=255)     cATAAAGTAAATTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGacac                                                                                                                                  >  2:401396/1‑148 (MQ=255)             aaaTTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCa                                                                                                                          <  2:177776/148‑1 (MQ=255)             aaaTTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCa                                                                                                                          <  2:203604/148‑1 (MQ=255)             aaaTTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCa                                                                                                                          >  2:171438/1‑148 (MQ=255)             aaaTTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCACTGAAATTATTTGTTCCATGGCGACACGTTCTTCa                                                                                                                          <  2:415879/148‑1 (MQ=255)                ttCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCt                                                                                                                       >  1:79399/1‑148 (MQ=255)                ttCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCt                                                                                                                       >  1:246172/1‑148 (MQ=255)                  cctcctACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTAt                                                                                                                     >  2:128510/1‑148 (MQ=255)                  cctcctACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTAt                                                                                                                     >  2:354626/1‑148 (MQ=255)                                 ggAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGc                                                                                                      >  2:345078/1‑148 (MQ=255)                                 ggAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGc                                                                                                      >  2:218754/1‑148 (MQ=255)                                            aaTATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAg                                                                                           <  1:340730/148‑1 (MQ=255)                                            aaTATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAg                                                                                           <  1:93948/148‑1 (MQ=255)                                            aaTATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCACGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAg                                                                                           <  2:358448/148‑1 (MQ=255)                                              tataAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTc                                                                                         >  2:208734/1‑148 (MQ=255)                                              tataAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTCTGTGCCGCTGCTTTAGTc                                                                                         <  1:116010/148‑1 (MQ=255)                                               ataAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCa                                                                                        <  1:407109/148‑1 (MQ=255)                                               ataAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCa                                                                                        <  1:304676/148‑1 (MQ=255)                                                taAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCAt                                                                                       <  2:97393/148‑1 (MQ=255)                                                taAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCAt                                                                                       <  2:149625/148‑1 (MQ=255)                                                        tatttatCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATAtt                                                                               >  1:407378/1‑148 (MQ=255)                                                        tatttatCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATAtt                                                                               >  1:361871/1‑148 (MQ=255)                                                        tatttatCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACCATAtt                                                                               >  2:339218/1‑148 (MQ=255)                                                           ttatCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATtctc                                                                            >  1:148320/1‑148 (MQ=255)                                                              tCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCTAc                                                                         >  2:52722/1‑148 (MQ=255)                                                              tCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCAc                                                                         <  1:7791/148‑1 (MQ=255)                                                              tCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCAc                                                                         <  1:313014/148‑1 (MQ=255)                                                               cTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACt                                                                        >  2:388759/1‑148 (MQ=255)                                                                tAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTc                                                                       <  2:283249/148‑1 (MQ=255)                                                                 aaaTAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTca                                                                      >  1:191838/1‑148 (MQ=255)                                                                       cGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATtat                                                                <  2:32460/148‑1 (MQ=255)                                                                         aCAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCTCGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTatat                                                              >  1:220433/1‑148 (MQ=255)                                                                                 ccTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCaacaac                                                      >  2:278932/1‑148 (MQ=255)                                                                                 ccTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCaacaac                                                      >  2:190768/1‑148 (MQ=255)                                                                                     tGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGa                                                  >  1:1268/1‑148 (MQ=255)                                                                                                  ataatCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAt                                     <  1:128510/148‑1 (MQ=255)                                                                                                  ataatCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAt                                     <  1:354626/148‑1 (MQ=255)                                                                                                   taatCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAtt                                    <  1:55065/148‑1 (MQ=255)                                                                                                   taatCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAtt                                    <  1:401396/148‑1 (MQ=255)                                                                                                   taatCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAtt                                    <  1:323781/148‑1 (MQ=255)                                                                                                   taatCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAtt                                    <  1:267954/148‑1 (MQ=255)                                                                                                      tCATATCAGTACCTGCACGTTTCATTGAAACCATTTGTTCTATGTCGTCACGTTCTTCATCTATCTATCCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTTTATGCAACTACTGGAACATTCGTATGAttttt                                 <  2:268365/148‑1 (MQ=255)                                                                                                      tCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAttttt                                 <  2:315274/148‑1 (MQ=255)                                                                                                      tCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGAttttt                                 <  2:36074/148‑1 (MQ=255)                                                                                                                       cGTTTCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTatat                <  1:241868/148‑1 (MQ=255)                                                                                                                          ttCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTATATCTa             >  1:85377/1‑148 (MQ=255)                                                                                                                          ttCATTGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTATATCTa             >  1:319258/1‑148 (MQ=255)                                                                                                                               tGAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTATATCTAAATaa        <  1:412514/148‑1 (MQ=255)                                                                                                                                gAAACCATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTAGATCTAAATAAc       >  2:304561/1‑148 (MQ=255)                                                                                                                                     cATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGCCCATTTGTAACCGTTGGATCGTCTGGAAGTATATCTAAATGACtgatg  >  1:185492/1‑144 (MQ=255)                                                                                                                                     cATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTATATCTAAATAACTTaga  >  1:392977/1‑148 (MQ=255)                                                                                                                                     cATTTGTTCCATGACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTATATCTAAATAACTTaga  >  1:228035/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                        ATCTCATAAAGTAAATTCCTCCTACAGTTTATGGAAAATCATAAATATAAAACCTTATTTATCTAAATAGCGACAAATGTCCTTTGCAAAATACGTAATAATCATATCAGCACCTGCACGTTTCATTGAAACCATTTGTTCCATAACGACACGTTCTTCATCTATCCAACCATTTTGTGCCGCTGCTTTAGTCATACTATATTCTCCACTCACATTATATGCAACAACTGGAACATTCGTATGATTTTTAACATCTCGAACTATATCTAAATAACTTAGA  >  CP000730/1769428‑1769707

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