Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,187,436 G→A Y264Y (TAC→TAT U712_10960 ← Prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,187,4360GA97.9% 142.0 / NA 48Y264Y (TAC→TATU712_10960Prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/30);  minor base T (0/1);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATGATCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATAATCCTTGGCCGTTTTAAAATAACGAAATAGATTTTCGGCATCTTCTTGTTCTACGTATGTCCGGATCGTCTCAATTTCACGAATCCACTCTTCAAAACGGTCCAAGATTTTATCTTTATTATGTAATAGAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAAT  >  CP006881/2187302‑2187555
                                                                                                                                      |                                                                                                                       
ggATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTTAATAAACGAAATAGATTTTCCGCATCTTCTTGTTGTACATATg                                                                                                                     >  2:312066/1‑139 (MQ=255)
ggATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATg                                                                                                                     <  1:392721/139‑1 (MQ=255)
 gatgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGGTCTACATATGt                                                                                                                    >  7:34575/1‑139 (MQ=255)
  atgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGta                                                                                                                   <  2:165845/139‑2 (MQ=255)
  atgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTGGGCGGTTTTTAAATAAAGAAATAGATTTTACGCATGTTGTTGTTGTAGATATGtt                                                                                                                   >  3:285795/1‑138 (MQ=255)
  atgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTAGGCGGTTTTAAAATTAAGAAATAGATTTTACCCATGTTCTTTTTCTACATATGta                                                                                                                   >  7:17636/1‑138 (MQ=255)
   tgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGtac                                                                                                                  <  8:161475/139‑3 (MQ=255)
                 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                                                                               <  7:274995/139‑1 (MQ=255)
                                      aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                               <  4:62982/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                                                                               <  4:437799/105‑1 (MQ=255)
                                                                        tCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                               >  3:437799/1‑105 (MQ=255)
                                                                           cgATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAatat                                          <  8:17636/139‑1 (MQ=255)
                                                                                   cTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCca                                   <  1:271882/138‑1 (MQ=255)
                                                                                      ggCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACAtt                               >  5:144000/1‑139 (MQ=255)
                                                                                      ggCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACAtt                               >  3:58536/1‑139 (MQ=255)
                                                                                           ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                               >  4:282807/1‑86 (MQ=255)
                                                                                           ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa                                                                               <  3:282807/86‑1 (MQ=255)
                                                                                           ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCca                                   >  4:267273/1‑130 (MQ=255)
                                                                                           ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCca                                   <  3:267273/130‑1 (MQ=255)
                                                                                           ttttAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCgg                          >  8:39927/1‑139 (MQ=255)
                                                                                                 aaTAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGgctgct                    <  4:367735/139‑1 (MQ=255)
                                                                                                 aaTAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGgctgct                    <  1:18126/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   <  5:399220/136‑1 (MQ=255)
                                                                                                                   ccgCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCTAAGATTTTATCTTTATTATGTAATAAAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTaa   <  6:379266/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)
                                                                                                                                      |                                                                                                                       
GGATGATCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATAATCCTTGGCCGTTTTAAAATAACGAAATAGATTTTCGGCATCTTCTTGTTCTACGTATGTCCGGATCGTCTCAATTTCACGAATCCACTCTTCAAAACGGTCCAAGATTTTATCTTTATTATGTAATAGAATATCCCGCCACATTGCCGGGCTGCTTGATGCAATCCTTGTAAT  >  CP006881/2187302‑2187555

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: