Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,671,661 G→A 29.2% 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.2% 79.6 / 40.4 55R204Q (CGG→CAG) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (16/22);  new base A (7/9);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGA  >  W3110S.gb/3671597‑3671721
                                                                |                                                            
ttCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCg                                                              <  1:1121278/65‑1 (MQ=255)
  ccgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGt                                                       >  1:244755/1‑70 (MQ=255)
   cgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTg                                                      <  1:3030077/70‑1 (MQ=255)
   cgccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTg                                                      <  1:2157654/70‑1 (MQ=255)
    gccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCTGGTACGTg                                                      <  1:2311666/69‑1 (MQ=255)
    gccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGaa                                                    <  1:4315532/71‑1 (MQ=255)
     ccgTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAAc                                                   <  1:1163992/71‑1 (MQ=255)
          tGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCg                                              >  1:3076658/1‑71 (MQ=255)
          tGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCg                                              >  1:2267598/1‑71 (MQ=255)
          tGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCg                                              >  1:2462812/1‑71 (MQ=255)
                cTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGa                                                     <  1:2018832/58‑1 (MQ=255)
                 tGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACg                                                        >  1:3310439/1‑54 (MQ=255)
                   cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAg                                        >  1:4175115/1‑68 (MQ=255)
                   cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAg                                        >  1:763092/1‑68 (MQ=255)
                   cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                      <  1:3025088/70‑1 (MQ=255)
                   cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                      <  1:146860/70‑1 (MQ=255)
                   cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGc                                       <  1:1937902/69‑1 (MQ=255)
                   cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGc                                       <  1:224891/69‑1 (MQ=255)
                   cAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGc                                       <  1:1294114/69‑1 (MQ=255)
                    aaaCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGaa                                                    >  1:2204182/1‑55 (MQ=255)
                    aaaCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                      <  1:4610903/69‑1 (MQ=255)
                    aaaCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGa                                    <  1:4313378/71‑1 (MQ=255)
                    aaaCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGa                                    <  1:3482841/71‑1 (MQ=255)
                     aaCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAAc                                                   >  1:4620052/1‑55 (MQ=255)
                     aaCCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGAt                                   >  1:2976961/1‑71 (MQ=255)
                        cAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGaa                                <  1:1994348/71‑1 (MQ=255)
                         aGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAAc                               <  1:4190319/71‑1 (MQ=255)
                          gcgcCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCa                                      <  1:3035332/63‑1 (MQ=255)
                            gcCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAAc                               <  1:3193312/68‑1 (MQ=255)
                              cTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCg                                          <  1:3303692/55‑1 (MQ=255)
                              cTCAACATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                             >  1:945581/1‑68 (MQ=255)
                               tCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAAc                               <  1:2611231/65‑1 (MQ=255)
                               tCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                             >  1:4141877/1‑67 (MQ=255)
                                 aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                        <  1:861380/70‑1 (MQ=255)
                                 aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                        >  1:1818504/1‑70 (MQ=255)
                                 aaCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                        >  1:4443547/1‑70 (MQ=255)
                                 aaCATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                        >  1:1483438/1‑70 (MQ=255)
                                 aaCATGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                        >  1:1211538/1‑70 (MQ=255)
                                  aCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGaa                      <  1:1959441/71‑1 (MQ=255)
                                  aCATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGaa                      >  1:4568385/1‑71 (MQ=255)
                                    aTGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAata                    >  1:2883804/1‑71 (MQ=255)
                                     tGCTGTTCTCCGCCACGCTTTCGTACCAGGTACGTg                                                      <  1:3115038/36‑1 (MQ=255)
                                      gCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAatat                   >  1:1618580/1‑70 (MQ=255)
                                       cTGTTCTCCGCCACTCATTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCg                        <  1:4593028/64‑1 (MQ=255)
                                          ttCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTg               >  1:2826228/1‑70 (MQ=255)
                                           gctcCGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGaa             <  1:2837735/70‑1 (MQ=255)
                                            ctcCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGAt                                   <  1:471907/48‑1 (MQ=255)
                                               cGCCACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTg          >  1:4102145/1‑70 (MQ=255)
                                                 ccACGCTTTCGTACCGGGTACGTGAACTGGCGTTCg                                          <  1:1989309/36‑1 (MQ=255)
                                                 ccACGCTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGa                                 <  1:3844257/45‑1 (MQ=255)
                                                     gCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCgg   <  1:1879369/71‑1 (MQ=255)
                                                      cTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                             >  1:183189/1‑44 (MQ=255)
                                                      cTTTCGTACCAGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGa  <  1:4173578/71‑1 (MQ=255)
                                                        ttCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGa              <  1:240348/57‑1 (MQ=255)
                                                          cGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACaa                             >  1:4627422/1‑40 (MQ=255)
                                                            tACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGaa                      <  1:253841/45‑1 (MQ=255)
                                                                |                                                            
TTCCGCCGTATGCCGCCTGCAAACCAGCGCCTCAACATGCTGTTCTCCGCCACGCTTTCGTACCGGGTACGTGAACTGGCGTTCGAGCAGATGAACAATGCCGAATATATTGAAGTGGAACCGGA  >  W3110S.gb/3671597‑3671721

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: