Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,829,870 T→A 17.2% L149H (CTC→CAC) ‡ pphB → serine/threonine‑specific protein phosphatase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,829,8700TA17.2% 112.8 / 15.9 58L149H (CTC→CAC) ‡pphBserine/threonine‑specific protein phosphatase 2
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base A (5/5);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

ATCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAA  >  minE/1829809‑1829936
                                                             |                                                                  
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:1000253/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:3354747/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:1210369/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:1689629/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:3181492/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:2535826/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:2730288/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:2863092/1‑68 (MQ=255)
aTCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTg                                                               >  1:2169174/1‑67 (MQ=255)
  cACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:3059921/1‑66 (MQ=255)
  cACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:2485230/1‑66 (MQ=255)
  cACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTgg                                                              >  1:3162262/1‑66 (MQ=255)
     tAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGa                                                        <  1:1245858/69‑1 (MQ=255)
     tAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGa                                                        <  1:61972/69‑1 (MQ=255)
     tAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGa                                                        <  1:1096971/69‑1 (MQ=255)
                         ccATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTa                                     <  1:1024594/68‑1 (MQ=255)
                         ccATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTa                                     <  1:1949555/68‑1 (MQ=255)
                         ccATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTa                                     <  1:1945561/68‑1 (MQ=255)
                         ccATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTa                                     <  1:1873970/68‑1 (MQ=255)
                         ccATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTa                                     <  1:1844049/68‑1 (MQ=255)
                                 caTGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGcc                             >  1:1147105/1‑68 (MQ=255)
                                 caTGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCt                            >  1:980497/1‑69 (MQ=255)
                                 caTGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCt                            >  1:1325913/1‑69 (MQ=255)
                                 caTGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCt                            >  1:3262612/1‑69 (MQ=255)
                                  aTGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTg                           <  1:944500/69‑1 (MQ=255)
                                         ttATCCGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCg                    >  1:1756966/1‑69 (MQ=255)
                                         ttATCCGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCg                    >  1:1768110/1‑69 (MQ=255)
                                         ttATCCGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCg                    >  1:307941/1‑69 (MQ=255)
                                         ttATCCGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCg                    >  1:1632986/1‑69 (MQ=255)
                                         ttATCCGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCg                    >  1:67493/1‑69 (MQ=255)
                                          tATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTgg                               <  1:3238235/57‑1 (MQ=255)
                                          tATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTgg                               <  1:1268467/57‑1 (MQ=255)
                                          tATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTgg                               <  1:3016252/57‑1 (MQ=255)
                                          tATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTgg                               <  1:3012315/57‑1 (MQ=255)
                                          tATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTgg                               <  1:2890198/57‑1 (MQ=255)
                                          tATCCGGGGAGTGAATATCTATTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTgg                               <  1:1107266/57‑1 (MQ=255)
                                             ccGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtgt                 <  1:3393192/68‑1 (MQ=255)
                                             ccGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtgt                 <  1:989294/68‑1 (MQ=255)
                                             ccGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtgt                 <  1:1395852/68‑1 (MQ=255)
                                             ccGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtgt                 <  1:1937453/68‑1 (MQ=255)
                                             ccGGGGAAGGAGTATCAGTTTGGCAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtgt                 <  1:2413416/68‑1 (MQ=255)
                                              cGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtg                  >  1:1780682/1‑66 (MQ=255)
                                              cGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtg                  >  1:3302957/1‑66 (MQ=255)
                                                  gAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGaa           <  1:685486/69‑1 (MQ=255)
                                                       aaTATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCg       >  1:1843075/1‑68 (MQ=255)
                                                       aaTATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCg       >  1:911809/1‑68 (MQ=255)
                                                       aaTATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCg       >  1:2532395/1‑68 (MQ=255)
                                                       aaTATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCg       >  1:2529794/1‑68 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:2687377/69‑1 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:2924565/69‑1 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:473681/69‑1 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:413812/69‑1 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    >  1:3336583/1‑69 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:320711/69‑1 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:1677325/69‑1 (MQ=255)
                                                         tatCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:1325010/69‑1 (MQ=255)
                                                          atCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:2826718/68‑1 (MQ=255)
                                                          atCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt    <  1:2809731/68‑1 (MQ=255)
                                                           tctcTTTGGTAAAGAAATAGCGGAGAGCGAATTAc                                    <  1:3385219/35‑1 (MQ=255)
                                                           tctcTTTGGTAAAGAAATAGCGGAGAGCGAATTAc                                    <  1:2086038/35‑1 (MQ=255)
                                                           tctcTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTaa  >  1:2516272/1‑69 (MQ=255)
                                                           tctcTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTaa  >  1:231469/1‑69 (MQ=255)
                                                             |                                                                  
ATCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAA  >  minE/1829809‑1829936

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