Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,436,937	A→T	D305E (GAT→GAA)	pdxB ←	erythronate‑4‑phosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,436,937	0	A	T	100.0%	183.9 / NA	54	D305E (GAT→GAA)	pdxB	erythronate‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/33);  total (21/33)

TATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACT  >  NC_000913/2436840‑2437025                                                                                                  |                                                                                         tATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATTCACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTc                                                                                        <  1:920476/100‑1 (MQ=255) tATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTc                                                                                        <  2:405140/100‑1 (MQ=255) tATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTc                                                                                        <  1:195083/100‑1 (MQ=255)  aTCGAACTGAGCCGGTATCCCGGCGACTTTAGGCAGCGGTGCGTCAGCGCGGCGCACATTATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                       <  2:409281/100‑1 (MQ=255)  aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                       >  1:495549/1‑100 (MQ=255)  aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                       <  1:385209/100‑1 (MQ=255)  aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                       <  2:750356/100‑1 (MQ=255)  aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                       <  1:288354/100‑1 (MQ=255)  aTCGAACTCACCCGGTATACCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                       >  2:209337/1‑100 (MQ=255)     gAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                    <  1:821621/100‑1 (MQ=255)     gAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                    >  2:717642/1‑100 (MQ=255)         tCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                <  1:115482/100‑1 (MQ=255)          cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                               >  1:650911/1‑100 (MQ=255)          cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                               >  1:630441/1‑100 (MQ=255)            cccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCAt                                                                             <  1:266976/100‑1 (MQ=255)              cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGc                                                                           >  1:694533/1‑100 (MQ=255)              cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGc                                                                           <  1:48914/100‑1 (MQ=255)                  aTCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCg                                                                       >  1:1230761/1‑100 (MQ=255)                   tCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCg                                                                       <  2:1204110/99‑1 (MQ=255)                     ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                    >  1:1190926/1‑100 (MQ=255)                     ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                    >  2:130170/1‑100 (MQ=255)                       ggCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATg                                                                  <  2:113862/100‑1 (MQ=255)                       ggCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCATaa                                                                    <  2:1111944/98‑1 (MQ=255)                        gCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGc                                                                 <  1:682897/100‑1 (MQ=255)                        gCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGc                                                                 <  1:1273662/100‑1 (MQ=255)                           aCTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAc                                                              <  1:712474/100‑1 (MQ=255)                                 cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                        >  1:723677/1‑100 (MQ=255)                                    aGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTg                                                     <  2:808210/100‑1 (MQ=255)                                      cGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgc                                                   <  1:781001/100‑1 (MQ=255)                                          gCGTCATCGCGGCGCACATCACACACCAAATGCACCACCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                               <  1:126306/100‑1 (MQ=255)                                           cGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcagg                                              >  1:370741/1‑100 (MQ=255)                                            gTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCGTGCAGCGTAATGCGACCAAACTCTGGCgcaggc                                             <  2:782473/100‑1 (MQ=255)                                            gTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggc                                             <  1:390327/100‑1 (MQ=255)                                            gTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggc                                             <  2:194480/100‑1 (MQ=255)                                            gTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggc                                             <  2:495545/100‑1 (MQ=255)                                             tCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggca                                            >  1:345657/1‑100 (MQ=255)                                               aTCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGt                                          >  1:398900/1‑100 (MQ=255)                                                   cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                      <  2:685262/100‑1 (MQ=255)                                                   cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                      >  2:410515/1‑100 (MQ=255)                                                     gcgcACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTa                                    <  1:410519/100‑1 (MQ=255)                                                            tcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgcg                             <  2:817367/100‑1 (MQ=255)                                                            tcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgcg                             <  2:666559/100‑1 (MQ=255)                                                             catACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGATCGGGCCATGCAGCGTAATGCGACCAAACTCTGTCGCAGTCAGTAATGTATCCAgcgc                            >  2:239966/1‑100 (MQ=255)                                                             catACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgcgc                            >  2:171285/1‑100 (MQ=255)                                                              atACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCa                           <  1:550490/100‑1 (MQ=255)                                                                 cacCAAATGCACCAGCCTTTTCAGCGTCGTTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACg                        <  2:783961/100‑1 (MQ=255)                                                                     aaaTGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCt                    >  2:1304321/1‑100 (MQ=255)                                                                          cACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCa               >  2:226330/1‑100 (MQ=255)                                                                            ccAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATg             >  1:1125490/1‑100 (MQ=255)                                                                             cAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGc            >  1:268256/1‑100 (MQ=255)                                                                             cAGCCTTTTCAGCGTCCGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTATGGCGCAGGCAGTGATGTACCCAGCGCAACGTGCtcgtca               >  2:815825/1‑92 (MQ=255)                                                                               gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc          <  1:302764/100‑1 (MQ=255)                                                                                  ttttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCaan       <  1:338574/100‑2 (MQ=255)                                                                                       aGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACt  <  2:474872/100‑1 (MQ=255)                                                                                                  |                                                                                         TATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACT  >  NC_000913/2436840‑2437025

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: ‑