Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 976,718 G→A 17.5% E216E (GAG→GAA mukE → protein involved in chromosome partitioning

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb976,7180GA17.5% 138.7 / 20.2 57E216E (GAG→GAAmukEprotein involved in chromosome partitioning
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base A (7/3);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCT  >  W3110S.gb/976651‑976782
                                                                    |                                                                
gCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAAt                                                               >  1:3116275/1‑71 (MQ=255)
   tGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAg                                                            >  1:4104930/1‑71 (MQ=255)
   tGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAg                                                            >  1:3967013/1‑71 (MQ=255)
    gATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAagag                                                                  <  1:348773/64‑1 (MQ=255)
    gATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGc                                                           >  1:3317753/1‑71 (MQ=255)
      ttCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAagag                                                                  <  1:989084/62‑1 (MQ=255)
      ttCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAagag                                                                  <  1:3832826/62‑1 (MQ=255)
      ttCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGc                                                           >  1:3400229/1‑69 (MQ=255)
          tGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGAt                                                      <  1:100262/71‑1 (MQ=255)
          tGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGAt                                                      <  1:2453496/71‑1 (MQ=255)
           caTGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAg                                                    <  1:3968258/70‑1 (MQ=255)
            aTGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAg                                                    <  1:2841881/70‑1 (MQ=255)
              ggCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCg                                                  >  1:5031962/1‑70 (MQ=255)
               gCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGa                                                       >  1:2837346/1‑64 (MQ=255)
                   aGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGgaga                                              >  1:4240271/1‑69 (MQ=255)
                      aaTGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGAt                                                      <  1:3700448/59‑1 (MQ=255)
                      aaTGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaa                                          <  1:4193299/71‑1 (MQ=255)
                      aaTGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaa                                          <  1:1476287/70‑1 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCg                                                  >  1:471639/1‑60 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaagaa                                       >  1:3550352/1‑71 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGaaaa                                                                >  1:309676/1‑45 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaagaa                                       >  1:2782022/1‑70 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaagaa                                       >  1:3033077/1‑70 (MQ=255)
                           cGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGaata                                     <  1:3325893/71‑1 (MQ=255)
                           cGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGaata                                     <  1:1654668/71‑1 (MQ=255)
                            gATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGaataat                                   <  1:4525721/71‑1 (MQ=255)
                              ttGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGa                                 >  1:2933224/1‑71 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATGAAACCGAagag                                                                  >  1:3799620/1‑37 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATGAAACCGAagag                                                                  >  1:1023370/1‑37 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATGAAACCGAagag                                                                  >  1:4704804/1‑37 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGAt                                >  1:2874946/1‑71 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGAt                                >  1:2482775/1‑71 (MQ=255)
                                   aaTCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                            >  1:860734/1‑71 (MQ=255)
                                   aaTCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGa                             <  1:1885533/70‑1 (MQ=255)
                                    atcatcTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAAc                           >  1:4785491/1‑71 (MQ=255)
                                    atcatcTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAAc                           >  1:2687563/1‑71 (MQ=255)
                                        tcTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCg                                                  <  1:2919911/44‑1 (MQ=255)
                                         cTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCg                                                  <  1:4372492/43‑1 (MQ=255)
                                          tGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCgg                                                 >  1:4166350/1‑43 (MQ=255)
                                          tGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaaga                                        >  1:157071/1‑52 (MQ=255)
                                            cAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:2199642/1‑71 (MQ=255)
                                            cAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:3820526/1‑71 (MQ=255)
                                            cAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:1176026/1‑71 (MQ=255)
                                              aCTCAGCGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAtt                 >  1:2752152/1‑71 (MQ=255)
                                                tCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTg                              <  1:3507937/57‑1 (MQ=255)
                                                 cAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:625100/1‑66 (MQ=255)
                                                 cAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAgg                                            >  1:2612029/1‑41 (MQ=255)
                                                 cAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGaataa                                    >  1:4332486/1‑49 (MQ=255)
                                                 cAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGaataa                                    >  1:3226657/1‑49 (MQ=255)
                                                  aaCGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                            >  1:4449011/1‑56 (MQ=255)
                                                  aaCGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                            >  1:1313782/1‑56 (MQ=255)
                                                  aaCGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaa                    <  1:782459/64‑1 (MQ=255)
                                                   aCGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCt            <  1:2020295/71‑1 (MQ=255)
                                                   aCGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCt            <  1:2700204/71‑1 (MQ=255)
                                                   aCGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCt            <  1:5002453/71‑1 (MQ=255)
                                                       tGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACt        <  1:1074223/71‑1 (MQ=255)
                                                          aaaCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGa      >  1:947797/1‑71 (MQ=255)
                                                          aaaCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGa      >  1:3099121/1‑71 (MQ=255)
                                                            aCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGaata                                     >  1:2382163/1‑38 (MQ=255)
                                                             ccGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTa                     <  1:5088035/53‑1 (MQ=255)
                                                             ccGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGc   <  1:3020548/71‑1 (MQ=255)
                                                              cGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGa                                 >  1:2394979/1‑40 (MQ=255)
                                                              cGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCg              <  1:495869/59‑1 (MQ=255)
                                                              cGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCt  <  1:2527354/71‑1 (MQ=255)
                                                                   agaATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCt  <  1:403884/66‑1 (MQ=255)
                                                                    gaATCAGCCAGATAGCGGAGAGGAAGAATAATGATTg                              <  1:2945192/37‑1 (MQ=255)
                                                                    |                                                                
GCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCT  >  W3110S.gb/976651‑976782

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: