Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,187,451 C→G E259D (GAG→GAC)  U712_10960 ← Prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,187,4510CG98.1% 167.9 / NA 53E259D (GAG→GAC) U712_10960Prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (22/30);  minor base T (1/0);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATAATCCTTGGCCGTTTTAAAATAACGAAATAGATTTTCGGCATCTTCTTGTTCTACGTATGTCCGGATCGTCTCAATTTCACGAATCCACTCTTCAAAACGGTCCAAGATTTTATCTTTATTATGTAATAGAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCCCTGAACCCTCCGGCAGCAAAACGT  >  CP006881/2187319‑2187583
                                                                                                                                    |                                                                                                                                    
cacaTATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAAt                                                                                                                                 <  4:360224/138‑1 (MQ=255)
 acaTATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAAttt                                                                                                                               <  6:327807/139‑1 (MQ=255)
        aaTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCAGTTTTAAAATAACGAAATAGATTTTCCGCATCTTATTGTTGTACATATGTACGGATCTTTTCAATTTAACCAAt                                                                                                                        >  7:41020/1‑139 (MQ=255)
        aaTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCGTGGCTGTTTTTAAATAACGAAATAGATTTTCCGCATCTTCTTGTTGTTCATATGGACGGATCTTGTCAATTTAACGAAt                                                                                                                        >  2:296322/1‑139 (MQ=255)
            atAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCa                                                                                                                     <  5:374190/137‑1 (MQ=255)
                    cAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTAGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGAATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaa                                                                                                            <  7:180464/139‑1 (MQ=255)
                     aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  4:402303/139‑1 (MQ=255)
                     aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  4:62982/139‑1 (MQ=255)
                     aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  7:274995/139‑1 (MQ=255)
                     aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  1:137171/139‑1 (MQ=255)
                     aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCCTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  3:250510/139‑1 (MQ=255)
                     aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTAGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  6:298871/139‑1 (MQ=255)
                                     tGCCGAAGCGGCAGCCCGTCGCGGTAAACCTTTGGGGTTTTAAAATAACGAAATAGATTTTACGGATCGTCTTGTTGTACATATGTACGGATGTTGTCAATTTCACAAATCCACTAATCAAAACCGACTAAAAttttt                                                                                            >  2:126366/1‑137 (MQ=255)
                                       ccGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTAt                                                                                           <  5:237864/137‑1 (MQ=255)
                                       ccGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGGTTTAAAATAACGAAATAGATTTTCCGCAGCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTAt                                                                                           <  2:74299/137‑1 (MQ=255)
                                          aaGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTtta                                                                                      <  4:285795/139‑1 (MQ=255)
                                            gCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTGTACATATGTTCGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTTTCtttttt                                                                                    >  5:325203/1‑139 (MQ=255)
                                             cGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTttatta                                                                                   <  2:376897/139‑1 (MQ=255)
                                                       tCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           >  3:437799/1‑105 (MQ=255)
                                                       tCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  4:437799/105‑1 (MQ=255)
                                                          cgATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAatat                                                                      <  8:17636/139‑1 (MQ=255)
                                                                  cTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCca                                                               <  1:271882/138‑1 (MQ=255)
                                                                     ggCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACAtt                                                           >  3:58536/1‑139 (MQ=255)
                                                                     ggCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACAtt                                                           >  5:144000/1‑139 (MQ=255)
                                                                          ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           <  3:282807/86‑1 (MQ=255)
                                                                          ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                                                           >  4:282807/1‑86 (MQ=255)
                                                                          ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCca                                                               <  3:267273/130‑1 (MQ=255)
                                                                          ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCca                                                               >  4:267273/1‑130 (MQ=255)
                                                                          ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCgg                                                      >  8:39927/1‑139 (MQ=255)
                                                                                aaTAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGgctgct                                                <  1:18126/139‑1 (MQ=255)
                                                                                aaTAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGgctgct                                                <  4:367735/139‑1 (MQ=255)
                                                                                    aCGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCTCTCATCAAAAGGTTCTAAGATTTTATCTTTATTATGTAGTAAAATATCCAGCCACATTGCCGGGAtgattgat                                            >  5:6836/1‑139 (MQ=255)
                                                                                    aCGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGAt                                            >  3:303839/1‑139 (MQ=255)
                                                                                       aaaTAGATTTTCCTCATCTTCTTGTTCTACTTATGTACTGATCTTGTCAATTTCTCGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCa                                         <  7:224085/139‑1 (MQ=255)
                                                                                        aaaaGATTTTCCGCATCTTCTTGATCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCaa                                        <  1:191320/136‑1 (MQ=255)
                                                                                            gaatttCCGCATCTTCTTGTTCTAAATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCt                                    <  3:402472/136‑1 (MQ=255)
                                                                                                  ccgCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTaa                               <  6:379266/136‑1 (MQ=255)
                                                                                                  ccgCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTaa                               <  5:399220/136‑1 (MQ=255)
                                                                                                   cgCATCTTCTTGTGCTACATATGTACGGATCTTGTAAATTTCACTAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTAATATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAat                              <  1:386309/137‑1 (MQ=255)
                                                                                                                  tACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCt                                    >  4:395136/1‑117 (MQ=255)
                                                                                                                  tACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCt                                    <  3:395136/117‑1 (MQ=255)
                                                                                                                   acatATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATTTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCTCTGAACCCGCCgg             >  2:437060/4‑139 (MQ=255)
                                                                                                                      tATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCTCTGAACCCGCCGgcag          >  2:71282/1‑139 (MQ=255)
                                                                                                                         gtacGGATCTTGTCAATTTCACGAATCCACTCATCAAACCGGTCTAATATTTTATCTTTATTATGTAAGAAAATATCCCGCCACATTGCCGGGCTGCTTGAGGCAAGCCAGGTAATATCTCTGAAACCGCCGgcagca        >  2:177696/4‑138 (MQ=255)
                                                                                                                         gtacGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGgctgct                                                <  3:366483/95‑1 (MQ=255)
                                                                                                                         gtacGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGgctgct                                                >  4:366483/4‑98 (MQ=255)
                                                                                                                         gtacGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCTCTGAACCCGCCGgcagca        >  3:329125/4‑138 (MQ=255)
                                                                                                                          tacGGATCTTGTCAATTTCACGAATCCACTCATCAAAACTGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGc                                          <  3:85168/101‑1 (MQ=255)
                                                                                                                          tacGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGc                                          >  4:85168/3‑103 (MQ=255)
                                                                                                                              gATCTTGTCAATTTCACGATTCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCTCTGAACCCGCCGGCAGCAAAACgc  >  5:433136/1‑138 (MQ=255)
                                                                                                                              gATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCTCTGAACCCGCCGGCAGCAAAACgc  >  8:413831/1‑138 (MQ=255)
                                                                                                                              gATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCTCTGAACCCGCCGGCAGCAAAACgc  >  6:392385/1‑138 (MQ=255)
                                                                                                                              gATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCTCTGAACCCGCCGGCAGCAAAACgc  >  4:219752/1‑138 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATAATCCTTGGCCGTTTTAAAATAACGAAATAGATTTTCGGCATCTTCTTGTTCTACGTATGTCCGGATCGTCTCAATTTCACGAATCCACTCTTCAAAACGGTCCAAGATTTTATCTTTATTATGTAATAGAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAATATCCCTGAACCCTCCGGCAGCAAAACGT  >  CP006881/2187319‑2187583

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑