Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,649,687 T→C 19.3% 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,6870TC19.3% 59.8 / 8.7 29Y14H (TAT→CAT) dicApredicted regulator for DicB
Reads supporting (aligned to +/- strand):  ref base T (13/10);  new base C (2/3);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTG  >  W3110S.gb/1649621‑1649756
                                                                  |                                                                     
gTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTAtc                                                                    <  1:553749/70‑1 (MQ=255)
       cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaa                                                            <  1:617289/71‑1 (MQ=255)
       cGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGaaa                                                            <  1:495112/71‑1 (MQ=255)
        gATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTTTCGGCGAACGCATCAGGTATCGTCGGaaaa                                                           <  1:299118/71‑1 (MQ=255)
                aTTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCaaa                                                   <  1:877731/71‑1 (MQ=255)
                 ttAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCa                                                     <  1:1853836/68‑1 (MQ=255)
                   aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAacac                                                >  1:2290456/1‑71 (MQ=255)
                   aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                                 >  1:1046486/1‑70 (MQ=255)
                   aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                                 >  1:434914/1‑70 (MQ=255)
                   aGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAaca                                                 >  1:2008907/1‑70 (MQ=255)
                       aaCTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTAt                                                                     >  1:299413/1‑46 (MQ=255)
                           aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                            <  1:1859500/51‑1 (MQ=255)
                           aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGaaa                                                            <  1:2271234/51‑1 (MQ=255)
                           aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGCATCGTCGGAAAAACCTCAAaca                                                 >  1:1194018/1‑62 (MQ=255)
                           aTGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGAATCGTCGGaaa                                                            <  1:935178/51‑1 (MQ=255)
                            tGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGt                                                                       <  1:1538053/39‑1 (MQ=255)
                            tGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                        <  1:537341/70‑1 (MQ=255)
                             ggAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAacacac                                              >  1:2167808/1‑63 (MQ=255)
                             ggAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAacacac                                              >  1:1570518/1‑63 (MQ=255)
                                         tttAACTATCGGCGAACGCATCAGGCATCGTCGGaaaa                                                           <  1:477983/38‑1 (MQ=255)
                                                 tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGaa                   <  1:395883/70‑1 (MQ=255)
                                                   ggCGAACGCATCAGGCATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAAt                >  1:2138557/1‑71 (MQ=255)
                                                        aCGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTcc            <  1:148720/70‑1 (MQ=255)
                                                               aGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                        >  1:1717513/1‑35 (MQ=255)
                                                               aGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                                        >  1:311300/1‑35 (MQ=255)
                                                               aGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGtctc                                   >  1:1531909/1‑40 (MQ=255)
                                                               aGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTctct                                  >  1:358078/1‑41 (MQ=255)
                                                               aGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTctct                                  >  1:1552347/1‑41 (MQ=255)
                                                               aGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTc    >  1:2156724/1‑71 (MQ=255)
                                                                 gTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTg  <  1:1207858/71‑1 (MQ=255)
                                                                  |                                                                     
GTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTG  >  W3110S.gb/1649621‑1649756

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