Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,402 A→C 21.2% T194T (ACA→ACC yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,4020AC21.2% 125.3 / 32.9 80T194T (ACA→ACCyadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (25/38);  new base C (7/10);  total (32/48)
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.83e-01

GCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCAT  >  minE/130336‑130467
                                                                  |                                                                 
gCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATc                                                                 <  1:1286490/69‑1 (MQ=255)
 cTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATc                                                                 >  1:1605101/1‑68 (MQ=255)
 cTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATc                                                                 >  1:3272445/1‑68 (MQ=255)
 cTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATc                                                                 <  1:2905317/68‑1 (MQ=255)
 cTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATc                                                                 >  1:2488748/1‑68 (MQ=255)
   aaCCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCgg                                                               <  1:2172833/68‑1 (MQ=255)
         tCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACgcgc                                                        >  1:2743308/1‑69 (MQ=255)
         tCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACgcgc                                                        >  1:1512310/1‑69 (MQ=255)
                atTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAAt                                                 >  1:1240565/1‑69 (MQ=255)
                atTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAAt                                                 >  1:2614673/1‑69 (MQ=255)
                atTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAAt                                                 >  1:3326804/1‑69 (MQ=255)
                atTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAAt                                                 >  1:3380479/1‑69 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAAt                                                 >  1:3391754/1‑68 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAAt                                                 >  1:1534354/1‑68 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATa                                                >  1:272956/1‑69 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATa                                                >  1:3235938/1‑69 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATa                                                >  1:2522210/1‑69 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATa                                                >  1:2346893/1‑69 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATa                                                >  1:349019/1‑69 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATa                                                >  1:1459874/1‑69 (MQ=255)
                 ttCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATa                                                >  1:1150690/1‑69 (MQ=255)
                         gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCAATAACCAGATg                                        >  1:1584683/1‑69 (MQ=255)
                         gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCAATAACCAGATg                                        >  1:1667882/1‑69 (MQ=255)
                         gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCAATAACCAGATg                                        >  1:1300801/1‑69 (MQ=255)
                         gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCAATAACCAGATg                                        >  1:508481/1‑69 (MQ=255)
                         gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCAATAACCAGATg                                        >  1:1243421/1‑69 (MQ=255)
                         gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCAATAACCAGATg                                        >  1:591785/1‑69 (MQ=255)
                         gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCAATAACCAGATg                                        >  1:728123/1‑69 (MQ=255)
                          acgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATgg                                       <  1:2420262/69‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:2343184/52‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:954500/52‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:2284677/52‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:495383/52‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:708823/52‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:864404/52‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:1872742/52‑1 (MQ=255)
                               gaaAACACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:1083121/52‑1 (MQ=255)
                                aaaaCACCCGTTTTCGCCATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:581207/51‑1 (MQ=255)
                                  aaCACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTg                               <  1:2026774/69‑1 (MQ=255)
                                  aaCACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTg                               <  1:3404537/69‑1 (MQ=255)
                                  aaCACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTg                               <  1:255573/69‑1 (MQ=255)
                                   acacCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTg                               <  1:2247838/68‑1 (MQ=255)
                                   acacCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTg                               <  1:1386298/68‑1 (MQ=255)
                                    cacCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTg                               <  1:3400682/67‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:1740171/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:1267252/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:1503066/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:1677130/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:17403/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:3217443/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:1076127/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:3096555/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:3089362/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:1875435/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:2184431/69‑1 (MQ=255)
                                     acCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGt                            <  1:2102066/69‑1 (MQ=255)
                                         gTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACa                        <  1:294025/69‑1 (MQ=255)
                                            ttCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAggg                     <  1:2616321/69‑1 (MQ=255)
                                            ttCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAggg                     <  1:2581403/69‑1 (MQ=255)
                                            ttCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAggg                     <  1:2170058/69‑1 (MQ=255)
                                              cGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAggg                     <  1:456703/67‑1 (MQ=255)
                                                ccATACTTCGACGACCACCTCGGTACGCGCTTTCa                                                   <  1:623705/35‑1 (MQ=255)
                                                     cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGa            >  1:1000150/1‑69 (MQ=255)
                                                     cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGa            >  1:822121/1‑69 (MQ=255)
                                                     cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGa            >  1:1087385/1‑69 (MQ=255)
                                                     cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGa            >  1:68552/1‑69 (MQ=255)
                                                     cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGa            >  1:1153258/1‑69 (MQ=255)
                                                     cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGa            >  1:348199/1‑69 (MQ=255)
                                                     cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGa            >  1:2263344/1‑69 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:350186/68‑1 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:2134426/68‑1 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:1621366/68‑1 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:1981798/68‑1 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:3130204/68‑1 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:2341435/68‑1 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:92007/68‑1 (MQ=255)
                                                       tCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:2965260/68‑1 (MQ=255)
                                                            gacCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAg     <  1:1468562/69‑1 (MQ=255)
                                                            gacCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAg     <  1:2668608/69‑1 (MQ=255)
                                                             acCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGac           <  1:2133113/62‑1 (MQ=255)
                                                               cacaTCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCAt  >  1:42336/1‑69 (MQ=255)
                                                                  |                                                                 
GCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCAT  >  minE/130336‑130467

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