Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,481,493 G→A 8.5% V84V (GTC→GTT)  rhlB ← ATP‑dependent RNA helicase RhlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,481,4930GA 8.5% 126.6 / 2.9 47V84V (GTC→GTT) rhlBATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base A (2/2);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GGCGTCCTGGCCACGCAGGGTGTAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCC  >  NC_002947/1481352‑1481640
                                                                                                                                             |                                                                                                                                                   
ggCGTCCTGGCCACGCAGGGTGTAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGaa                                                                                                                                             <  2:348816/150‑1 (MQ=255)
               cAGGGTGTAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCca                                                                                                                              >  2:631283/1‑150 (MQ=255)
                         cccAGGCCCTGCGCTTGGATCGGCGTGCAGTACGGGGAGGCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGcc                                                                                                                    <  2:561602/150‑1 (MQ=255)
                          ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                   >  1:354335/1‑150 (MQ=255)
                            aGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCtgt                                                                                                                 >  2:389624/1‑150 (MQ=255)
                             ggACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCtgtg                                                                                                                <  1:631283/150‑1 (MQ=255)
                              gACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTc                                                                                                                                       <  2:1122434/126‑1 (MQ=255)
                              gACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTc                                                                                                                                       >  1:1122434/1‑126 (MQ=255)
                                ccTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggct                                                                                                             <  2:926425/150‑1 (MQ=255)
                                 cTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCTGTggaag                                                                                                            <  1:50032/150‑4 (MQ=255)
                                     gcCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgggc                                                                                                        <  1:417118/150‑1 (MQ=255)
                                               gCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGCGGCTGGTCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGGTTGGGCTTGcaatcga                                                                                              >  1:385039/1‑144 (MQ=255)
                                                       tACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGt                                                                                      >  1:1033894/1‑150 (MQ=255)
                                                                cccAGGTCGTGGATGGCGTGCATTAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTg                                                                             >  2:1095780/1‑150 (MQ=255)
                                                                cccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTg                                                                             >  1:614030/1‑150 (MQ=255)
                                                                  cAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGc                                                                           >  1:259864/1‑150 (MQ=255)
                                                                          ggATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCTGTgg                                                                                                               >  1:1106491/1‑106 (MQ=255)
                                                                          ggATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCTGTgg                                                                                                               <  2:1106491/106‑1 (MQ=255)
                                                                           gATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATc                                                                  <  2:91127/150‑1 (MQ=255)
                                                                              ggCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGc                                                                                                                     >  2:985235/1‑96 (MQ=255)
                                                                              ggCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGc                                                                                                                     <  1:985235/96‑1 (MQ=255)
                                                                                 gTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCgg                                                                                                                                     >  2:800845/1‑77 (MQ=255)
                                                                                 gTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCgg                                                                                                                                     <  1:800845/77‑1 (MQ=255)
                                                                                  tGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGt                                                           <  2:854004/150‑1 (MQ=255)
                                                                                                    aGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGt                                         >  1:670414/1‑150 (MQ=255)
                                                                                                    aGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCAGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGACTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGt                                         >  1:285857/1‑150 (MQ=255)
                                                                                                    aGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCCGCGGGTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGt                                         >  1:285840/1‑150 (MQ=255)
                                                                                                       ttGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                   <  1:85266/73‑1 (MQ=255)
                                                                                                       ttGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                   >  2:85266/1‑73 (MQ=255)
                                                                                                       ttGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCAACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGc                                      >  1:122665/1‑150 (MQ=255)
                                                                                                         gAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTg                                                          >  1:718980/1‑128 (MQ=255)
                                                                                                         gAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTg                                                          <  2:718980/128‑1 (MQ=255)
                                                                                                          aaGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCAACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCa                                   <  2:122665/150‑1 (MQ=255)
                                                                                                                  gAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGGGTTTGGGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTggg                           >  2:116792/1‑150 (MQ=255)
                                                                                                                    aGCGGGTTTTGTCTTCCGGCGGCTGGACCCCGACGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGTCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTggggg                         >  2:617157/1‑150 (MQ=255)
                                                                                                                        ggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCggg                     >  2:1020989/1‑150 (MQ=255)
                                                                                                                         gTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgggg                    >  2:1026635/1‑150 (MQ=255)
                                                                                                                         gTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgggg                    <  2:614030/150‑1 (MQ=255)
                                                                                                                          ttttGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGt                   <  2:354335/150‑1 (MQ=255)
                                                                                                                           tttGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGtt                  <  1:1026635/150‑1 (MQ=255)
                                                                                                                                  tCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCTACCGGTGGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATgg           <  2:285840/150‑1 (MQ=255)
                                                                                                                                  tCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATgg           <  2:285857/150‑1 (MQ=255)
                                                                                                                                     tGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGt        >  2:46056/1‑150 (MQ=255)
                                                                                                                                        ggCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGATTCCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCg     >  1:991452/1‑150 (MQ=255)
                                                                                                                                        ggCTCAACCACGCAGTCTTCGGGCTTCCACAGGCTgg                                                                                                                      <  1:1069079/37‑1 (MQ=255)
                                                                                                                                        ggCTCAACCACGAAGTCTTCGGGCTTCCACAGGCTgg                                                                                                                      >  2:1069079/1‑37 (MQ=255)
                                                                                                                                           tCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGTGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGGCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTCCGGcc  >  1:136348/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GGCGTCCTGGCCACGCAGGGTGTAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCC  >  NC_002947/1481352‑1481640

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑