Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,677,674 T→A 28.5% intergenic (‑124/+25) ndk ← / ← pbpC multifunctional nucleoside diphosphate kinase, apyrimidinic endonuclease, and 3'‑phosphodiesterase/fused transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,677,6740TA28.5% 90.8 / 42.1 70intergenic (‑124/+25)ndk/pbpCmultifunctional nucleoside diphosphate kinase, apyrimidinic endonuclease, and 3'‑phosphodiesterase/fused transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (43/7);  new base A (9/11);  total (52/18)
Fisher's exact test for biased strand distribution p-value = 8.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCG  >  minE/1677608‑1677726
                                                                  |                                                    
ttATAATGGGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGtttt                                                      >  1:499679/1‑67 (MQ=255)
ttATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGtttt                                                      >  1:501978/1‑67 (MQ=255)
ttATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTa                                                     >  1:242154/1‑68 (MQ=255)
ttATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTa                                                     >  1:87949/1‑68 (MQ=255)
ttATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTa                                                     >  1:1700574/1‑68 (MQ=255)
ttATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTa                                                     >  1:3213578/1‑68 (MQ=255)
ttATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTa                                                     >  1:1493715/1‑68 (MQ=255)
   taATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGca                                                  >  1:2412538/1‑68 (MQ=255)
   taATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGca                                                  >  1:1807691/1‑68 (MQ=255)
     aTGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGcaaca                                               >  1:2605840/1‑69 (MQ=255)
     aTGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGcaaca                                               >  1:486758/1‑69 (MQ=255)
                 gCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACaa                                    >  1:1874340/1‑68 (MQ=255)
                 gCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACaa                                    >  1:3197506/1‑68 (MQ=255)
                  ccGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACaaaa                                  <  1:1715733/69‑1 (MQ=255)
                       gACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACaaaa                                  >  1:174993/1‑64 (MQ=255)
                       gACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACaaaa                                  >  1:103068/1‑64 (MQ=255)
                       gACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACaaaa                                  >  1:1276253/1‑64 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCa                                           >  1:502673/1‑50 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:2571313/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:2608655/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:2361510/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:2611400/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:2631369/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:484968/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:578824/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:591352/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:89549/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:924977/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:998326/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:1739876/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:1347888/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:1244947/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:1516320/1‑69 (MQ=255)
                            ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:114106/1‑69 (MQ=255)
                            ttGATCAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTg                        >  1:1768997/1‑69 (MQ=255)
                              gatgaAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGc                       >  1:2911595/1‑68 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTg                        >  1:2634608/1‑66 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCa                      >  1:3067123/1‑68 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCAt                     >  1:3040194/1‑69 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCAt                     >  1:1947377/1‑69 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCAt                     >  1:877961/1‑69 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCAt                     >  1:2708636/1‑69 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCAt                     >  1:738721/1‑69 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCAt                     >  1:1129687/1‑69 (MQ=255)
                               atgaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCAt                     >  1:455723/1‑69 (MQ=255)
                                 gaAGGTAAGATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:1304734/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAATAATCAGGGACAAAATAGACTATTGCATg                    <  1:1373284/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:1442507/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:1396696/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:1170142/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:2270984/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:519358/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:2712709/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:2477835/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:1041968/68‑1 (MQ=255)
                                 gaAGGTAACATTTTTTTAAAATAAACCCAGTTTAAGCAACAATCAGGGACAAAATAGACTATTGCATg                    <  1:2385201/68‑1 (MQ=255)
                                        aCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCa                                           >  1:3133769/1‑38 (MQ=255)
                                        aCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAg                                          >  1:1030701/1‑39 (MQ=255)
                                        aCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAAttt            >  1:319035/1‑69 (MQ=255)
                                        aCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAAttt            >  1:289673/1‑69 (MQ=255)
                                        aCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAAttt            >  1:1729967/1‑69 (MQ=255)
                                        aCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAAtt             >  1:504789/1‑68 (MQ=255)
                                          aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTc           <  1:3331267/68‑1 (MQ=255)
                                          aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTc           <  1:3123394/68‑1 (MQ=255)
                                          aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTc           <  1:2766049/68‑1 (MQ=255)
                                          aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTc           <  1:2064139/68‑1 (MQ=255)
                                          aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCa          >  1:2747902/1‑69 (MQ=255)
                                          aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCa          >  1:2146993/1‑69 (MQ=255)
                                          aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTAATGCATGACAAATTTc           <  1:1362472/68‑1 (MQ=255)
                                           tttttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCAc         >  1:2811241/1‑69 (MQ=255)
                                           tttttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCAc         >  1:2076212/1‑69 (MQ=255)
                                                   aaataaaatGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTcgcg  <  1:656741/68‑1 (MQ=255)
                                                                  |                                                    
TTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCG  >  minE/1677608‑1677726

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: