Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,302,805 G→A 58.4% R92R (CGC→CGT dacD ← D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,302,8050GA58.4% 36.8 / 80.6 108R92R (CGC→CGTdacDD‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base A (41/22);  total (60/48)
Fisher's exact test for biased strand distribution p-value = 3.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CGCCCTCTTTCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATCCACGACATAA  >  minE/1302741‑1302871
                                                                |                                                                  
cGCCCTCTTTCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCgcgc                                                                 >  1:1612955/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:7140/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:2861069/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:2201462/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:2157544/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:1847119/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:52326/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:1865194/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:3091270/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:1051135/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTg                                                        >  1:1020003/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGt                                                         >  1:1906185/1‑67 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGt                                                         >  1:2358889/1‑67 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:3081581/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:3168443/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:1929198/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:1949236/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:3097894/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:2025466/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:1534216/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:1018223/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:1467781/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:1070985/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:2512460/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:2761074/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:815699/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:516362/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:532910/1‑68 (MQ=255)
         tCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTg                                                        >  1:756496/1‑68 (MQ=255)
             aaaCATCTGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:1507360/69‑1 (MQ=255)
             aaaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:1166150/69‑1 (MQ=255)
             aaaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:744318/69‑1 (MQ=255)
             aaaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:2773233/69‑1 (MQ=255)
             aaaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:3327293/69‑1 (MQ=255)
             aaaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:2276565/69‑1 (MQ=255)
              aaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:3030750/68‑1 (MQ=255)
              aaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:2494976/68‑1 (MQ=255)
              aaCATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAat                                                   <  1:27083/68‑1 (MQ=255)
                          gaaCCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGt                                       <  1:1545112/68‑1 (MQ=255)
                          gaaCCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGt                                       <  1:475226/68‑1 (MQ=255)
                          gaaCCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGt                                       <  1:140143/68‑1 (MQ=255)
                           aaCCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGt                                       <  1:356334/67‑1 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGccc                                                               >  1:3259432/1‑38 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:3034995/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:3297092/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:353362/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:2911754/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:2906355/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:387979/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:493327/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:929615/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:2737739/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:929756/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:1572598/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:1073580/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:2533361/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:2675106/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:1834497/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:1182438/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGt                                       >  1:2107710/1‑62 (MQ=255)
                                aCAAACACCGGATTATCTTTCGCCCACGCATCACGCCCCACGGTGACAATATCGTCTGGCGTaaaa                                   >  1:351581/1‑64 (MQ=255)
                                 cAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGAt                               <  1:1968603/69‑1 (MQ=255)
                                 cAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGAt                               <  1:612563/69‑1 (MQ=255)
                                 cAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGAt                               <  1:3095336/69‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTgg                                          >  1:2586185/1‑54 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACt                           >  1:779196/1‑69 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACt                           >  1:1961064/1‑69 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACt                           >  1:1891834/1‑69 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACt                           >  1:1869505/1‑69 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACt                           >  1:836300/1‑69 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1036666/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1147038/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1499204/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1325374/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:907613/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1204274/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1085639/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:7870/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:334913/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1088311/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1155526/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1627614/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:2537478/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:2776396/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:2323352/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:2153435/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:2108072/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:3169538/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:3416928/35‑1 (MQ=255)
                                     cacCGGATTATCTTTCGCCCACGCATCACGCCCCa                                                             <  1:1773889/35‑1 (MQ=255)
                                          gATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGt                                       <  1:1426847/52‑1 (MQ=255)
                                              aTCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAgc                  >  1:3277459/1‑69 (MQ=255)
                                                     gCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATc            <  1:300817/68‑1 (MQ=255)
                                                     gCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATc            <  1:821312/68‑1 (MQ=255)
                                                     gCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATc            <  1:2360729/68‑1 (MQ=255)
                                                     gCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATc            <  1:260965/68‑1 (MQ=255)
                                                              tCGCGCCCCACGGTGACAATATCGTTTGGCGTAATGCGATGACTATCGATAGCGCGATCCACGACATaa  <  1:2288012/69‑1 (MQ=255)
                                                              tCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATCCACGACATaa  <  1:368124/69‑1 (MQ=255)
                                                               cgcgCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATCCACGACATaa  <  1:2977341/68‑1 (MQ=255)
                                                                gcgcCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATc                        <  1:1803630/45‑1 (MQ=255)
                                                                |                                                                  
CGCCCTCTTTCAAAAACATCAGTGAAGAACCGACAAACACCGGATTATCTTTCGCCCACGCATCGCGCCCCACGGTGACAATATCGTCTGGCGTAATGCGATGACTATCGATAGCGCGATCCACGACATAA  >  minE/1302741‑1302871

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: