Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,074,351 T→C 39.1% intergenic (+68/‑40) valW → / → ydhR tRNA‑Val/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,074,3510TC39.1% 41.8 / 62.2 92intergenic (+68/‑40)valW/ydhRtRNA‑Val/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base C (14/22);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATTTTCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCATT  >  minE/1074283‑1074415
                                                                    |                                                                
aGATTTTCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGAtt                                                                  <  1:885730/69‑1 (MQ=255)
aGATTTTCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGAtt                                                                  <  1:1780195/69‑1 (MQ=255)
aGATTTTCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGAtt                                                                  <  1:3328095/69‑1 (MQ=255)
aGATTTTCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGAtt                                                                  <  1:3110548/69‑1 (MQ=255)
   ttttCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCTATACGTTAAAAGATTTac                                                               <  1:1829752/69‑1 (MQ=255)
   ttttCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTac                                                               <  1:2593520/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:2576866/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:1962460/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:3014092/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:3084795/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:3137702/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:3154931/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:3177151/69‑1 (MQ=255)
             cTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:1742786/69‑1 (MQ=255)
              tGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:2574746/68‑1 (MQ=255)
                  cttctCCTTTTTCCCTCTGGTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTa                                                   <  1:1625415/66‑1 (MQ=255)
                    tctcCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCg                                                     <  1:2071892/62‑1 (MQ=255)
                     ctcCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGtt                                             >  1:2506192/1‑69 (MQ=255)
                     ctcCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGtt                                             >  1:477130/1‑69 (MQ=255)
                     ctcCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGtt                                             >  1:577582/1‑69 (MQ=255)
                     ctcCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGtt                                             >  1:2018348/1‑69 (MQ=255)
                     ctcCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGt                                              >  1:2329736/1‑68 (MQ=255)
                     ctcCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACCTTAAAAGATTTACACTGTc                                                         >  1:1474127/1‑57 (MQ=255)
                                  ctGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGgag                                  >  1:3069455/1‑67 (MQ=255)
                                  ctGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGaa                                >  1:2885141/1‑69 (MQ=255)
                                  ctGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGaa                                >  1:2530917/1‑69 (MQ=255)
                                      ttcttcTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAgg                                    >  1:1034227/1‑61 (MQ=255)
                                      ttcttcTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAgg                                    >  1:89487/1‑61 (MQ=255)
                                      atcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAgg                                    >  1:2112470/2‑61 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGgag                                  <  1:2933058/62‑1 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            <  1:971014/68‑1 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            <  1:1019812/68‑1 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            <  1:2283110/68‑1 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            <  1:1506126/68‑1 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            <  1:3026203/68‑1 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            >  1:2990476/1‑68 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            <  1:1715047/68‑1 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCt                           >  1:1988495/1‑69 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCt                           >  1:109248/1‑69 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCt                           >  1:399649/1‑69 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCt                           >  1:2891496/1‑69 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCt                           >  1:2757869/1‑69 (MQ=255)
                                       tcttctCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCt                           >  1:2300358/1‑69 (MQ=255)
                                       tcttctCTGTATACAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGc                            <  1:1173585/68‑1 (MQ=255)
                                             ctGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGc                      <  1:506065/68‑1 (MQ=255)
                                             ctGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGc                      <  1:862783/68‑1 (MQ=255)
                                             ctGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGc                      <  1:2788923/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:2551055/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:3171057/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:3120636/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:2903414/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:486393/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:2812587/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:2508125/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:1935516/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:2674439/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:1650129/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:1692132/68‑1 (MQ=255)
                                                    cAATACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:1789583/68‑1 (MQ=255)
                                                     aaTACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAg                                     >  1:2381739/1‑45 (MQ=255)
                                                     aaTACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGt              >  1:1457635/1‑68 (MQ=255)
                                                     aaTACGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTg               <  1:1282892/67‑1 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGtt                                             >  1:1916749/1‑35 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCt                           >  1:2590787/1‑53 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTa            >  1:1993116/1‑68 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           <  1:2188913/69‑1 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           >  1:404065/1‑69 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           >  1:992689/1‑69 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           >  1:1309782/1‑69 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           <  1:2833016/69‑1 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           >  1:1556656/1‑69 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           >  1:2576172/1‑69 (MQ=255)
                                                       tACGTTAAAAGATTTACACTGTCTTCGTATGCGTGATCAGAAGGAGAATCGCTATGGCAACTTTGTTAc           >  1:1569665/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:385569/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:3271196/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:2505514/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:957189/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:1852590/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:2353264/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:521385/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:7699/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:1059194/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:563023/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:129017/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:686135/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:699785/1‑69 (MQ=255)
                                                        aCGTTAAAAGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACa          >  1:712478/1‑69 (MQ=255)
                                                                aGATCTACTCTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:2718475/69‑1 (MQ=255)
                                                                aGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:928638/69‑1 (MQ=255)
                                                                aGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:2551597/69‑1 (MQ=255)
                                                                aGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:536994/69‑1 (MQ=255)
                                                                aGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:533294/69‑1 (MQ=255)
                                                                aGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:116621/69‑1 (MQ=255)
                                                                aGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:3376339/69‑1 (MQ=255)
                                                                aGATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:1473530/69‑1 (MQ=255)
                                                                 gATCTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCAtt  <  1:269581/68‑1 (MQ=255)
                                                                    |                                                                
AGATTTTCTTAATCTGGTCTTCTCCTTTTTCCCTCTGTTTCTTCTCTGTATCCAATACGTTAAAAGATTTACACTGTCTTCGTATGCGTTATCAGAAGGAGAATCGCTATGGCAACTTTGTTACAACTTCATT  >  minE/1074283‑1074415

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: