Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,671,661 G→A 29.1% 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,6610GA29.1% 84.4 / 40.0 55R204Q (CGG→CAG) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (25/14);  new base A (9/7);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGAAAA  >  W3110S.gb/3671596‑3671729
                                                                 |                                                                    
gTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTAc                                                                 >  1:1458007/1‑71 (MQ=255)
 ttCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTAc                                                                 <  1:3969709/70‑1 (MQ=255)
   ccgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGt                                                               >  1:1544529/1‑70 (MQ=255)
    cgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTg                                                              <  1:631305/70‑1 (MQ=255)
         tATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTg                                                         <  1:2254440/70‑1 (MQ=255)
         tATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTg                                                         <  1:2929413/70‑1 (MQ=255)
             ccgccTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTa                                                                  >  1:808671/1‑57 (MQ=255)
             ccgccTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTa                                                                  >  1:1613832/1‑57 (MQ=255)
             ccgccTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTa                                                                  >  1:2754994/1‑57 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                              <  1:3618124/70‑1 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                              <  1:1441179/70‑1 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAg                                                >  1:4085141/1‑68 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAg                                                >  1:549788/1‑68 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAg                                                >  1:3632328/1‑68 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAg                                                >  1:3288307/1‑68 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAg                                                >  1:1739300/1‑68 (MQ=255)
                    cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCa                                              <  1:1505287/70‑1 (MQ=255)
                       aCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTg                                                              <  1:2108654/51‑1 (MQ=255)
                       aCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTg                                                              <  1:1907075/51‑1 (MQ=255)
                        ccAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGaa                                                            >  1:1658277/1‑52 (MQ=255)
                        ccAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGaa                                                            >  1:3556278/1‑52 (MQ=255)
                         cAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGaa                                                            >  1:972566/1‑51 (MQ=255)
                             gcCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAAt                                    <  1:110812/71‑1 (MQ=255)
                             gcCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAAt                                    <  1:1778153/71‑1 (MQ=255)
                                tCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                                     >  1:25770/1‑67 (MQ=255)
                                tCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                                     >  1:1694087/1‑67 (MQ=255)
                                tCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                                     >  1:3675772/1‑67 (MQ=255)
                                  aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                              >  1:337070/1‑56 (MQ=255)
                                  aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                              >  1:869155/1‑56 (MQ=255)
                                      tGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGaa                                                            >  1:503820/1‑38 (MQ=255)
                                      tGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAata                            >  1:1523245/1‑70 (MQ=255)
                                        cTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAatatat                         >  1:1332914/1‑71 (MQ=255)
                                           ttCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTg                       <  1:1986258/70‑1 (MQ=255)
                                             ctcCGCGACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGATCa                                      >  1:521316/1‑53 (MQ=255)
                                             ctcCGCCACGCTTTCGTATCGGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                                     >  1:1296119/1‑54 (MQ=255)
                                             ctcCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGa                      <  1:889354/69‑1 (MQ=255)
                                             ctcCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAg                    >  1:3868121/1‑71 (MQ=255)
                                             ctcCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAg                    >  1:1620213/1‑71 (MQ=255)
                                                cGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGaa                     >  1:359536/1‑67 (MQ=255)
                                                 gcgaCGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTgg                 <  1:2140536/67‑1 (MQ=255)
                                                 gCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGa                <  1:166163/71‑1 (MQ=255)
                                                   cACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACa                                      >  1:3656107/1‑47 (MQ=255)
                                                   cACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACa                                      >  1:998929/1‑47 (MQ=255)
                                                    aCGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAAcc             >  1:105675/1‑71 (MQ=255)
                                                      gCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAAt                                    <  1:2427545/46‑1 (MQ=255)
                                                         ttCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGa                      <  1:520690/57‑1 (MQ=255)
                                                           cGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAg                    >  1:3800381/1‑57 (MQ=255)
                                                             tACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaa    <  1:423268/71‑1 (MQ=255)
                                                             tACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaa    <  1:142412/71‑1 (MQ=255)
                                                             tACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaa    <  1:1471425/71‑1 (MQ=255)
                                                             tACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaa    <  1:2977498/71‑1 (MQ=255)
                                                               ccGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaaaa  >  1:1222672/1‑71 (MQ=255)
                                                               ccGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaaaa  >  1:3736708/1‑71 (MQ=255)
                                                               ccGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaaaa  >  1:2941830/1‑71 (MQ=255)
                                                               ccGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGaaaa  >  1:2527463/1‑71 (MQ=255)
                                                                 |                                                                    
GTTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGAACAGAAAA  >  W3110S.gb/3671596‑3671729

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: