Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,973,768 T→C 28.6% I23I (ATT→ATC metK → methionine adenosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,973,7680TC28.6% 146.6 / 65.5 143I23I (ATT→ATCmetKmethionine adenosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (55/47);  new base C (18/23);  total (73/70)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 150 total aligned reads displayed.

ATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTACG  >  minE/1973700‑1973835
                                                                    |                                                                   
atGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAAtt                                                                     >  1:588236/1‑69 (MQ=255)
   gCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCt                                                                  <  1:699678/69‑1 (MQ=255)
   gCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCt                                                                  <  1:3204500/69‑1 (MQ=255)
           ccTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:2355623/1‑69 (MQ=255)
           catttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:70289/3‑69 (MQ=255)
             tttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:2266304/1‑67 (MQ=255)
             tttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:1815616/1‑67 (MQ=255)
             tttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:3327634/1‑67 (MQ=255)
             tttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:2510381/1‑67 (MQ=255)
             tttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:642115/1‑67 (MQ=255)
             tttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:2893627/1‑67 (MQ=255)
             tttttACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt                                                          >  1:1223672/1‑67 (MQ=255)
                 tACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTa                                                      <  1:1844241/67‑1 (MQ=255)
                 tACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTa                                                      <  1:2281351/67‑1 (MQ=255)
                 tACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTa                                                      <  1:229833/67‑1 (MQ=255)
                  aCGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTa                                                      <  1:925903/66‑1 (MQ=255)
                      ccGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGa                                               <  1:2284340/69‑1 (MQ=255)
                      ccGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGa                                               <  1:2442029/69‑1 (MQ=255)
                           tCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAAtt                                                                     >  1:1824732/1‑42 (MQ=255)
                           tCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCt                                           >  1:1916566/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:95750/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:2930873/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:1095228/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:876137/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:1987504/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:2074468/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:278101/1‑68 (MQ=255)
                             cGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCg                                         >  1:3148810/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:3225088/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:3361734/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:1024408/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:2686667/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:2823263/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:3000177/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:38807/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:1747159/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:1362765/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:2212255/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:470250/1‑68 (MQ=255)
                                ctctGAAGGGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:1618708/1‑68 (MQ=255)
                                ctctGAAGCGCATCCTGACAAAATTGCTGACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAAc                                      >  1:2867260/1‑68 (MQ=255)
                                       gggCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATcc                               >  1:2627890/1‑68 (MQ=255)
                                       gggCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATcc                               >  1:3005971/1‑68 (MQ=255)
                                       gggCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATcc                               >  1:2849832/1‑68 (MQ=255)
                                       gggCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATcc                               >  1:1454779/1‑68 (MQ=255)
                                          cATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGaa                            >  1:1403816/1‑68 (MQ=255)
                                          cATCCTGACAAAATTCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAg                          <  1:3062740/69‑1 (MQ=255)
                                           aTCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAg                          <  1:1483942/69‑1 (MQ=255)
                                           aTCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAg                          <  1:2813108/69‑1 (MQ=255)
                                             ccTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCa                        <  1:2421599/69‑1 (MQ=255)
                                             ccTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCa                        <  1:2522463/69‑1 (MQ=255)
                                             ccTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCa                        <  1:3220731/69‑1 (MQ=255)
                                             ccTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCa                        <  1:1622244/69‑1 (MQ=255)
                                             ccTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCa                        <  1:1227847/69‑1 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCa                        <  1:2931699/68‑1 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:2134227/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:303619/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:29371/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:980592/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:1736184/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:2300289/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:908518/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:1310350/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:436680/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:2722737/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:3185096/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:137876/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:2358381/1‑69 (MQ=255)
                                              cTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCAc                       >  1:259157/1‑69 (MQ=255)
                                                            gACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:1039355/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:2884534/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:1050898/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:1198649/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:1228223/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:923532/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:125719/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:129852/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:14410/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:1511076/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:164042/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:2095563/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:3429920/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:2422131/69‑1 (MQ=255)
                                                            gACCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:2472261/69‑1 (MQ=255)
                                                             aCCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:901139/68‑1 (MQ=255)
                                                             aCCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:1419366/68‑1 (MQ=255)
                                                             aCCAAATCTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGaa         <  1:3286913/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAcc      <  1:2938150/69‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:638202/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:1536613/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:251356/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:2526594/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:467010/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:2710996/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:2789506/68‑1 (MQ=255)
                                                               cAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAAc       <  1:3384725/68‑1 (MQ=255)
                                                                 aaTTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTa    <  1:839291/69‑1 (MQ=255)
                                                                 aaTTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTa    <  1:2583916/69‑1 (MQ=255)
                                                                    ttCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTACg  >  1:616074/1‑68 (MQ=255)
                                                                    ttCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTACg  >  1:125457/1‑68 (MQ=255)
                                                                    |                                                                   
ATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTAGACGCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTACG  >  minE/1973700‑1973835

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: