Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,649,687 T→C 16.8% Y14H (TAT→CAT)  dicA → predicted regulator for DicB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,649,6870TC16.8% 117.1 / 16.2 51Y14H (TAT→CAT) dicApredicted regulator for DicB
Reads supporting (aligned to +/- strand):  ref base T (18/24);  new base C (4/4);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTG  >  W3110S.gb/1649626‑1649752
                                                             |                                                                 
tACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGa                                                          <  1:755450/71‑1 (MQ=255)
tACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGa                                                          <  1:2521253/71‑1 (MQ=255)
  cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaa                                                        <  1:1617909/71‑1 (MQ=255)
  cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaa                                                         >  1:3532584/1‑70 (MQ=255)
  cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaa                                                         >  1:2925816/1‑70 (MQ=255)
  cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaa                                                         >  1:2060152/1‑70 (MQ=255)
  cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaa                                                         >  1:1967933/1‑70 (MQ=255)
   gATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaa                                                       <  1:903733/71‑1 (MQ=255)
   gATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaa                                                       <  1:2989158/71‑1 (MQ=255)
   gATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaa                                                       <  1:2828864/71‑1 (MQ=255)
   gATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGAATCGTCGGaaaa                                                       <  1:4578152/71‑1 (MQ=255)
      ttAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAc                                                     >  1:2114444/1‑70 (MQ=255)
      ttAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAc                                                     >  1:1370379/1‑70 (MQ=255)
         gCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCa                                                 <  1:4108874/71‑1 (MQ=255)
           aTTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCaaa                                               <  1:1676512/71‑1 (MQ=255)
           aTTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCaaa                                               <  1:412318/71‑1 (MQ=255)
            ttAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGt                                                                   <  1:1381423/50‑1 (MQ=255)
            ttAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAac                                              <  1:1830106/71‑1 (MQ=255)
              aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                             >  1:3083074/1‑70 (MQ=255)
              aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                             >  1:3323972/1‑70 (MQ=255)
              aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAaca                                             <  1:4117716/70‑1 (MQ=255)
                  aaCTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACAcc                                         >  1:2202627/1‑70 (MQ=255)
                      aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTAtc                                                                <  1:2215617/43‑1 (MQ=255)
                      aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTAtc                                                                <  1:3179795/43‑1 (MQ=255)
                      aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCaa                                      <  1:3453017/69‑1 (MQ=255)
                      aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                    >  1:1482540/1‑71 (MQ=255)
                         gAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaaa                                                      >  1:802521/1‑50 (MQ=255)
                         gAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaaa                                                      >  1:4022156/1‑50 (MQ=255)
                         gAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGt                                  >  1:3114016/1‑70 (MQ=255)
                         gAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGtc                                 <  1:1256996/71‑1 (MQ=255)
                          aaaCAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGtct                                >  1:3702460/1‑71 (MQ=255)
                               aaaaaTTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTg                            <  1:521302/70‑1 (MQ=255)
                                 aaaTTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTa                         <  1:4593559/71‑1 (MQ=255)
                                   aTTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                    >  1:1840639/1‑58 (MQ=255)
                                   aTTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTaaa                       <  1:527657/71‑1 (MQ=255)
                                    tttAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAg                      <  1:3335671/71‑1 (MQ=255)
                                     ttAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAcc                                                    <  1:2006756/40‑1 (MQ=255)
                                      tAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCa                                       <  1:1208541/52‑1 (MQ=255)
                                        aCTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGccc                   <  1:316686/70‑1 (MQ=255)
                                        aCTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGcc                    <  1:280839/69‑1 (MQ=255)
                                         cTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAac                                              >  1:4622858/1‑42 (MQ=255)
                                         cTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAac                                              >  1:2521936/1‑42 (MQ=255)
                                           aTCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCt                          <  1:2031141/60‑1 (MQ=255)
                                            tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                             >  1:1524920/1‑40 (MQ=255)
                                            tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGa                <  1:989484/69‑1 (MQ=255)
                                               gCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAAt            <  1:1791905/70‑1 (MQ=255)
                                                    cGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTaaa                       >  1:683068/1‑54 (MQ=255)
                                                      cATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCAtg    >  1:3131828/1‑71 (MQ=255)
                                                      cATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCAtg    >  1:1197924/1‑71 (MQ=255)
                                                       aTCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCAtgt   >  1:4223354/1‑71 (MQ=255)
                                                        tCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCAtgtg  <  1:2804431/71‑1 (MQ=255)
                                                        tCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCAtgtg  <  1:1089198/71‑1 (MQ=255)
                                                             |                                                                 
TACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTG  >  W3110S.gb/1649626‑1649752

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: