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→GAA) mukEprotein 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

GCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCT  >  W3110S.gb/976651‑976782
                                                                    |                                                                
gCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAAt                                                               >  1:3116275/1‑71 (MQ=255)
   tGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAg                                                            >  1:4104930/1‑71 (MQ=255)
   tGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAg                                                            >  1:3967013/1‑71 (MQ=255)
    gATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAagag                                                                  <  1:348773/64‑1 (MQ=255)
    gATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGc                                                           >  1:3317753/1‑71 (MQ=255)
      ttCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAagag                                                                  <  1:989084/62‑1 (MQ=255)
      ttCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAagag                                                                  <  1:3832826/62‑1 (MQ=255)
      ttCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGc                                                           >  1:3400229/1‑69 (MQ=255)
          tGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGAt                                                      <  1:100262/71‑1 (MQ=255)
          tGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGAt                                                      <  1:2453496/71‑1 (MQ=255)
           caTGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAAAATCAGCCAGATAg                                                    <  1:3968258/70‑1 (MQ=255)
            aTGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAg                                                    <  1:2841881/70‑1 (MQ=255)
              ggCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAAAATCAGCCAGATAGCg                                                  >  1:5031962/1‑70 (MQ=255)
               gCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAAAATCAGCCAGa                                                       >  1:2837346/1‑64 (MQ=255)
                   aGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGgaga                                              >  1:4240271/1‑69 (MQ=255)
                      aaTGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGAt                                                      <  1:3700448/59‑1 (MQ=255)
                      aaTGCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaa                                          <  1:4193299/71‑1 (MQ=255)
                      aaTGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaa                                          <  1:1476287/70‑1 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCg                                                  >  1:471639/1‑60 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaagaa                                       >  1:3550352/1‑71 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGaaaa                                                                >  1:309676/1‑45 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaagaa                                       >  1:2782022/1‑70 (MQ=255)
                         gCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaagaa                                       >  1:3033077/1‑70 (MQ=255)
                           cGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGaata                                     <  1:3325893/71‑1 (MQ=255)
                           cGATTGAAAATCATCTGCAACTCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGaata                                     <  1:1654668/71‑1 (MQ=255)
                            gATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGaataat                                   <  1:4525721/71‑1 (MQ=255)
                              ttGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGa                                 >  1:2933224/1‑71 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATG‑AAACCGAagag                                                                  >  1:3799620/1‑37 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATG‑AAACCGAagag                                                                  >  1:1023370/1‑37 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATG‑AAACCGAagag                                                                  >  1:4704804/1‑37 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGAt                                >  1:2874946/1‑71 (MQ=255)
                               tGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGAt                                >  1:2482775/1‑71 (MQ=255)
                                   aaTCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                            >  1:860734/1‑71 (MQ=255)
                                   aaTCATCTGCAACTCAACGATG‑AAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGa                             <  1:1885533/70‑1 (MQ=255)
                                    atcatcTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAAc                           >  1:4785491/1‑71 (MQ=255)
                                    atcatcTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAAc                           >  1:2687563/1‑71 (MQ=255)
                                        tcTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCg                                                  <  1:2919911/44‑1 (MQ=255)
                                         cTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCg                                                  <  1:4372492/43‑1 (MQ=255)
                                          tGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCgg                                                 >  1:4166350/1‑43 (MQ=255)
                                          tGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGgaaga                                        >  1:157071/1‑52 (MQ=255)
                                            cAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:2199642/1‑71 (MQ=255)
                                            cAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:3820526/1‑71 (MQ=255)
                                            cAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:1176026/1‑71 (MQ=255)
                                              aCTCAGCGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAtt                 >  1:2752152/1‑71 (MQ=255)
                                                tCAACGATGAAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTg                              <  1:3507937/57‑1 (MQ=255)
                                                 cAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                   >  1:625100/1‑66 (MQ=255)
                                                 cAACGATG‑AAACCGAAGAAAATCAGCCAGATAGCGGAGAgg                                            >  1:2612029/1‑41 (MQ=255)
                                                 cAACGATG‑AAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGaataa                                    >  1:4332486/1‑49 (MQ=255)
                                                 cAACGATG‑AAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGaataa                                    >  1:3226657/1‑49 (MQ=255)
                                                  aaCGATG‑AAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                            >  1:4449011/1‑56 (MQ=255)
                                                  aaCGATG‑AAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                            >  1:1313782/1‑56 (MQ=255)
                                                  aaCGATG‑AAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaa                    <  1:782459/64‑1 (MQ=255)
                                                   aCGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCt            <  1:2020295/71‑1 (MQ=255)
                                                   aCGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCt            <  1:2700204/71‑1 (MQ=255)
                                                   aCGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCt            <  1:5002453/71‑1 (MQ=255)
                                                       tG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACt        <  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)
                                                                    |                                                                
GCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATG‑AAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCT  >  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: ‑