Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 11,087 C→G 22.1% R90R (CGG→CGC
A86A (GCC→GCG
yaaW ←
htgA →
conserved hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE11,0870CG22.1% 168.6 / 40.7 95R90R (CGG→CGC
A86A (GCC→GCG
yaaW
htgA
conserved hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (52/22);  new base G (8/13);  total (60/35)
Fisher's exact test for biased strand distribution p-value = 1.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCT  >  minE/11020‑11153
                                                                   |                                                                  
tttCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGcct                                                                   >  1:1463419/1‑68 (MQ=255)
tttCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGccc                                                                   >  1:623581/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:3344729/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:2300272/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:2261246/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:2208120/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:3434077/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:878052/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:3160892/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGa                                                                 >  1:3189057/1‑54 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:201185/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:2446105/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:2768854/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:2195901/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:347534/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:2178091/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:3008372/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1477391/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1034772/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1936005/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1130934/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:778581/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:3130352/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1658533/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1648658/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1626836/1‑69 (MQ=255)
                 tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:1545198/1‑69 (MQ=255)
                 tctTTGTCAGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:3032722/1‑69 (MQ=255)
                                aGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTttgttg                                   >  1:1527849/1‑69 (MQ=255)
                                 gCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTttgttg                                   <  1:3170035/68‑1 (MQ=255)
                                    tCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCTCGCAGTTTGTTGGCGa                               <  1:1548328/69‑1 (MQ=255)
                                    tCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               <  1:3389002/69‑1 (MQ=255)
                                    tCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               <  1:724731/69‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAt                              <  1:2329885/67‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAt                              <  1:667563/67‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAt                              <  1:2985085/67‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAt                              <  1:2869183/67‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAt                              <  1:2217514/67‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAt                              <  1:2682520/67‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATa                             <  1:2145746/68‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATa                             <  1:410873/68‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATa                             <  1:3007267/68‑1 (MQ=255)
                                       aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATa                             <  1:2717714/68‑1 (MQ=255)
                                       aTCGATTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAt                              <  1:1570602/67‑1 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTttg                                      >  1:511755/1‑58 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:219174/1‑69 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:637469/1‑69 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:3405747/1‑69 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:2634487/1‑69 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:2611307/1‑69 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:2241975/1‑69 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:403724/1‑69 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAc                            <  1:2189579/68‑1 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAc                            <  1:2136886/68‑1 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAc                            <  1:2866756/68‑1 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           <  1:2394751/69‑1 (MQ=255)
                                        tCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           <  1:2547929/69‑1 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGGCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTAt                         >  1:2472262/1‑69 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAAt                                                           >  1:1809744/1‑35 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCgtgt                                                  >  1:354954/1‑44 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:2526924/1‑67 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:2827182/1‑67 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:438900/1‑67 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:1666371/1‑67 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:539507/1‑67 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:1471584/1‑67 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:2341797/1‑67 (MQ=255)
                                          gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                           >  1:3368706/1‑67 (MQ=255)
                                             ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               >  1:3424408/1‑60 (MQ=255)
                                             ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               >  1:1106263/1‑60 (MQ=255)
                                             ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               >  1:3040338/1‑60 (MQ=255)
                                             ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               >  1:2907090/1‑60 (MQ=255)
                                             ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               >  1:2805360/1‑60 (MQ=255)
                                             ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               >  1:808649/1‑60 (MQ=255)
                                             ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGa                               >  1:2268151/1‑60 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATTCTATCGCCACCaa                <  1:1295399/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:3295027/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:1764769/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:9577/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:594187/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:2426419/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:2471979/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:253044/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:346843/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:2809044/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:2812658/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:768063/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:2390216/69‑1 (MQ=255)
                                                   cATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:747576/69‑1 (MQ=255)
                                                    aTCGAGCAAAATGGCGCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                <  1:980074/68‑1 (MQ=255)
                                                       gAGCAAAATGGCCCGATACAATTTAGCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATg            >  1:1503687/1‑69 (MQ=255)
                                                       gAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATg            >  1:3318932/1‑69 (MQ=255)
                                                       gAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATg            >  1:3086162/1‑69 (MQ=255)
                                                       gAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATg            >  1:2712594/1‑69 (MQ=255)
                                                       gAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATg            >  1:838456/1‑69 (MQ=255)
                                                       gAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATg            >  1:2056371/1‑69 (MQ=255)
                                                                 gCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtct  >  1:1964505/1‑69 (MQ=255)
                                                                   |                                                                  
TTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCT  >  minE/11020‑11153

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: