Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 976,718 G→A 34.6% 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,7180GA34.6% 26.5 / 21.5 26E216E (GAG→GAAmukEprotein involved in chromosome partitioning
Reads supporting (aligned to +/- strand):  ref base G (6/11);  new base A (5/4);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

ATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGA  >  W3110S.gb/976663‑976784
                                                       |                                                                  
aTGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGc                                                     <  1:1235173/71‑1 (MQ=255)
     gAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGGATCAGCCAGATAGCGgaga                                                <  1:22063/71‑1 (MQ=255)
                  ttGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGa                                   <  1:1034670/71‑1 (MQ=255)
                  ttGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGa                                   >  1:311949/1‑71 (MQ=255)
                       aaTCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                              <  1:1345507/71‑1 (MQ=255)
                        atcatcTGCAACTCAACGATGAAACCGAAGAAAATCAGCCa                                                           <  1:655287/41‑1 (MQ=255)
                           atcTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGgaga                                                <  1:1251575/49‑1 (MQ=255)
                            tcTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCg                                                    <  1:322288/44‑1 (MQ=255)
                               gCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGAtt                                 >  1:1627881/1‑60 (MQ=255)
                                cAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGt                        <  1:575962/68‑1 (MQ=255)
                                 aaCTCAACGTTGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAt                    >  1:1149577/1‑71 (MQ=255)
                                 aaCTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAt                    <  1:350762/71‑1 (MQ=255)
                                  aCTCAACGAGGAAACCGCAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGACCGCGGTAAAtt                   <  1:1307492/71‑1 (MQ=255)
                                     cAACGATGAAACCGAAGAAAATCAGCCAGATAGCg                                                    >  1:1679509/1‑35 (MQ=255)
                                     cAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCg                >  1:1013676/1‑71 (MQ=255)
                                     cAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCg                >  1:1718019/1‑71 (MQ=255)
                                     cAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCg                >  1:497976/1‑71 (MQ=255)
                                      aaCGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTc                 <  1:657875/69‑1 (MQ=255)
                                        cGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTc             <  1:921989/71‑1 (MQ=255)
                                         gATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCa            <  1:366728/71‑1 (MQ=255)
                                         gATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCa            <  1:830969/71‑1 (MQ=255)
                                           tGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACt          >  1:458258/1‑71 (MQ=255)
                                            gAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCg                >  1:1748968/1‑64 (MQ=255)
                                            gAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACt          >  1:710379/1‑70 (MQ=255)
                                               aCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAttt                  >  1:1383731/1‑59 (MQ=255)
                                                   aaGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGAGGTAAATTTCGCTCACTGACGCTGa  <  1:600839/71‑1 (MQ=255)
                                                     gagaATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa                     >  1:1605557/1‑50 (MQ=255)
                                                       |                                                                  
ATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGA  >  W3110S.gb/976663‑976784

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