Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,264,919	C→T	D402D (GAC→GAT)	hemA →	glutamyl tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,264,919	0	C	T	100.0%	204.0 / NA	64	D402D (GAC→GAT)	hemA	glutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/35);  total (29/35)

AGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGC  >  NC_000913/1264786‑1265048                                                                                                                                      |                                                                                                                                  aGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGa                                                                                                                              >  2:1623791/1‑139 (MQ=255)        cGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATaacgaa                                                                                                                       <  1:1677884/139‑1 (MQ=255)          ccAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATaacgaacg                                                                                                                     <  2:885411/139‑1 (MQ=255)           cAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGc                                                                                                                    >  1:1012852/1‑139 (MQ=255)           cAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGc                                                                                                                    >  1:1524859/1‑139 (MQ=255)                           cTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGACCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATaacgaacg                                                                                                                     >  2:1376935/1‑122 (MQ=255)                           cTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGACCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATaacgaacg                                                                                                                     <  1:1376935/122‑1 (MQ=255)                                       ggcgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGc                                                                                        <  2:1524859/139‑1 (MQ=255)                                       ggcgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGc                                                                                        <  1:1623791/139‑1 (MQ=255)                                            cgcgCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAg                                                                                   <  1:1115027/139‑1 (MQ=255)                                                    ccATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTa                                                                           <  2:1012852/139‑1 (MQ=255)                                                               gATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCt                                                                >  1:1121021/1‑139 (MQ=255)                                                                      cATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTCACAGTACATCATTTTCTTTTTTTa                                                         <  2:1121021/139‑1 (MQ=255)                                                                        tGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                       <  1:1460572/139‑1 (MQ=255)                                                                        tGGAAACGGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACTAACGCCTGCATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                       <  1:1042387/139‑1 (MQ=255)                                                                               tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGAc                                                                                                  <  2:786020/89‑1 (MQ=255)                                                                               tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGAc                                                                                                  >  1:786020/1‑89 (MQ=255)                                                                               tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCt                                                                                             <  2:95624/94‑1 (MQ=255)                                                                               tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCt                                                                                             >  1:95624/1‑94 (MQ=255)                                                                                 aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              <  1:766608‑M1/139‑5 (MQ=255)                                                                                 aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              <  2:121321‑M1/139‑5 (MQ=255)                                                                                 aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              <  2:1067743‑M1/139‑5 (MQ=255)                                                                                     aCCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggagga                                          >  1:975785‑M1/1‑131 (MQ=255)                                                                                      ccGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaa                                         >  1:899856‑M1/1‑130 (MQ=255)                                                                                          ttGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataat                                     <  1:713544‑M1/139‑14 (MQ=255)                                                                                           tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                      <  2:428891/121‑1 (MQ=255)                                                                                           tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                      >  1:428891/1‑121 (MQ=255)                                                                                           tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaata                                       >  1:487306‑M1/1‑125 (MQ=255)                                                                                           tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaata                                       <  2:487306‑M1/136‑12 (MQ=255)                                                                                           tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              <  1:1746432‑M1/129‑5 (MQ=255)                                                                                           tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              >  2:1746432‑M1/1‑125 (MQ=255)                                                                                           tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTCgaggaggaataatg                                    >  2:1129965‑M1/1‑125 (MQ=255)                                                                                            gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                   <  2:848934‑M1/139‑16 (MQ=255)                                                                                            gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                   <  1:1129965‑M1/139‑16 (MQ=255)                                                                                            gATCCACGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                   <  1:1707206‑M1/139‑16 (MQ=255)                                                                                                 aTGCGCCAACGAAATCATTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTaac                                                                               >  1:1050171/1‑88 (MQ=255)                                                                                                 aTGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTaac                                                                               <  2:1050171/90‑3 (MQ=255)                                                                                                  tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGt                                                                                  >  2:1491260/1‑86 (MQ=255)                                                                                                  tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGt                                                                                  <  1:1491260/86‑1 (MQ=255)                                                                                                        aaCGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCt                                                                                             >  1:201827/1‑69 (MQ=255)                                                                                                        aaCGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCt                                                                                             <  2:201827/69‑1 (MQ=255)                                                                                                        aaCGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCt                                                                                       <  1:101957/75‑1 (MQ=255)                                                                                                        aaCGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCt                                                                                       >  2:101957/1‑75 (MQ=255)                                                                                                            aaaTCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTgg                                                                                     >  1:148152/1‑73 (MQ=255)                                                                                                            aaaTCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTgg                                                                                     <  2:148152/73‑1 (MQ=255)                                                                                                              aTCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTgg                                                                                     <  2:340352/71‑1 (MQ=255)                                                                                                              aTCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTgg                                                                                     >  1:340352/1‑71 (MQ=255)                                                                                                                    tCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCt                                                                                             >  1:1161156/1‑57 (MQ=255)                                                                                                                    tCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCt                                                                                             <  2:1161156/57‑1 (MQ=255)                                                                                                                     cAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagga                                             >  2:1386682‑M1/1‑99 (MQ=255)                                                                                                                     cAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagga                                             <  1:1386682‑M1/104‑6 (MQ=255)                                                                                                                           gccgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGa                                                                                                   >  1:1250322/1‑44 (MQ=255)                                                                                                                           gccgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGa                                                                                                   <  2:1250322/44‑1 (MQ=255)                                                                                                                           gccgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataag    >  1:1385201‑M1/1‑93 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGt                                                                                  >  1:205490/1‑59 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGt                                                                                  <  2:205490/59‑1 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc  >  1:794716‑M1/1‑91 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc  >  1:1352734‑M1/1‑91 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              >  1:231137‑M1/1‑91 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              >  1:1749068‑M1/1‑91 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              <  2:1749068‑M1/95‑5 (MQ=255)                                                                                                                             cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                              <  2:231137‑M1/95‑5 (MQ=255)                                                                                                                               cccGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGAc                                                                                                  <  2:1477509/41‑1 (MQ=255)                                                                                                                               cccGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGAc                                                                                                  >  1:1477509/1‑41 (MQ=255)                                                                                                                                      |                                                                                                                                  AGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGC  >  NC_000913/1264786‑1265048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑