Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 856,437 T→C 39.4% P311P (CCA→CCG trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE856,4370TC39.4% 29.9 / 47.4 71P311P (CCA→CCGtrpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (30/13);  new base C (22/6);  total (52/19)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCC  >  minE/856371‑856496
                                                                  |                                                           
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:510205/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:2032535/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:3137055/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:1588489/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:210237/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:3110377/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:1553472/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:1549479/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:1032013/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:447538/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:2370596/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:2383016/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:680150/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:693961/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:2821426/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:1309331/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:2956775/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:2971863/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:2979920/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:1158312/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCATGTTGCgg                                                           >  1:981557/2‑69 (MQ=255)
caGGGCTTCATCGTCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTGAGATACGCAAGTTGCgg                                                           >  1:851856/2‑69 (MQ=255)
   ggCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGcc                                                        >  1:3029041/1‑69 (MQ=255)
   ggCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGcc                                                        >  1:3136076/1‑69 (MQ=255)
   ggCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGcc                                                        >  1:1809277/1‑69 (MQ=255)
   ggCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGcc                                                        >  1:455075/1‑69 (MQ=255)
   ggCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGcc                                                        >  1:155735/1‑69 (MQ=255)
                  aaTAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATc                                         <  1:3158044/69‑1 (MQ=255)
                  aaTAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATc                                         <  1:3413117/69‑1 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:1995100/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:1725885/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:1647311/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:3340487/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:1168585/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:1235462/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:2322799/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:2464905/1‑69 (MQ=255)
                      gACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt                                     >  1:1335781/1‑69 (MQ=255)
                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGc                                                         >  1:2273187/1‑44 (MQ=255)
                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCcgg                                 >  1:2976005/1‑68 (MQ=255)
                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCcgg                                 >  1:1940657/1‑68 (MQ=255)
                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCcgg                                 >  1:1269648/1‑68 (MQ=255)
                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCcgg                                 >  1:543814/1‑68 (MQ=255)
                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCcgg                                 >  1:3395623/1‑68 (MQ=255)
                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCcg                                  >  1:3160716/1‑67 (MQ=255)
                              aTCAGCGCGTCCAGTGCTGTTGAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGga                             <  1:1969739/69‑1 (MQ=255)
                              aTCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGga                             <  1:2204543/69‑1 (MQ=255)
                              aTCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGga                             <  1:1090500/69‑1 (MQ=255)
                              aTCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGga                             <  1:1541834/69‑1 (MQ=255)
                              aTCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGga                             <  1:2030156/69‑1 (MQ=255)
                                          aGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGa                 <  1:3424642/69‑1 (MQ=255)
                                          aGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGa                 <  1:895556/69‑1 (MQ=255)
                                          aGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGa                 <  1:1365888/69‑1 (MQ=255)
                                          aGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGa                 <  1:1616193/69‑1 (MQ=255)
                                          aGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGa                 <  1:31462/69‑1 (MQ=255)
                                          aGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGa                 <  1:194927/69‑1 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATc                                         >  1:3331445/1‑41 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATc                                         >  1:1533137/1‑41 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATc                                         >  1:1490126/1‑41 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGAttc              >  1:3104848/1‑68 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGAttc              >  1:948753/1‑68 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGAttc              >  1:2147201/1‑68 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGAttc              >  1:1303112/1‑68 (MQ=255)
                                              cTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGAttc              >  1:1401377/1‑68 (MQ=255)
                                                      aTACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTg     >  1:2100700/1‑69 (MQ=255)
                                                         cgcaTGTTGCGTGCCGACGGACGTGAAATCCAGCCCGGCAGAAATGGAGTAAGACTCTTCAATTTGccc  <  1:2091001/69‑1 (MQ=255)
                                                         cgcaTGTTGCGGGCCGACGGACGGGAAATCCAGCCCGGCATAAATGGAGTAAGACTCTTCAATTTGccc  <  1:572059/69‑1 (MQ=255)
                                                         cgcaTGTTGCGGGCCGACGGACGGGAAATCCAGCCCGGCAGAAATGGAGTAAGACTCTTCAATTTGccc  <  1:1350278/69‑1 (MQ=255)
                                                         cgcaTGTTGCGGGCCGACGGACGGGAAATCCAGCCCGGCAGAAATGGAGTAAGACTCTTCAATTTGccc  <  1:1416063/69‑1 (MQ=255)
                                                         cgcaTGTTGCGGGCCGACGGACGGGAAATCCAGCCCGGCAGAAATGGAGTAAGACTCTTCAATTTGccc  <  1:2991093/69‑1 (MQ=255)
                                                         cgcaTGTTGCGGGCCGACGGACGGGAAATCCAGCCCGGCAGAAATGGAGTAAGACTCTTCAATTTGccc  <  1:1557336/69‑1 (MQ=255)
                                                                  |                                                           
AAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCC  >  minE/856371‑856496

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