Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,628 T→G 42.8% A375A (GCT→GCG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6280TG42.8% 10.5 / 56.4 56A375A (GCT→GCGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (10/22);  new base G (14/10);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 5.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCAC  >  minE/189566‑189696
                                                              |                                                                    
gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                                <  1:2247072/69‑1 (MQ=255)
gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                                <  1:951077/69‑1 (MQ=255)
gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                                <  1:3082972/69‑1 (MQ=255)
gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                                <  1:1707793/69‑1 (MQ=255)
gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                                <  1:2872832/69‑1 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:2907357/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:2192034/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:2230900/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:2366258/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:2527156/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:63278/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:2369701/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:3084641/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:3182008/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:3309554/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:357579/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:536430/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:1823945/2‑68 (MQ=255)
  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa                                                               >  1:1711881/2‑68 (MQ=255)
     aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCca                                                           <  1:214569/69‑1 (MQ=255)
     aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCca                                                           <  1:1480924/69‑1 (MQ=255)
         gTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGAGGCTTCGCTGATCCACATTa                                                      <  1:2105721/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCTCTGATCCACATTaa                                                     <  1:1236701/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:2843994/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:35972/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:1076258/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:3085861/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:858484/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:2914234/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:550792/69‑1 (MQ=255)
           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:2319149/69‑1 (MQ=255)
                           tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     <  1:77306/69‑1 (MQ=255)
                           tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     <  1:1027204/69‑1 (MQ=255)
                           tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     <  1:3369657/69‑1 (MQ=255)
                           tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     <  1:1110660/69‑1 (MQ=255)
                           tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     <  1:2703795/69‑1 (MQ=255)
                           tCGACCCACAAAATCCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     <  1:1825723/69‑1 (MQ=255)
                                            ggcggcGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGc                    <  1:3247339/69‑1 (MQ=255)
                                            ggcggcGTTATCGCAGGCGTCGCTGATCCACATTaa                                                     <  1:2886920/36‑1 (MQ=255)
                                                  gTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTat              <  1:2280492/69‑1 (MQ=255)
                                                     aTCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     >  1:2957148/1‑43 (MQ=255)
                                                     aTCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     >  1:2867877/1‑43 (MQ=255)
                                                     aTCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                     >  1:2468252/1‑43 (MQ=255)
                                                       cGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg          >  1:1469601/1‑68 (MQ=255)
                                                       cGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg          >  1:1420048/1‑68 (MQ=255)
                                                       cGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg          >  1:2620747/1‑68 (MQ=255)
                                                       cGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg          >  1:2319810/1‑68 (MQ=255)
                                                       cGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg          >  1:722172/1‑68 (MQ=255)
                                                       cGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg          >  1:2492393/1‑68 (MQ=255)
                                                       cGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgat           >  1:1074023/1‑67 (MQ=255)
                                                          aGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCAt       <  1:540195/68‑1 (MQ=255)
                                                          aGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCAt       <  1:2625063/68‑1 (MQ=255)
                                                          aGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCAt       <  1:1351244/68‑1 (MQ=255)
                                                          aGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCAt       <  1:1305529/68‑1 (MQ=255)
                                                          aGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCAt       <  1:941116/68‑1 (MQ=255)
                                                          aGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCAt       <  1:1283640/68‑1 (MQ=255)
                                                              ttCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccac  <  1:2580684/69‑1 (MQ=255)
                                                              ttCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccac  <  1:2075475/69‑1 (MQ=255)
                                                              |                                                                    
GCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCAC  >  minE/189566‑189696

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