Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,649,687 T→C 43.1% 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,6870TC43.1% 17.0 / 65.1 51Y14H (TAT→CAT) dicApredicted regulator for DicB
Reads supporting (aligned to +/- strand):  ref base T (24/5);  new base C (10/12);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 7.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

AATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTGT  >  W3110S.gb/1649624‑1649757
                                                               |                                                                      
aaTACGATTTAGCATTAGCTAACTATTGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCAtcgtcg                                                                 >  1:33511/1‑71 (MQ=255)
aaTACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCAtcgtcg                                                                 >  1:415883/1‑71 (MQ=255)
   aCGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGa                                                               >  1:175315/1‑70 (MQ=255)
    cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaa                                                              >  1:1550538/1‑70 (MQ=255)
    cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaa                                                              >  1:3487099/1‑70 (MQ=255)
    cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaa                                                              >  1:3314963/1‑70 (MQ=255)
    cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:3741383/71‑1 (MQ=255)
        ttAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCAtcgtcg                                                                 >  1:3121468/1‑63 (MQ=255)
        ttAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaaaa                                                           >  1:792491/1‑69 (MQ=255)
        ttAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAAcc                                                         >  1:3156627/1‑71 (MQ=255)
            cATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:1258946/63‑1 (MQ=255)
            cATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:764237/63‑1 (MQ=255)
            cATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:426841/63‑1 (MQ=255)
            cATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:4044536/63‑1 (MQ=255)
            cATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:3334575/63‑1 (MQ=255)
            cATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:2448773/63‑1 (MQ=255)
             aTTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:134797/62‑1 (MQ=255)
             aTTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                             <  1:1698843/62‑1 (MQ=255)
              ttAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaaaa                                                           >  1:1713885/1‑63 (MQ=255)
                aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                                  <  1:2466545/70‑1 (MQ=255)
                aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                                  <  1:130775/70‑1 (MQ=255)
                aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                                  >  1:3961647/1‑70 (MQ=255)
                aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                                  >  1:1927283/1‑70 (MQ=255)
                  cTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaa                                                             <  1:2384957/57‑1 (MQ=255)
                      cTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaa                                                             <  1:3097648/53‑1 (MQ=255)
                       tATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCa                                            >  1:1977363/1‑69 (MQ=255)
                        aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaaa                                                           >  1:657704/1‑53 (MQ=255)
                        aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaaa                                                           >  1:2586262/1‑53 (MQ=255)
                        aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCaa                                           >  1:1397198/1‑69 (MQ=255)
                        aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                         >  1:2193851/1‑71 (MQ=255)
                        aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                         >  1:1071665/1‑71 (MQ=255)
                        aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                         >  1:428281/1‑71 (MQ=255)
                         tGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaaaa                                                           >  1:3848024/1‑52 (MQ=255)
                         tGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAc                                                          >  1:2894645/1‑53 (MQ=255)
                            aaaCAAAAAATTTAACTATCGGCGAACGAATCAGGTAtcg                                                                    >  1:3918376/1‑40 (MQ=255)
                                 aaaaaTTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGc                                >  1:3437249/1‑71 (MQ=255)
                                 aaaaaTTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGc                                >  1:3441823/1‑71 (MQ=255)
                                      tttAACTATCGGCGAACGCATCAGGTATCGTCGGaaaa                                                            >  1:4071210/1‑38 (MQ=255)
                                      tttAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACCCCCaaa                                          >  1:2132076/1‑56 (MQ=255)
                                      tttAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCaaa                                          >  1:1515582/1‑56 (MQ=255)
                                      tttAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCaaa                                          >  1:3163289/1‑56 (MQ=255)
                                      tttAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAg                           >  1:2177832/1‑71 (MQ=255)
                                              tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGaa                    <  1:1284847/70‑1 (MQ=255)
                                              tCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGaaa                   <  1:205671/71‑1 (MQ=255)
                                              tCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGaaa                   <  1:729072/71‑1 (MQ=255)
                                              tCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGaaa                   <  1:752801/71‑1 (MQ=255)
                                                ggCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGaaaa                  >  1:58674/1‑70 (MQ=255)
                                                      cGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTcc             >  1:2962591/1‑69 (MQ=255)
                                                            aGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTc     >  1:1695353/1‑71 (MQ=255)
                                                            aGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTc     >  1:2305254/1‑71 (MQ=255)
                                                            aGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTc     >  1:3825474/1‑71 (MQ=255)
                                                               tATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTGt  <  1:2855152/71‑1 (MQ=255)
                                                               |                                                                      
AATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTGT  >  W3110S.gb/1649624‑1649757

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: