Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,514,597 G→A 41.3% G556G (GGC→GGT recG ← ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,514,5970GA41.3% 42.2 / 74.7 104G556G (GGC→GGTrecGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base A (18/25);  total (46/58)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 106 total aligned reads displayed.

TAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATCACCGC  >  minE/2514531‑2514665
                                                                  |                                                                    
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:80186/69‑1 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:769063/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:1777714/69‑1 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:2012991/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:2220974/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:2256756/69‑1 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:3266814/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:2812299/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:3192899/69‑1 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:879179/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:296816/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:921718/69‑1 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    >  1:1011731/1‑69 (MQ=255)
tAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATAGACTGGCGTTAGGCACATCAACGCCg                                                                   >  1:2954117/1‑69 (MQ=255)
 aaCTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:2705857/68‑1 (MQ=255)
 aaCTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:3296166/68‑1 (MQ=255)
 aaCTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:1018653/68‑1 (MQ=255)
         cagacccagacGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAAt                                                           <  1:1129149/69‑1 (MQ=255)
         cagacccagacGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAAt                                                           <  1:2197133/69‑1 (MQ=255)
         cagacccagacGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAAt                                                           <  1:2183784/69‑1 (MQ=255)
         cagacccagacGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAAt                                                           <  1:1068086/69‑1 (MQ=255)
         cagacccagacGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAAt                                                           <  1:1678317/69‑1 (MQ=255)
                      tCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCa                                              <  1:2996902/69‑1 (MQ=255)
                      tCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCa                                              <  1:54155/69‑1 (MQ=255)
                      tCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCa                                              <  1:614129/69‑1 (MQ=255)
                      tCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCa                                              <  1:1515328/69‑1 (MQ=255)
                      tCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCa                                              <  1:1457361/69‑1 (MQ=255)
                           gTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACcag                                         <  1:3103115/69‑1 (MQ=255)
                           gTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACcag                                         <  1:2770274/69‑1 (MQ=255)
                           gTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACcag                                         <  1:436851/69‑1 (MQ=255)
                           gTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACcag                                         <  1:1606735/69‑1 (MQ=255)
                              ttCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGt                                                     <  1:1220642/54‑1 (MQ=255)
                              ttCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGt                                                     <  1:549214/54‑1 (MQ=255)
                              ttCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGt                                                     <  1:2134221/54‑1 (MQ=255)
                                cGATAATCATCAGACTGGCGTTAGGCACATCAACGcc                                                                    <  1:3078967/37‑1 (MQ=255)
                                 gATAATCATCAGACTGGCGGTAGGCACATCAACGCCGACTTCAAt                                                           <  1:2195404/45‑1 (MQ=255)
                                   tAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGgtgt                                  <  1:2128540/68‑1 (MQ=255)
                                   tAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGgtgt                                  <  1:43133/68‑1 (MQ=255)
                                   tAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGgtgt                                  <  1:2608189/68‑1 (MQ=255)
                                   tAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGgtgt                                  <  1:2316055/68‑1 (MQ=255)
                                    aaTCATCAGACTGGGGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1993005/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGTCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1861983/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1576716/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:882937/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:782423/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:2945623/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:73290/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1899270/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:81118/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1258622/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:278865/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:653986/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1277835/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1355578/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:2669736/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:329277/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1597566/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:2597497/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1421200/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:2085045/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:1579426/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:791456/50‑3 (MQ=255)
                                    aaTCATCAGACTGGCGTTAGGCACATCAACACCGACTTCAATAACGGtag                                                   <  1:2255798/50‑3 (MQ=255)
                                        atcaGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCa                            >  1:944491/1‑69 (MQ=255)
                                           aGACTGGCGTTAGGCCCATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCt                         >  1:3172694/1‑69 (MQ=255)
                                           aGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCt                         >  1:2877837/1‑69 (MQ=255)
                                           aGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCt                         >  1:2873913/1‑69 (MQ=255)
                                           aGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCt                         >  1:2585643/1‑69 (MQ=255)
                                           aGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCt                         >  1:1577278/1‑69 (MQ=255)
                                           aGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCt                         >  1:1967582/1‑69 (MQ=255)
                                                      aGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAacga              >  1:2205856/1‑69 (MQ=255)
                                                      aGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAacga              >  1:458976/1‑69 (MQ=255)
                                                      aGGCACATCAACGCCGACTTCAATAACGGTTGTGGAACCAGCAGGTGTAGCTCACCTTGTTTAAacgac             >  1:203564/1‑69 (MQ=255)
                                                           caTCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGcc          >  1:2323382/1‑68 (MQ=255)
                                                             tCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATc       >  1:1342670/1‑69 (MQ=255)
                                                             tCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATc       >  1:1237275/1‑69 (MQ=255)
                                                             tCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATc       >  1:1146665/1‑69 (MQ=255)
                                                             tCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATc       >  1:905446/1‑69 (MQ=255)
                                                             tCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATc       >  1:2846337/1‑69 (MQ=255)
                                                             tCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCCTGTTTAAACGACGCCATc       >  1:3217944/1‑69 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:1077283/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:1715280/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:850220/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:852239/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:856969/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:1104387/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:764688/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:186713/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:71997/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:598103/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:2183273/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:535066/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:2424070/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:414900/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:3370996/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:2662940/1‑56 (MQ=255)
                                                              cAACACCGACTTCAATAACGGTAGTGGCAACCAGCAGGTGTAGATc                             >  1:1618887/1‑46 (MQ=255)
                                                              cAACACCGACTTCAATAAAGGTAGTGGCAACCAGCAGGTGTAGTTCACCTTGTTTa                   >  1:3093716/1‑56 (MQ=255)
                                                                  gCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATCACCGc  >  1:1194439/1‑69 (MQ=255)
                                                                  gCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATCACCGc  >  1:272011/1‑69 (MQ=255)
                                                                  |                                                                    
TAACTGCGCCAGACCCAGACGCTCCGGGTTTTCGATAATCATCAGACTGGCGTTAGGCACATCAACGCCGACTTCAATAACGGTTGTGGCAACCAGCAGGTGTAGCTCACCTTGTTTAAACGACGCCATCACCGC  >  minE/2514531‑2514665

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