Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,355,216 T→C 5.7% Y328C (TAT→TGT)  yihE ← predicted kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,355,2160TC 5.7% 248.0 / 6.1 87Y328C (TAT→TGT) yihEpredicted kinase
Reads supporting (aligned to +/- strand):  ref base T (36/46);  new base C (2/3);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AACGCTAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTGT  >  minE/2355152‑2355284
                                                                |                                                                    
aacgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACAt                                                                  <  1:2159972/69‑1 (MQ=255)
 acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAgg                                                               >  1:3333082/1‑71 (MQ=255)
 acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCAATTAACACCTatg                                                               <  1:3848179/71‑3 (MQ=255)
 acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCAATTAACACCTatg                                                               <  1:1484006/71‑3 (MQ=255)
 acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCAATTAACACCTatg                                                               <  1:780513/71‑3 (MQ=255)
  cgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGgt                                                              <  1:158712/71‑1 (MQ=255)
       aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:3674646/69‑1 (MQ=255)
       aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:275768/69‑1 (MQ=255)
       aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:3124976/69‑1 (MQ=255)
       aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTaa                                                         <  1:1216852/71‑1 (MQ=255)
       aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTaa                                                         <  1:3833444/71‑1 (MQ=255)
       aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTaa                                                         <  1:543450/71‑1 (MQ=255)
          cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:3553002/66‑1 (MQ=255)
          cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:2381166/66‑1 (MQ=255)
          cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:3166888/66‑1 (MQ=255)
          cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:2773391/66‑1 (MQ=255)
          cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:2092275/66‑1 (MQ=255)
          cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                           <  1:2087683/66‑1 (MQ=255)
             aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                    >  1:3778880/1‑70 (MQ=255)
             aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                    <  1:1972541/70‑1 (MQ=255)
             aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                    >  1:767628/1‑70 (MQ=255)
             aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                    >  1:1531569/1‑70 (MQ=255)
                     agcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGtt                                           >  1:2889840/1‑71 (MQ=255)
                     agcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGt                                            >  1:1308963/1‑70 (MQ=255)
                     agcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGt                                            >  1:2832942/1‑70 (MQ=255)
                      gcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAAt                                                                      <  1:1683641/43‑1 (MQ=255)
                      gcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAAt                                                                      <  1:1465501/43‑1 (MQ=255)
                      gcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATa                                                                     <  1:928088/44‑1 (MQ=255)
                      gcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATa                                                                     <  1:3833950/44‑1 (MQ=255)
                          cagccaAATCTTTTTCATGATCTACTCTCTCCGATTAATa                                                                     <  1:3785720/40‑1 (MQ=255)
                          cagccaAATCTTTTTCATGATCTACTCTCTCCGATTAATa                                                                     <  1:3417672/40‑1 (MQ=255)
                          cagccaAATCTTTTTCATGATCTACTCTCTCCGATTAATa                                                                     <  1:3581211/40‑1 (MQ=255)
                          cagccaAATCTTTTTCATGATCTACTCTCTCCGATTAATa                                                                     <  1:3413106/40‑1 (MQ=255)
                          cagccaAATCTTTTTCATGATCTACTCTCTCCGATTAATa                                                                     <  1:3703468/40‑1 (MQ=255)
                            gccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                   >  1:727820/1‑40 (MQ=255)
                            gccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                   >  1:1592025/1‑40 (MQ=255)
                            gccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                   >  1:3384802/1‑40 (MQ=255)
                             ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATATGTGTTAATTGCa                                                    <  1:2432538/54‑1 (MQ=255)
                             ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                    <  1:2510367/54‑1 (MQ=255)
                             ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                    <  1:2424189/54‑1 (MQ=255)
                             ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGa                                   <  1:1509867/71‑1 (MQ=255)
                              caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                  >  1:3400854/1‑71 (MQ=255)
                              caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                  >  1:2253/1‑71 (MQ=255)
                              caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                  >  1:2595720/1‑71 (MQ=255)
                              caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                  >  1:2560214/1‑71 (MQ=255)
                              caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                  >  1:2555435/1‑71 (MQ=255)
                              caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                  >  1:2296033/1‑71 (MQ=255)
                              caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                  >  1:105668/1‑71 (MQ=255)
                                aaTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                    <  1:783016/51‑1 (MQ=255)
                                  tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                              >  1:473590/1‑71 (MQ=255)
                                  tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                              >  1:1327012/1‑71 (MQ=255)
                                  tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                               >  1:2543942/1‑70 (MQ=255)
                                   cTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                                >  1:107521/1‑36 (MQ=255)
                                   cTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                                >  1:1359769/1‑36 (MQ=255)
                                   cTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGa                                   >  1:3080129/1‑65 (MQ=255)
                                   cTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                              >  1:3808099/1‑70 (MQ=255)
                                   cTTTTTCATGATCTACTCTCTCCGATTAACAGCTATGTATTAAtcg                                                      >  1:109306/1‑44 (MQ=16)
                                   cTTTTTCATGATCTACTCTCTCCGATTAACACCTATGTATTAAtcg                                                      >  1:153995/1‑44 (MQ=18)
                                      tttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGcc                          >  1:509583/1‑71 (MQ=255)
                                          aTGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGt                       <  1:4061963/70‑1 (MQ=255)
                                             aTCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                       <  1:1086580/51‑1 (MQ=255)
                                             aTCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                       <  1:3447282/51‑1 (MQ=255)
                                             aTCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                       <  1:1648194/51‑1 (MQ=255)
                                                  ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAg                                    <  1:3367610/49‑1 (MQ=255)
                                                  ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAg                                    <  1:3719996/49‑1 (MQ=255)
                                                  ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAg                                    <  1:3074027/49‑1 (MQ=255)
                                                  ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATaaaa              <  1:1339758/71‑1 (MQ=255)
                                                    ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATaaaa              >  1:2529207/1‑69 (MQ=255)
                                                    ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt            >  1:331933/1‑71 (MQ=255)
                                                    ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt            >  1:3856149/1‑71 (MQ=255)
                                                    ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt            >  1:379008/1‑71 (MQ=255)
                                                    ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt            >  1:3035723/1‑71 (MQ=255)
                                                    ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCCATAAAAGt            >  1:3105047/1‑71 (MQ=255)
                                                    ctctcCGATTAATACATAGGTGTTAATAGCAAAGGGGGTTCTTGTAGAACTTTTGcc                          >  1:1722492/1‑57 (MQ=255)
                                                     tctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAg             <  1:3209840/69‑1 (MQ=255)
                                                     tctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAg             <  1:65672/69‑1 (MQ=255)
                                                     tctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAg             <  1:2705431/69‑1 (MQ=255)
                                                      ctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtttt                             <  1:2147930/52‑1 (MQ=255)
                                                      ctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtttt                             <  1:1800615/52‑1 (MQ=255)
                                                       tcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtttt                             <  1:3861975/51‑1 (MQ=255)
                                                         cGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCg       >  1:601610/1‑71 (MQ=255)
                                                           aTTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                              >  1:270491/1‑46 (MQ=255)
                                                             tAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTg   >  1:889286/1‑71 (MQ=255)
                                                             tAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTg   >  1:4099036/1‑71 (MQ=255)
                                                               aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACt                                <  1:3280183/40‑1 (MQ=255)
                                                               aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACt                                <  1:3740095/40‑1 (MQ=255)
                                                               aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACt                                <  1:3935261/40‑1 (MQ=255)
                                                               aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTGt  <  1:3990424/70‑1 (MQ=255)
                                                               aTACATAGGTGTTAATTGCAAAGGGGGTTCTTCTAGAACt                                <  1:1645190/40‑1 (MQ=255)
                                                                |                                                                    
AACGCTAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTGT  >  minE/2355152‑2355284

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: ‑