Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,962,938 C→T 29.6% V140I (GTT→ATT)  yjjI ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,962,9380CT29.6% 114.1 / 49.9 111V140I (GTT→ATT) yjjIconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/50);  new base T (16/17);  total (44/67)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 112 total aligned reads displayed.

GCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCA  >  minE/2962881‑2963000
                                                         |                                                              
gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGt                                                     <  1:1858086/69‑1 (MQ=255)
gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGt                                                     <  1:2054924/69‑1 (MQ=255)
gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGt                                                     <  1:617621/69‑1 (MQ=255)
gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGt                                                     <  1:2634843/69‑1 (MQ=255)
gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGt                                                     <  1:2570027/69‑1 (MQ=255)
gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGt                                                     <  1:2944346/69‑1 (MQ=255)
 cATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGt                                                     <  1:3077389/68‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:1646507/69‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:3365211/69‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:3215588/69‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:382773/69‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:2075993/69‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:2656723/69‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:2803831/69‑1 (MQ=255)
          gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:876942/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:1218571/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:241189/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:1263922/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:710063/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:1106575/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:2660797/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:10826/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:1413689/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:1538879/69‑1 (MQ=255)
           tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:2684492/69‑1 (MQ=255)
           tCTGGCAGGGTGCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:1030378/69‑1 (MQ=255)
            cTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:2838773/68‑1 (MQ=255)
            cTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGa                                          <  1:3082437/68‑1 (MQ=255)
             tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:2012836/66‑1 (MQ=255)
               gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAAc                                                                >  1:149410/1‑43 (MQ=255)
               gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAAtt                                       >  1:3138774/1‑68 (MQ=255)
               gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAAtt                                       >  1:3346051/1‑68 (MQ=255)
               gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAAtt                                       >  1:1351721/1‑68 (MQ=255)
               gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAAtt                                       >  1:129554/1‑68 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:2267403/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:2713673/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:874004/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:806458/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:2937609/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:2299544/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:2179747/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:580365/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:1596447/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:3146151/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTaa                                   >  1:1369427/1‑69 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTa                                    >  1:1254495/1‑68 (MQ=255)
                  gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTa                                    >  1:2020539/1‑68 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGtt                                             <  1:2435188/57‑1 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGtt                                             <  1:794644/57‑1 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGtt                                             <  1:267780/57‑1 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACa                                 >  1:1899696/1‑69 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACa                                 >  1:7027/1‑69 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACa                                 >  1:1138534/1‑69 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACa                                 >  1:2937906/1‑69 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACa                                 >  1:106911/1‑69 (MQ=255)
                    gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACa                                 >  1:2788604/1‑69 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:2809577/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:3318400/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:948483/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:1702911/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:1065270/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:1068277/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:1323599/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:1629527/1‑58 (MQ=255)
                     ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:1800784/1‑58 (MQ=255)
                     atCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           >  1:1591355/2‑58 (MQ=255)
                       cTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAg                                           <  1:3324676/56‑1 (MQ=255)
                            gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGt                           <  1:2802325/67‑1 (MQ=255)
                            gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGt                           <  1:2351996/67‑1 (MQ=255)
                            gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGt                           <  1:1100964/67‑1 (MQ=255)
                            gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGt                           <  1:421414/67‑1 (MQ=255)
                            gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGt                           <  1:611890/67‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTTTTCGAACATCGATCTCGTCTTgtgt                                              <  1:1868292/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2688674/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:141799/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:611559/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2718463/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:746300/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:1087056/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2754559/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:880967/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:90964/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:3302388/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2682086/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2398183/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2343295/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:1853796/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:3121511/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:205734/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2199830/42‑1 (MQ=255)
                                  aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTgtgt                                              <  1:2965608/42‑1 (MQ=255)
                                    gCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACg                 >  1:1110371/1‑69 (MQ=255)
                                    gCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACg                 >  1:42504/1‑69 (MQ=255)
                                           aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAg          >  1:2934124/1‑69 (MQ=255)
                                                ttATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGc      >  1:2109209/1‑68 (MQ=255)
                                                ttATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGc      >  1:135297/1‑68 (MQ=255)
                                                    tCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa  <  1:1109187/68‑1 (MQ=255)
                                                    tCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa  <  1:2702833/68‑1 (MQ=255)
                                                    tCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa  <  1:1159104/68‑1 (MQ=255)
                                                    tCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa  <  1:2521343/68‑1 (MQ=255)
                                                         |                                                              
GCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCA  >  minE/2962881‑2963000

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: