Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,212,559 T→C 8.7% N42D (AAC→GAC)  yobG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,212,5590TC 8.7% 296.3 / 17.6 109N42D (AAC→GAC) yobGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (60/39);  new base C (6/3);  total (66/43)
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.21e-01

Only 100 of 112 total aligned reads displayed.

ATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTGCGCCCAAAG  >  minE/1212490‑1212628
                                                                     |                                                                     
aTTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGtt                                                                       >  1:843831/1‑70 (MQ=255)
   cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                                                   >  1:3754337/1‑71 (MQ=255)
   cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                                                   >  1:2755446/1‑71 (MQ=255)
   cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                                                   >  1:459682/1‑71 (MQ=255)
       cACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGca                                                                >  1:1027310/1‑70 (MQ=255)
       cACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGca                                                                >  1:1238479/1‑70 (MQ=255)
        aCTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGcaca                                                              <  1:1818166/71‑1 (MQ=255)
         ctgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaa                                                             <  1:3040541/71‑1 (MQ=255)
         ctgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaa                                                             <  1:2796192/71‑1 (MQ=255)
         ctgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTCAATg                                                                   <  1:3120063/65‑1 (MQ=255)
            ctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                                                            <  1:2944166/69‑1 (MQ=255)
                agagGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGc                                                      >  1:1716379/1‑71 (MQ=255)
                 gagGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCt                                                     >  1:45164/1‑71 (MQ=255)
                  agGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCt                                                     >  1:3573506/1‑70 (MQ=255)
                    gaagaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGa                                                   >  1:3004142/1‑70 (MQ=255)
                      agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaa                                                 >  1:2262437/1‑70 (MQ=255)
                      agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaa                                                 >  1:1100052/1‑70 (MQ=255)
                       gaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa                                               <  1:624636/71‑1 (MQ=255)
                        aaaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTcc                                                        <  1:3833301/61‑1 (MQ=255)
                        aaaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAAt                                              >  1:2173117/1‑71 (MQ=255)
                        aaaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAAt                                              >  1:2985073/1‑71 (MQ=255)
                        aaaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAAt                                              >  1:3138892/1‑71 (MQ=255)
                        aaaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTCAATGGCACAAATTCCGCTGAAAAAt                                              >  1:2185178/1‑71 (MQ=255)
                         aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAAtt                                             >  1:3802708/1‑71 (MQ=255)
                         aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAAtt                                             >  1:2480224/1‑71 (MQ=255)
                          aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGGACAAATTCCGCTGAAAAATTg                                            <  1:3396858/71‑1 (MQ=255)
                          aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                                                            >  1:3835399/1‑55 (MQ=255)
                          aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg                                            >  1:1993000/1‑71 (MQ=255)
                          aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg                                            <  1:1448947/71‑1 (MQ=255)
                          aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg                                            >  1:255084/1‑71 (MQ=255)
                          aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg                                            >  1:1125906/1‑71 (MQ=255)
                            tCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaa                                                             <  1:1639566/52‑1 (MQ=255)
                             cTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAATATTGTa                                          >  1:3852623/1‑70 (MQ=255)
                             cTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTa                                          >  1:1285447/1‑70 (MQ=255)
                                gTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTa                                          <  1:3149195/67‑1 (MQ=255)
                                gTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATc                                      <  1:1113920/71‑1 (MQ=255)
                                gTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATc                                      <  1:4002793/71‑1 (MQ=255)
                                   cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAGTTCCGCTGAAAAATTGTACATCCTg                                   <  1:2382627/71‑1 (MQ=255)
                                   cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCt                                    <  1:3686029/70‑1 (MQ=255)
                                     gAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa                                               <  1:3250439/57‑1 (MQ=255)
                                      aaaaaaTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATc                                >  1:1683965/1‑71 (MQ=255)
                                      aaaaaaTGATATCACCACGGGATAAACTGGTCAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATc                                >  1:1626892/1‑71 (MQ=255)
                                        aaaaTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTcc                                                        >  1:2299059/1‑45 (MQ=255)
                                        aaaaTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa                                                >  1:1210500/1‑53 (MQ=255)
                                        aaaaTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCg                               >  1:3045170/1‑70 (MQ=255)
                                         aaaTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTa                                          >  1:1620727/1‑58 (MQ=255)
                                           aTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGcaca                           >  1:1729016/1‑71 (MQ=255)
                                            tGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTAc                                         >  1:3634402/1‑56 (MQ=255)
                                            tGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGAt                                 >  1:2576007/1‑64 (MQ=255)
                                            tGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCAcat                          <  1:1428088/71‑1 (MQ=255)
                                             gATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCAcat                          >  1:383431/1‑70 (MQ=255)
                                              atatCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa                                                <  1:3957764/47‑1 (MQ=255)
                                              atatCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa                                                <  1:3627375/47‑1 (MQ=255)
                                              atatCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa                                                <  1:1700022/47‑1 (MQ=255)
                                              atatCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATc                                >  1:661011/1‑63 (MQ=255)
                                              atatCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCAcatca                        <  1:458234/71‑1 (MQ=255)
                                              atatCACCACGGGATAAACTGGTCAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCAcatca                        >  1:2405934/1‑71 (MQ=255)
                                               tatCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCAcatcat                       <  1:3022782/71‑1 (MQ=255)
                                               tatCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCAcatca                        >  1:534531/1‑70 (MQ=255)
                                                 tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa                                               <  1:687856/45‑1 (MQ=255)
                                                 tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa                                               <  1:2030535/45‑1 (MQ=255)
                                                 tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTg                                   <  1:1096776/57‑1 (MQ=255)
                                                 tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATg                      <  1:1302854/70‑1 (MQ=255)
                                                 tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATg                      <  1:2081494/70‑1 (MQ=255)
                                                 tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATg                      <  1:3885303/70‑1 (MQ=255)
                                                 tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGt                     >  1:217198/1‑71 (MQ=255)
                                                 tCACCACGGGATAAACTGGTGAATGGCACAAATTCCGCTGaaa                                                 <  1:3993331/43‑1 (MQ=255)
                                                  caccacGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCAcatca                        >  1:2199140/1‑67 (MQ=255)
                                                   accacGGGATAAACTGGTTAATGGCACAAATTCCGc                                                      >  1:4102236/1‑36 (MQ=255)
                                                   accacGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACa                                        >  1:294280/1‑50 (MQ=255)
                                                   accacGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGca                             >  1:3522441/1‑61 (MQ=255)
                                                   accacGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTg                   <  1:2123548/71‑1 (MQ=255)
                                                   accacGGGATAAACTGGTCAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGtt                    >  1:600615/1‑70 (MQ=255)
                                                    ccacGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTg                                   <  1:20758/54‑1 (MQ=255)
                                                    ccacGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTg                                   <  1:2412364/54‑1 (MQ=255)
                                                    ccacGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTg                                   <  1:3023836/54‑1 (MQ=255)
                                                    ccacGGGATAAACTGGTCAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGc                              >  1:3377738/1‑59 (MQ=255)
                                                    ccacGGGATAAACTGGTCAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGc                              >  1:2756834/1‑59 (MQ=255)
                                                    ccacGGGATAAACTGGTCAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGa                  <  1:539954/71‑1 (MQ=255)
                                                       cGGGATAAACTTGTTAATGGCACAAATTCCGCTGAAAAAt                                              <  1:2446136/40‑1 (MQ=255)
                                                       cGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGAt                                 >  1:632202/1‑53 (MQ=255)
                                                       cGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGAt                                 >  1:3767901/1‑53 (MQ=255)
                                                       cGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGAt                                 >  1:3353607/1‑53 (MQ=255)
                                                       cGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAAt                >  1:1112857/1‑70 (MQ=255)
                                                         ggATAAACTGGTCAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCTCACATCATGTTg                   <  1:1043749/65‑1 (MQ=255)
                                                          gATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTg                                   <  1:4022989/48‑1 (MQ=255)
                                                           aTAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACa                                        >  1:825404/1‑42 (MQ=255)
                                                           aTAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACa                                        >  1:2830599/1‑42 (MQ=255)
                                                            tAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGaa                 >  1:3957207/1‑64 (MQ=255)
                                                             aaaCTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGtt                    >  1:2437937/1‑60 (MQ=255)
                                                             aaaCTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTgcg         >  1:509154/1‑71 (MQ=255)
                                                             aaaCTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTgcg         >  1:2705776/1‑71 (MQ=255)
                                                              aaCTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCg                               <  1:2316177/48‑1 (MQ=255)
                                                              aaCTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTgcgc        >  1:599839/1‑71 (MQ=255)
                                                              aaCTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTgcgc        >  1:1853949/1‑71 (MQ=255)
                                                              aaCTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTgcgc        >  1:261919/1‑71 (MQ=255)
                                                                 tGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCt            >  1:2179651/1‑64 (MQ=255)
                                                                   gTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTGCGCCCaa    <  1:3777772/70‑1 (MQ=255)
                                                                    ttAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGca                             <  1:1571571/44‑1 (MQ=255)
                                                                    ttAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTGCGCCCAAag  <  1:317889/71‑1 (MQ=255)
                                                                     |                                                                     
ATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACATCATGTTGAATACCTGCGCCCAAAG  >  minE/1212490‑1212628

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: