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	95.6%	235.4 / 1.5	68	D305E (GAT→GAA)	pdxB	erythronate‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (31/34);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTAT  >  NC_000913/2436841‑2437031                                                                                                 |                                                                                               aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             >  2:1254650/1‑100 (MQ=255) aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             <  1:1190297/100‑1 (MQ=255) aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             >  2:1190290/1‑100 (MQ=255)  tCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGa                                                                                            >  2:476909/1‑100 (MQ=255)    gAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                          <  1:1393210/100‑1 (MQ=255)    gAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                          >  1:862288/1‑100 (MQ=255)      aCTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCgg                                                                                        >  1:269989/1‑100 (MQ=255)       cTCACCCTGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCggg                                                                                       <  2:188764/100‑1 (MQ=255)       cTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCggg                                                                                       <  1:355189/100‑1 (MQ=255)       cTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCggg                                                                                       <  2:1521553/100‑1 (MQ=255)        tCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                      >  1:143252/1‑100 (MQ=255)         cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                     >  1:706866/1‑100 (MQ=255)         cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                     >  2:1467483/1‑100 (MQ=255)         cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                     >  2:405448/1‑100 (MQ=255)            ccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                      <  2:1485088/96‑1 (MQ=255)             cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGTGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGTGCCATGc                                                                                 <  2:1084048/100‑1 (MQ=255)             cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGc                                                                                 >  1:904680/1‑100 (MQ=255)                  tCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGt                                                                            <  2:1089108/100‑1 (MQ=255)                  tCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGt                                                                            <  2:862281/100‑1 (MQ=255)                   cccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTa                                                                           >  2:489551/1‑100 (MQ=255)                    ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                          >  1:8729/1‑100 (MQ=255)                    ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                          >  2:59119/1‑100 (MQ=255)                          aCTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAc                                                                    <  2:1427706/100‑1 (MQ=255)                          aCTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAc                                                                    <  1:1286865/100‑1 (MQ=255)                           cTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAcc                                                                   <  2:269988/100‑1 (MQ=255)                           cTTTACACAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAcc                                                                   <  1:622843/100‑1 (MQ=255)                              tACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCaaa                                                                >  2:1339563/1‑100 (MQ=255)                                cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                              >  1:1353908/1‑100 (MQ=255)                                cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                              >  2:1100419/1‑100 (MQ=255)                                cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                              >  2:1315853/1‑100 (MQ=255)                                    gCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTgg                                                          <  1:1443279/100‑1 (MQ=255)                                    gCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTgg                                                          <  1:414874/100‑1 (MQ=255)                                      ggTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcg                                                        <  1:1449012/100‑1 (MQ=255)                                      ggTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcg                                                        >  1:1078597/1‑100 (MQ=255)                                       gTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcgc                                                       >  2:1050535/1‑100 (MQ=255)                                        tGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgca                                                      <  2:380732/100‑1 (MQ=255)                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                     <  2:934790/100‑1 (MQ=255)                                             cATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggcag                                                 >  2:364220/1‑100 (MQ=255)                                                cgcgGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTaa                                              >  2:1226278/1‑100 (MQ=255)                                                 gcggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAAt                                             <  2:1498503/100‑1 (MQ=255)                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  1:1254658/100‑1 (MQ=255)                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  1:462361/100‑1 (MQ=255)                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  1:1172463/100‑1 (MQ=255)                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAAt                                             <  1:213067/99‑1 (MQ=255)                                                    gcgcACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTa                                          <  1:1483770/100‑1 (MQ=255)                                                          atcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgc                                    >  1:1221435/1‑100 (MQ=255)                                                           tcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgcg                                   >  1:551510/1‑100 (MQ=255)                                                             atACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCa                                 >  1:1445925/1‑100 (MQ=255)                                                                cacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACg                              >  1:638918/1‑100 (MQ=255)                                                                   cAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                           >  1:840811/1‑100 (MQ=255)                                                                   cAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                           <  2:461200/100‑1 (MQ=255)                                                                         cACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCa                     <  2:551507/100‑1 (MQ=255)                                                                         cACCAGCCTTTACAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCa                     >  2:824519/1‑100 (MQ=255)                                                                           ccAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATg                   >  1:63057/1‑100 (MQ=255)                                                                           ccAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATg                   <  2:930403/100‑1 (MQ=255)                                                                             aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                 >  1:338467/1‑100 (MQ=255)                                                                             aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                 <  1:1226285/100‑1 (MQ=255)                                                                              gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                <  2:889684/100‑1 (MQ=255)                                                                              gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                <  2:1321873/100‑1 (MQ=255)                                                                              gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                >  1:65539/1‑100 (MQ=255)                                                                                  tttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATa            >  1:106367/1‑100 (MQ=255)                                                                                  tttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATa            <  1:186193/100‑1 (MQ=255)                                                                                     cAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAAc         >  2:1161001/1‑100 (MQ=255)                                                                                      aGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACt        <  2:1161357/100‑1 (MQ=255)                                                                                         gTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGc     <  2:1080016/100‑1 (MQ=255)                                                                                         gTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGc     <  2:904672/100‑1 (MQ=255)                                                                                          tCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCt    <  1:59119/100‑1 (MQ=255)                                                                                            ggTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCtat  >  2:801957/1‑100 (MQ=255)                                                                                                 |                                                                                               ATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTAT  >  NC_000913/2436841‑2437031

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