Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,073,474 T→C D73G (GAC→GGC)  epd ← D‑erythrose 4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,073,4740TC82.3% 146.9 / 18.0 62D73G (GAC→GGC) epdD‑erythrose 4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base C (28/23);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTAATTTCCGC  >  NC_000913/3073345‑3073610
                                                                                                                                 |                                                                                                                                        
ccGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACcaa                                                                                                                                 >  2:395566/1‑139 (MQ=255)
                    gCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTcc                                                                                                             >  2:1680807/1‑139 (MQ=255)
                        tGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTc                                                                                                         >  2:1536930/1‑139 (MQ=255)
                          cTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGt                                                                                                       <  1:1002566/139‑1 (MQ=255)
                             gcgGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                <  1:751377/127‑1 (MQ=255)
                             gcgGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                >  2:751377/1‑127 (MQ=255)
                                       tataCGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaa                                                                                          <  1:1680807/139‑1 (MQ=255)
                                           cGCCGGTGCAGTCGAGGACTACATCAACTCCAAGCTCACGCCAGGGGAGGGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCTCGTTCCTGTCGTACTTCCCATGCAAAACg                                                                                      >  2:1119564/1‑139 (MQ=255)
                                               ggTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCa                                                                                  >  2:1617405/1‑139 (MQ=255)
                                                 tGCAGTCGAGGAGTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATg                                                                                >  2:454693/1‑139 (MQ=255)
                                                 tGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTTCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATg                                                                                >  2:278235/1‑139 (MQ=255)
                                                 tGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATg                                                                                >  1:1469550/1‑139 (MQ=255)
                                                 tGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAAGGGAGAGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTCGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATg                                                                                >  1:1159269/1‑139 (MQ=255)
                                                  gCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATgg                                                                               <  1:1536930/139‑1 (MQ=255)
                                                             cTACATCAACGCCAAGTTCACGCCAGGGGAGCGAATGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCata                                                                    >  2:1694652/1‑139 (MQ=255)
                                                                cATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATAttt                                                                 >  1:442195/1‑139 (MQ=255)
                                                                   cAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTcaa                                                              <  1:395566/139‑1 (MQ=255)
                                                                     aCGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaaaa                                                                                        <  1:936208/111‑1 (MQ=255)
                                                                     aCGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaaaa                                                                                        >  2:936208/1‑111 (MQ=255)
                                                                     aCGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTcaaca                                                            >  2:808855/1‑139 (MQ=255)
                                                                        ccAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCatat                                                                   <  1:1422406/129‑1 (MQ=255)
                                                                        ccAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCatat                                                                   >  2:1422406/1‑129 (MQ=255)
                                                                        ccAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGACCATGGCTGGTGTCATATTTCAACAAAt                                                         >  1:990257/1‑139 (MQ=255)
                                                                          aaGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATgc                                                       <  1:808855/139‑1 (MQ=255)
                                                                           aGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATgcg                                                      <  2:442195/139‑1 (MQ=255)
                                                                           aGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGc                                                                              <  2:1569089/115‑1 (MQ=255)
                                                                           aGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGc                                                                              >  1:1569089/1‑115 (MQ=255)
                                                                           aGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATgcg                                                      >  1:853631/1‑139 (MQ=255)
                                                                             ttCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGcc                                                    <  1:278235/139‑1 (MQ=255)
                                                                                    cAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                <  2:440635/72‑1 (MQ=255)
                                                                                    cAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                >  1:440635/1‑72 (MQ=255)
                                                                                       gggAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                              >  1:1199986/1‑71 (MQ=255)
                                                                                       gggAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                              <  2:1199986/71‑1 (MQ=255)
                                                                                       gggAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCa                                                                                               <  2:300373/86‑1 (MQ=255)
                                                                                       gggAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCa                                                                                               >  1:300373/1‑86 (MQ=255)
                                                                                       gggAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCgcag                                          >  2:1300728/1‑139 (MQ=255)
                                                                                            cGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                              <  2:837591/66‑1 (MQ=255)
                                                                                            cGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                              >  1:837591/1‑66 (MQ=255)
                                                                                            cGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATcg                                     >  1:1379265/1‑138 (MQ=255)
                                                                                              aTTGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGc                                   <  1:1461988/139‑1 (MQ=255)
                                                                                                tGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                <  1:1559754/60‑1 (MQ=255)
                                                                                                tGCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                >  2:1559754/1‑60 (MQ=255)
                                                                                                 gCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                              <  1:785750/61‑1 (MQ=255)
                                                                                                 gCAGTGAACGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                              >  2:785750/1‑61 (MQ=255)
                                                                                                   aGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCgcagc                                         <  1:343079/128‑1 (MQ=255)
                                                                                                   aGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCgcagc                                         >  2:343079/1‑128 (MQ=255)
                                                                                                         cGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaaa                                                                                         >  1:1270669/1‑74 (MQ=255)
                                                                                                         cGTTCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaaa                                                                                         <  2:1270669/74‑1 (MQ=255)
                                                                                                           ttCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGAtt                      <  1:1300728/139‑1 (MQ=255)
                                                                                                            tCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                >  1:1638877/1‑48 (MQ=255)
                                                                                                            tCATGCAATACGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGt                                                                                                                <  2:1638877/48‑1 (MQ=255)
                                                                                                             cATTCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCGCGCAGCATCCGCCAGTTCGTTGATTGc                    <  1:1119564/139‑1 (MQ=255)
                                                                                                                  aaTACGCGGATGGCGCCATCACCAACAAAAAGGTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGCCTGGTGTCATATTTCAACAAATGCGCGATGCCCGCAGCATCCGCCAGTGCCTTGATTGccacca               <  2:1159269/139‑1 (MQ=255)
                                                                                                                     aCGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCg                            >  2:1722090/1‑123 (MQ=255)
                                                                                                                     aCGCGGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCg                            <  1:1722090/123‑1 (MQ=255)
                                                                                                                      cgcgGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGt                                                                                                          <  1:1477987/44‑1 (MQ=255)
                                                                                                                      cgcgGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGt                                                                                                          >  2:1477987/1‑44 (MQ=255)
                                                                                                                      cgcgGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaaa                                                                                         <  2:1677860/61‑1 (MQ=255)
                                                                                                                      cgcgGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaaa                                                                                         >  1:1677860/1‑61 (MQ=255)
                                                                                                                      cgcgGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGcc                                                                                   <  2:703156/67‑1 (MQ=255)
                                                                                                                      cgcgGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGcc                                                                                   >  1:703156/1‑67 (MQ=255)
                                                                                                                      cgcgGATGGCGCCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGt           >  2:1494354/1‑139 (MQ=255)
                                                                                                                              gcgcCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTc                                                                       <  1:600580/71‑1 (MQ=255)
                                                                                                                              gcgcCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGccacc                >  1:1502561/1‑126 (MQ=255)
                                                                                                                              gcgcCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGccacc                <  2:1502561/126‑1 (MQ=255)
                                                                                                                               cGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTAATTTCCGc  <  1:1617405/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                        
CCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTAATTTCCGC  >  NC_000913/3073345‑3073610

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