Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,732,532 A→G D299D (GAT→GAC tyrA ← bifunctional chorismate mutase/prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,732,5320AG97.9% 163.4 / ‑3.7 47D299D (GAT→GACtyrAbifunctional chorismate mutase/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (25/21);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

TGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGATCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGGCGCAGT  >  NZ_CP009273/2732396‑2732667
                                                                                                                                        |                                                                                                                                       
tGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTc                                                                                                                                       >  1:1183383/1‑139 (MQ=255)
tGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTc                                                                                                                                       >  2:308731/1‑139 (MQ=255)
  cGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCt                                                                                                                                     >  2:865028/1‑139 (MQ=255)
                 aaCGCCTGCTTATCGCCCTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGcc                                                                                                                      <  1:865028/139‑1 (MQ=255)
                             tCGCCCTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGc                                                                                                            >  1:102823/1‑137 (MQ=255)
                                  cTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCa                                                                                                     <  1:1288270/139‑1 (MQ=255)
                                  cTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCa                                                                                                     <  1:983272/139‑1 (MQ=255)
                                   tGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTc                                                                                                      >  2:1288270/1‑137 (MQ=255)
                                    gCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGcc                                                                                                                      >  2:180477/1‑120 (MQ=255)
                                    gCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGcc                                                                                                                      <  1:180477/120‑1 (MQ=255)
                                     cTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCaa                                                                                                    <  2:433072/137‑1 (MQ=255)
                                          gCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGt                                                                                              >  1:689759/1‑138 (MQ=255)
                                            aaCTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTaaa                                                                                           <  1:308731/139‑1 (MQ=255)
                                                       ctcgccGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAagag                                                                                <  2:515817/139‑1 (MQ=255)
                                                                       taGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCt                                                                <  2:284837/139‑1 (MQ=255)
                                                                            aCGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGc                                                           >  1:1212530/1‑139 (MQ=255)
                                                                            aCGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGc                                                           >  2:727190/1‑139 (MQ=255)
                                                                            aCGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGc                                                           >  1:1200151/1‑139 (MQ=255)
                                                                             cGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCt                                                          >  1:621004/1‑139 (MQ=255)
                                                                              gTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCt                                                          >  2:774490/1‑138 (MQ=255)
                                                                                ttGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGaa                                                       >  1:913648/1‑139 (MQ=255)
                                                                                        cGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTttct                                               >  1:968354/1‑139 (MQ=255)
                                                                                              aTTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGcc                                         <  1:807316/139‑1 (MQ=255)
                                                                                              aTTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGcc                                         <  1:774490/139‑1 (MQ=255)
                                                                                              aTTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGcc                                         <  2:1040035/139‑1 (MQ=255)
                                                                                              aTTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGATCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGcc                                         <  2:307425/139‑1 (MQ=255)
                                                                                                     tctGACGACATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTg                                    >  1:441179/1‑137 (MQ=255)
                                                                                                        gacgacATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGcc                               >  1:991977/1‑139 (MQ=255)
                                                                                                         acgacATAATGATGTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGccc                              >  2:431492/1‑139 (MQ=255)
                                                                                                                 atgatgTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCa                       >  2:140826/1‑138 (MQ=255)
                                                                                                                   gatgTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCaa                      >  1:304401/1‑137 (MQ=255)
                                                                                                                     tgTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTttctt                                              <  1:708394/111‑1 (MQ=255)
                                                                                                                     tgTCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTttctt                                              >  2:708394/1‑111 (MQ=255)
                                                                                                                       tCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAg                >  2:295210/1‑139 (MQ=255)
                                                                                                                       tCGGCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAg                >  1:987268/1‑139 (MQ=255)
                                                                                                                         ggCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCa              <  1:727190/139‑1 (MQ=255)
                                                                                                                         ggCATAAAGCTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCa              <  1:295210/139‑1 (MQ=255)
                                                                                                                                  cTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGgcgc     <  2:1200151/139‑1 (MQ=255)
                                                                                                                                  cTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGgcgc     <  2:437837/139‑1 (MQ=255)
                                                                                                                                  cTGCGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGgcgc     <  1:325967/139‑1 (MQ=255)
                                                                                                                                     cGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCa              <  2:777257/127‑1 (MQ=255)
                                                                                                                                     cGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCa              >  1:777257/1‑127 (MQ=255)
                                                                                                                                     cGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGgcg      <  1:756926/135‑1 (MQ=255)
                                                                                                                                     cGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGgcg      >  2:756926/1‑135 (MQ=255)
                                                                                                                                     cGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGGCGCa    <  1:431492/137‑1 (MQ=255)
                                                                                                                                     cGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGGCGCa    <  2:621004/137‑1 (MQ=255)
                                                                                                                                     cGGGTCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGGCGCAGt  >  2:267256/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTCCAGCAACTCAATCGCCTCGCCGAAACGCTTATAGTAACGTTTGATTAACGCCAGATTACGCTCTGACGACATAATGATGTCGGCATAAAGCTGCGGATCCTGAGCAAACAGTCGCCCGACCATCGCCAGCTCAAGGCGGTAAATCGGCGAAGAGAGCGCCAGAAGTTGCTCAAGCTGAACATTTTCTTCTGCCAGGTGCAGCCCGTAAGCAAAAGTAGCAAAGTGGCGCAGT  >  NZ_CP009273/2732396‑2732667

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