Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 982,670 A→C 16.2% N39T (AAT→ACT)  ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE982,6700AC16.2% 151.4 / 27.1 74N39T (AAT→ACT) ydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/37);  new base C (4/8);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAA  >  minE/982607‑982736
                                                               |                                                                  
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                             >  1:4718072/1‑71 (MQ=255)
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                             >  1:578358/1‑71 (MQ=255)
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                             >  1:5012760/1‑71 (MQ=255)
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc                                                              <  1:1170512/70‑1 (MQ=255)
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc                                                              <  1:1448562/70‑1 (MQ=255)
  aaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcg                                                           <  1:4637559/71‑1 (MQ=255)
  aaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcg                                                           <  1:5188924/71‑1 (MQ=255)
  aaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCg                                                           >  1:4722547/1‑71 (MQ=255)
      aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga                                                       <  1:3936469/71‑2 (MQ=255)
      aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga                                                       <  1:3247907/71‑2 (MQ=255)
         aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg                                                    <  1:3627447/71‑1 (MQ=255)
          cGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                    >  1:4481924/1‑54 (MQ=255)
          cGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATg                                                                  >  1:670407/1‑56 (MQ=255)
             cACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                                                   <  1:724862/52‑1 (MQ=255)
             aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGatggttgg                                                <  1:3552788/70‑1 (MQ=255)
               cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                                 >  1:5157372/1‑68 (MQ=255)
               cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                                 >  1:3223200/1‑68 (MQ=255)
               cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAAATGGTGCGCGCCgatggttg                                                 >  1:3871693/1‑68 (MQ=255)
               cGTTCACCGCAATGACGTATACCAATAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                             <  1:8980/56‑1 (MQ=255)
                   cACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                              >  1:5012133/1‑66 (MQ=255)
                   cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                           <  1:4180364/70‑1 (MQ=255)
                   cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                           <  1:2573547/70‑1 (MQ=255)
                   cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                           <  1:4847846/70‑1 (MQ=255)
                      ctcAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACt                                       >  1:2670830/3‑71 (MQ=255)
                      cGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTg                                                               <  1:3073809/47‑1 (MQ=255)
                        ccaTGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTg                                                                  >  1:5532609/3‑42 (MQ=37)
                          aTGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                             <  1:1546652/61‑1 (MQ=255)
                          aTGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                             <  1:5006571/61‑1 (MQ=255)
                           tGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                                                   <  1:2734298/38‑1 (MQ=255)
                           tGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                                                   <  1:5265494/38‑1 (MQ=255)
                           tGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                  <  1:2549657/71‑1 (MQ=255)
                           tGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                  <  1:3351675/71‑1 (MQ=255)
                            gACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                    >  1:2678744/1‑36 (MQ=255)
                            gACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                    >  1:5386829/1‑36 (MQ=255)
                             aCGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                    >  1:811884/1‑35 (MQ=255)
                             aCGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                    >  1:4665001/1‑35 (MQ=255)
                             aCGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                                                                 >  1:2599585/1‑38 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                         <  1:5663943/41‑1 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                         <  1:5434580/41‑1 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                         <  1:3813525/41‑1 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAg                                    <  1:1387796/62‑1 (MQ=255)
                                    ccAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgagtga                            <  1:3495523/68‑1 (MQ=255)
                                    ccAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGatat                         >  1:2952438/1‑71 (MQ=255)
                                      aaaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATg                        >  1:1078776/1‑70 (MQ=255)
                                       aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa                      >  1:3235177/1‑71 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa                                              <  1:1216691/46‑1 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa                                              <  1:4209019/46‑1 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAt                     >  1:4786817/1‑71 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAt                     >  1:3058826/1‑71 (MQ=255)
                                          gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                               >  1:4030641/1‑59 (MQ=255)
                                          gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                               >  1:824947/1‑59 (MQ=255)
                                          gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                               >  1:2541855/1‑59 (MQ=255)
                                          gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                    <  1:5535686/38‑1 (MQ=16)
                                          gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                    <  1:138582/38‑1 (MQ=16)
                                          gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                    <  1:904566/38‑1 (MQ=16)
                                           cAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTc                  >  1:3982701/1‑71 (MQ=255)
                                            aaaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGatat                         <  1:1819954/63‑1 (MQ=255)
                                            aaaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGatat                         <  1:1281384/63‑1 (MQ=255)
                                             aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTc                                      <  1:4594842/49‑1 (MQ=255)
                                             aaaCGTTTAAAACCGAAAATGGTGCCCGGCgctggctg                                                 <  1:5362432/38‑1 (MQ=255)
                                              aaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                                         >  1:1557849/1‑45 (MQ=255)
                                              aaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                           <  1:4373710/43‑1 (MQ=255)
                                               aCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                               >  1:848033/1‑38 (MQ=255)
                                                cGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTc                  <  1:2853523/66‑1 (MQ=255)
                                                 gTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGacac             >  1:1360133/1‑70 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGTCTGga                                               <  1:1823361/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                               <  1:3863878/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                               <  1:52607/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                               <  1:4712045/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                               <  1:462759/35‑1 (MQ=255)
                                                   ttAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAAt          >  1:231596/1‑71 (MQ=255)
                                                    tAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                                         <  1:2010744/39‑1 (MQ=255)
                                                     aaaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAgg                                   >  1:2804462/1‑44 (MQ=255)
                                                      aaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCa        <  1:2700682/70‑1 (MQ=255)
                                                        aCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAgg                                   >  1:4197133/1‑41 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTc                                      >  1:3943956/1‑35 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAg                                    >  1:3515714/1‑37 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGaaa  <  1:1496753/71‑1 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGa    <  1:4864487/69‑1 (MQ=255)
                                                            aaaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                  <  1:3779239/38‑1 (MQ=255)
                                                            aaaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                  <  1:3168484/38‑1 (MQ=255)
                                                              aaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGata                          >  1:2644295/1‑44 (MQ=255)
                                                               aTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                               <  1:2108554/38‑1 (MQ=255)
                                                               aTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                               <  1:4265683/38‑1 (MQ=255)
                                                               |                                                                  
CAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAA  >  minE/982607‑982736

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: