Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,769,572	A→G	V295V (GTT→GTC)	hemB ←	Delta‑aminolevulinic acid dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,769,572	0	A	G	100.0%	230.6 / NA	62	V295V (GTT→GTC)	hemB	Delta‑aminolevulinic acid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/29);  total (33/29)

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