Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 887,947 G→A A35T (GCT→ACT)  ybjK → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb887,9470GA97.8% 128.0 / ‑5.7 45A35T (GCT→ACT) ybjKpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/20);  minor base C (0/1);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACT  >  W3110S.gb/887882‑888014
                                                                 |                                                                   
ttATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTAcc                                                                <  1:1212419/71‑1 (MQ=255)
        gCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCgg                                                        >  1:4562765/1‑71 (MQ=255)
          cACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCgggg                                                      >  1:1157175/1‑71 (MQ=255)
          cACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCgggg                                                      >  1:1952529/1‑71 (MQ=255)
                    gCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTggg                                            <  1:4633339/71‑1 (MQ=255)
                    gCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTggg                                            <  1:19703/71‑1 (MQ=255)
                        tGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTAcc                                                                <  1:2478501/47‑1 (MQ=255)
                         gAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGa                                       >  1:819402/1‑71 (MQ=255)
                         gAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGa                                       >  1:1515817/1‑71 (MQ=255)
                         gAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGa                                       >  1:2240984/1‑71 (MQ=255)
                         gAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGa                                       >  1:3947992/1‑71 (MQ=255)
                         gAAACTTTACGGAATACATGCTGTTACGCAACGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGa                                       >  1:573689/1‑71 (MQ=255)
                           aaCTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATg                                     <  1:828309/71‑1 (MQ=255)
                           aaCTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATg                                     <  1:4100161/71‑1 (MQ=255)
                           aaCTTTACGGAATACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATg                                     <  1:4008538/71‑1 (MQ=255)
                                       tACATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATc                                         <  1:45230/55‑1 (MQ=255)
                                         cATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTAttttt                        >  1:2802230/1‑70 (MQ=255)
                                         cATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTAttttt                        >  1:4485465/1‑70 (MQ=255)
                                         cATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTAttttt                        >  1:1714651/1‑70 (MQ=255)
                                         cATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTAttttt                        >  1:569754/1‑70 (MQ=255)
                                         cATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTAttttt                        >  1:3173818/1‑70 (MQ=255)
                                         cATGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTAtttt                         >  1:324546/1‑69 (MQ=255)
                                           tGCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAg                     >  1:4532164/1‑71 (MQ=255)
                                            gCTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAg                     >  1:658262/1‑70 (MQ=255)
                                             cTGTTACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAgg                    >  1:678939/1‑70 (MQ=255)
                                                ttATGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAAtt                <  1:2989748/71‑1 (MQ=255)
                                                ttACGCACCGCAAAATTCCTACCCTTGCCGGTGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAAtt                <  1:3029511/71‑1 (MQ=255)
                                                 tACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAAtt                >  1:304045/1‑70 (MQ=255)
                                                 tACGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAAtt                >  1:4591654/1‑70 (MQ=255)
                                                  aCGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGaa                  >  1:2155409/1‑67 (MQ=255)
                                                  aCGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAAtt                <  1:1503306/69‑1 (MQ=255)
                                                  aCGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATtga              >  1:3908944/1‑71 (MQ=255)
                                                  aCGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATtg               <  1:4258456/70‑1 (MQ=255)
                                                   cGCACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATtgat             <  1:695956/71‑1 (MQ=255)
                                                     cACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATtgatga           <  1:1223116/71‑1 (MQ=255)
                                                     cACCGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATtgatga           <  1:2358902/71‑1 (MQ=255)
                                                       ccGCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCg                                        <  1:1403801/40‑1 (MQ=255)
                                                         gCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttg       >  1:4164591/1‑71 (MQ=255)
                                                         gCAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttg       >  1:2248564/1‑71 (MQ=255)
                                                          cAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttgt      <  1:2112958/71‑1 (MQ=255)
                                                          cAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttgt      <  1:3939501/71‑1 (MQ=255)
                                                          cAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttgt      <  1:3062470/71‑1 (MQ=255)
                                                          cAAAATTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttgt      <  1:2453061/71‑1 (MQ=255)
                                                           aaaaTTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttgt      <  1:3443336/70‑1 (MQ=255)
                                                              aTTACTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACt  >  1:3420365/1‑71 (MQ=255)
                                                                 |                                                                   
TTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACT  >  W3110S.gb/887882‑888014

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: