Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 399,522 T→A 25.7% D209V (GAT→GTT)  ddlA ← D‑alanine‑D‑alanine ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb399,5220TA25.7% 121.7 / 44.0 66D209V (GAT→GTT) ddlAD‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base T (33/16);  new base A (11/6);  total (44/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGC  >  W3110S.gb/399464‑399589
                                                          |                                                                   
ctcACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGCCGGCAATTg                                                         >  1:682670/1‑71 (MQ=255)
ctcACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTg                                                         >  1:376563/1‑71 (MQ=255)
ctcACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTg                                                         >  1:2910517/1‑71 (MQ=255)
    cGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTa                                                     <  1:3129575/71‑1 (MQ=255)
    cGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTa                                                     <  1:1095637/71‑1 (MQ=255)
    cGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTa                                                     <  1:3102114/71‑1 (MQ=255)
    cGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTa                                                     <  1:2138847/71‑1 (MQ=255)
    cGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTa                                                     <  1:2921525/71‑1 (MQ=255)
      aCCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCaa                                                            >  1:1438317/1‑62 (MQ=255)
      aCCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACt                                                   >  1:1365547/1‑71 (MQ=255)
      aCCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACt                                                   >  1:436447/1‑71 (MQ=255)
      aCCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACt                                                   >  1:27948/1‑71 (MQ=255)
             aTCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGttctt                                             <  1:1664043/70‑1 (MQ=255)
             aTCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGttctt                                             <  1:1649021/70‑1 (MQ=255)
             aTCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGttctt                                             <  1:2234671/70‑1 (MQ=255)
             aTCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGttctt                                             <  1:3166801/70‑1 (MQ=255)
             aTCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGttctt                                             <  1:2387059/70‑1 (MQ=255)
             aTCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGttctt                                             <  1:1583832/70‑1 (MQ=255)
                ccTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACt                                         >  1:1603504/1‑71 (MQ=255)
                ccTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACt                                         >  1:2375584/1‑71 (MQ=255)
                ccTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACt                                         >  1:2809415/1‑71 (MQ=255)
                ccTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACt                                         >  1:2198369/1‑71 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGc                              >  1:2058307/1‑70 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             >  1:1537025/1‑71 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             >  1:309908/1‑71 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             >  1:1926837/1‑71 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:599618/71‑1 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:838909/71‑1 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:2064708/71‑1 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:113838/71‑1 (MQ=255)
                            aTCACTTTATGATCGAACTCGAACGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:1760438/71‑1 (MQ=255)
                              cACTTTATGATCGAACTCGAACGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:1893464/69‑1 (MQ=255)
                                   tATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCAGTGGTTACTTTGCTGACAcc                       >  1:1927107/1‑70 (MQ=255)
                                   tATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACAcc                       >  1:1084287/1‑70 (MQ=255)
                                   tATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACAcc                       >  1:697863/1‑70 (MQ=255)
                                   tATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACAcc                       >  1:35520/1‑70 (MQ=255)
                                   tATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACAcc                       >  1:524703/1‑70 (MQ=255)
                                   tATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACAcc                       >  1:1282736/1‑70 (MQ=255)
                                   tATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACAcc                       >  1:140402/1‑70 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCTCTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:1913070/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACt                                         >  1:114795/1‑40 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCa                      >  1:758132/1‑59 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:1597430/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:2837687/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:1438686/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:433678/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:463330/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:1910915/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:213041/1‑71 (MQ=255)
                                               cgaacGCCAGAACGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc          >  1:1914799/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGc                              >  1:2041622/1‑43 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:806838/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:70632/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:1064935/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:1055776/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:3045562/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:2719302/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:2614970/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:2512951/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:2310461/1‑71 (MQ=255)
                                                       aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGc  >  1:2113485/1‑71 (MQ=255)
                                                        gATCGACTGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:1198336/43‑1 (MQ=255)
                                                        gATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:3068109/43‑1 (MQ=255)
                                                        gATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:2802274/43‑1 (MQ=255)
                                                        gATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:1810356/43‑1 (MQ=255)
                                                        gATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCt                             <  1:699589/43‑1 (MQ=255)
                                                          tCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTg                            >  1:3328007/1‑42 (MQ=255)
                                                          tCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGa                           >  1:1594107/1‑43 (MQ=255)
                                                          tCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGa                           >  1:1636043/1‑43 (MQ=255)
                                                          tCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGa                           >  1:3266588/1‑43 (MQ=255)
                                                          tCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGa                           >  1:976471/1‑43 (MQ=255)
                                                          |                                                                   
CTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGC  >  W3110S.gb/399464‑399589

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: