Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,671,661 G→A 15.9% R204Q (CGG→CAG)  rhlB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,671,6610GA15.9% 109.6 / 12.8 44R204Q (CGG→CAG) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (2/5);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GGCTGTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGAA  >  W3110S.gb/3671592‑3671727
                                                                     |                                                                  
ggctgTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCgg                                                                   >  1:1388745/1‑71 (MQ=255)
 gctgTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCggg                                                                  >  1:2351917/1‑71 (MQ=255)
 gctgTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCggg                                                                  <  1:1640735/71‑1 (MQ=255)
     ttCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTa                                                                <  1:1270166/69‑1 (MQ=255)
     ttCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTa                                                                <  1:1339983/69‑1 (MQ=255)
       ccgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGt                                                             >  1:700711/1‑70 (MQ=255)
       ccgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGt                                                             >  1:668663/1‑70 (MQ=255)
        cgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGTGTACGTGa                                                           <  1:2365284/71‑1 (MQ=255)
        cgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGa                                                           <  1:1974246/71‑1 (MQ=255)
        cgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTg                                                            <  1:422140/70‑1 (MQ=255)
                 ccgccTGCAAACGAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCgg                                                                   >  1:2360981/1‑54 (MQ=255)
                        cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAg                                              >  1:38903/1‑68 (MQ=255)
                        cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAg                                              >  1:2366249/1‑68 (MQ=255)
                        cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAg                                              >  1:523019/1‑68 (MQ=255)
                        cAAACCAGCGCCTCAACATGCAGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAg                                              >  1:579970/1‑68 (MQ=255)
                           aCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGa                                                           <  1:2414010/52‑1 (MQ=255)
                           aCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGa                                                           <  1:1102735/52‑1 (MQ=255)
                              aGCGCCTCAACATGCTGTTCTCCTCCACGCTTTCGTACCGGGTACGTGAACTGTCGTTCGAGCAGATGAAc                                     <  1:298752/71‑1 (MQ=255)
                              aGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGaa                                      >  1:132113/1‑70 (MQ=255)
                              aGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAAc                                     <  1:1259138/71‑1 (MQ=255)
                               gcgcCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAAc                                     <  1:1351502/70‑1 (MQ=255)
                                 gcCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTTAACTGGCGTTCGAGCAGATGAACaa                                   <  1:1654807/70‑1 (MQ=255)
                                    tCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAg                                              <  1:810970/56‑1 (MQ=255)
                                     cAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACa                                    >  1:2355303/1‑65 (MQ=255)
                                      aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATg                                 >  1:503360/1‑67 (MQ=255)
                                      aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                              >  1:2054964/1‑70 (MQ=255)
                                      aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                              >  1:1358583/1‑70 (MQ=255)
                                      aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                              <  1:374005/70‑1 (MQ=255)
                                       aCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGaa                            <  1:721239/71‑1 (MQ=255)
                                        cATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGaa                            <  1:1466644/70‑1 (MQ=255)
                                        cATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGaa                            <  1:1017903/70‑1 (MQ=255)
                                         aTGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAata                          <  1:675031/71‑1 (MQ=255)
                                         aTGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAata                          <  1:17160/71‑1 (MQ=255)
                                         aTGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAata                          <  1:1704255/71‑1 (MQ=255)
                                         aTGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAata                          >  1:753707/1‑71 (MQ=255)
                                         aTGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAata                          <  1:608178/71‑1 (MQ=255)
                                         aTGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTg                                                            <  1:2029209/37‑1 (MQ=255)
                                         aTGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGaa                                                          >  1:1917177/1‑39 (MQ=255)
                                           gCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAatata                        >  1:395238/1‑71 (MQ=255)
                                           gCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGaa                                                          >  1:1643231/1‑37 (MQ=255)
                                               ttCTCCGCCACGCTTTCGTACCGGGTACGTGATCTGGCGTTCGAGCAGATGAACAATGCCGAATATATTg                     <  1:2765686/70‑1 (MQ=255)
                                               ttCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTg                     <  1:1529742/70‑1 (MQ=255)
                                                tctcCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTg                     >  1:1709625/1‑69 (MQ=255)
                                                 ctcCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAg                  >  1:129481/1‑71 (MQ=255)
                                                        aCGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAAc                                     <  1:2671507/45‑1 (MQ=255)
                                                                  aCCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAAc            >  1:26592/1‑60 (MQ=255)
                                                                  aCCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAAc            >  1:171604/1‑60 (MQ=255)
                                                                  aCCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaa  <  1:2340186/70‑1 (MQ=255)
                                                                     |                                                                  
GGCTGTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGAA  >  W3110S.gb/3671592‑3671727

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: