Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,750 A→T 41.5% intergenic (‑243/+427) yobF ← / ← yebO hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,211,7500AT41.5% 37.4 / 80.7 82intergenic (‑243/+427)yobF/yebOhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/17);  new base T (12/22);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

AGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGT  >  minE/1211683‑1211814
                                                                   |                                                                 
aGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAg                                                                  <  1:789371/69‑1 (MQ=255)
            tACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCAGGCCAGTGGGAGAAGttat                                                     <  1:1100563/69‑1 (MQ=255)
             aCTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttat                                                     <  1:3179985/69‑1 (MQ=255)
             aCTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttat                                                     <  1:723131/69‑1 (MQ=255)
              cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAAAAGttat                                                     >  1:2321156/1‑68 (MQ=255)
              cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAAAAGttat                                                     >  1:303124/1‑68 (MQ=255)
              cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAAAAGttat                                                     >  1:2714167/1‑68 (MQ=255)
              cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAAAAGttat                                                     >  1:1226399/1‑68 (MQ=255)
              cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAAAAGttat                                                     >  1:484342/1‑68 (MQ=255)
              cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAAAAGttat                                                     >  1:1149396/1‑68 (MQ=255)
              cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAAAAGttat                                                     >  1:739421/1‑68 (MQ=255)
                tAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta                                                   >  1:129461/1‑68 (MQ=255)
                tAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta                                                   >  1:2675031/1‑68 (MQ=255)
                tAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta                                                   >  1:1069975/1‑68 (MQ=255)
                tAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta                                                   >  1:858711/1‑68 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACa                                                 <  1:2184990/67‑1 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACa                                                 <  1:2686456/67‑1 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACa                                                 <  1:3401592/67‑1 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACa                                                 <  1:138169/67‑1 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACa                                                 <  1:578794/67‑1 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATa                                               >  1:3089134/1‑69 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATa                                               >  1:420432/1‑69 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATa                                               >  1:2344348/1‑69 (MQ=255)
                   aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATa                                               >  1:1050952/1‑69 (MQ=255)
                    cTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACa                                                 <  1:634800/66‑1 (MQ=255)
                            ttttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta                                                   >  1:3187705/1‑56 (MQ=255)
                            ttttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta                                                   >  1:3183499/1‑56 (MQ=255)
                            ttttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta                                                   >  1:66663/1‑56 (MQ=255)
                              ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCaa                                     >  1:3140789/1‑68 (MQ=255)
                              ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCaa                                     >  1:1991689/1‑68 (MQ=255)
                              ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                    >  1:621614/1‑69 (MQ=255)
                              ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                    >  1:3229862/1‑69 (MQ=255)
                              ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                    >  1:345977/1‑69 (MQ=255)
                              ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                    >  1:1097373/1‑69 (MQ=255)
                                        gTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtttt                         >  1:1504185/1‑69 (MQ=255)
                                        gTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtttt                         >  1:697338/1‑69 (MQ=255)
                                        gTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtttt                         >  1:2469407/1‑69 (MQ=255)
                                        gTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtttt                         >  1:1079616/1‑69 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATatt                                             >  1:3267992/1‑45 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:2704187/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:2218451/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:329168/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:3410275/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:981109/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:1893826/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:2046657/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:2223069/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:2362422/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:2375024/1‑53 (MQ=255)
                                           cAATAAGTTCAGTTGTTCCAGGCCTGTGGGAGAAGTTATTACATAgtgcgtgc                                       >  1:2465531/1‑53 (MQ=255)
                                              tAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCt                    <  1:1192016/68‑1 (MQ=255)
                                               aaGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCt                    <  1:2040984/67‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2353245/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2695043/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:919057/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2707008/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2327263/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:629234/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2350875/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:1159031/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:395411/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:387035/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:1972384/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:180345/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2676686/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:1776012/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:1873650/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2565252/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2515493/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:3022711/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:3394012/68‑1 (MQ=255)
                                                      gttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:2802547/68‑1 (MQ=255)
                                                      gttgttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATg            <  1:1614977/68‑1 (MQ=255)
                                                       ttgttCCAGGCCTGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:1864936/67‑1 (MQ=255)
                                                       ttgttCCAGGCCTCTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTCTACTATGCAATg            <  1:940108/67‑1 (MQ=255)
                                                       ttgttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTTGCTATGCAATg            <  1:2303536/68‑1 (MQ=255)
                                                            ccAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATg            >  1:372052/1‑62 (MQ=255)
                                                            ccAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGaa          >  1:2240240/1‑64 (MQ=255)
                                                            ccAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATaaaa     >  1:734747/1‑69 (MQ=255)
                                                               ggCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGt  <  1:247401/69‑1 (MQ=255)
                                                               ggCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGt  <  1:1070278/69‑1 (MQ=255)
                                                               ggCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGt  <  1:100767/69‑1 (MQ=255)
                                                                   |                                                                 
AGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGT  >  minE/1211683‑1211814

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: