Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,481,007 C→A 9.5% V246V (GTG→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,0070CA 9.5% 126.1 / 3.5 42V246V (GTG→GTT) rhlBATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base A (2/2);  total (21/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 = 6.66e-01

CCAGGTTCATCACATCGTCGGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAGCTGCTTGTCGAAGTCCATGCCACCGACGAAGCT  >  NC_002947/1480871‑1481150
                                                                                                                                        |                                                                                                                                               
ccAGGTTCATCACATCGTCGGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGa                                                                                                                                    >  2:511601/1‑150 (MQ=255)
     ttCATCACATCGTCGGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATg                                                                                                                               >  1:279438/1‑150 (MQ=255)
             aTCGTCGGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGAACCATCACCTCGACCATGTCCAGGTg                                                                                                                       >  2:397545/1‑150 (MQ=255)
               cgtcgGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCTACCATACGGTCTGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCa                                                                                                                     <  1:511601/150‑1 (MQ=255)
                gtcgGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCAc                                                                                                                    <  2:279438/150‑1 (MQ=255)
                gtcgGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCAc                                                                                                                    >  2:920461/1‑150 (MQ=255)
                gtcgGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCAc                                                                                                                    >  2:319531/1‑150 (MQ=255)
                   ggTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTc                                                                                                                 <  1:1196418/150‑1 (MQ=255)
                             gCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGc                                                                                                                                                 >  2:838205/1‑108 (MQ=255)
                             gCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGc                                                                                                                                                 <  1:838205/108‑1 (MQ=255)
                                       gcagGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGAACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCca                                                                                             <  1:397545/150‑1 (MQ=255)
                                                     tcgctTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCAc                                                                                                                    >  1:70922/1‑113 (MQ=255)
                                                     tcgctTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCAc                                                                                                                    <  2:70922/113‑1 (MQ=255)
                                                     tcgctTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTg                                                                                                        >  2:1139340/1‑125 (MQ=255)
                                                     tcgctTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTg                                                                                                        <  1:1139340/125‑1 (MQ=255)
                                                     gcgctTTGAGGCGGGGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTGAGATGCGTGAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGc                                                                               <  2:41148/149‑1 (MQ=255)
                                                               gCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGccacca                                                                     <  1:319531/150‑1 (MQ=255)
                                                               gCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGccacca                                                                     <  1:920461/150‑1 (MQ=255)
                                                                  gtgtCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGaa                                                                  <  2:1014871/150‑1 (MQ=255)
                                                                           gAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGACAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGt                                                         >  2:709204/1‑150 (MQ=255)
                                                                            aaTGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTg                                                        >  1:799878/1‑150 (MQ=255)
                                                                                     gCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCAt                                                                                                                                <  2:488508/69‑1 (MQ=255)
                                                                                     gCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCAt                                                                                                                                >  1:488508/1‑69 (MQ=255)
                                                                                     gCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATc                                                              >  1:241523/1‑135 (MQ=255)
                                                                                     gCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATc                                                              <  2:241523/135‑1 (MQ=255)
                                                                                                   tGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAg                                            <  2:1000526/139‑1 (MQ=255)
                                                                                                   tGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAg                                            >  1:1000526/1‑139 (MQ=255)
                                                                                                      aGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCCACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAgctgct                              <  1:434987/150‑1 (MQ=255)
                                                                                                         ccATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAgctgc                               >  2:993526/1‑146 (MQ=255)
                                                                                                         ccATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAgctgc                               <  1:993526/146‑1 (MQ=255)
                                                                                                         ccATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACATCGACCATGTCCAGGTGCACTTCACCAAGCTGGTTGAAGTCCAGCAGGAGGCACGGCGTAGACACCAGAATATTGAGGTGGAGCGCTTCTAGGGATTTGAGGTGGTTGt                           >  2:436725/1‑150 (MQ=255)
                                                                                                             gTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAGCTGCTTGTCGaa                       >  1:1123606/1‑150 (MQ=255)
                                                                                                                gAGCATACGGTCGGCTTCGGCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCACAGTGGCGCGCTTCCAGGGCTTTGAGCTGCTTGTCGAAGTc                    >  2:81226/1‑150 (MQ=255)
                                                                                                                 aGCATACGGTCGGCTTCGTCAAGAACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAgc                                  >  1:965884/1‑135 (MQ=255)
                                                                                                                 aGCATACGGTCGGCTTCGTCAAGAACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAgc                                  <  2:965884/135‑1 (MQ=255)
                                                                                                                  gCATACGGTCGGCTTCGTCAAGCACCATCCCCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAANATCGCAGTGGCGCGCTTCCAGGGCTTTGAGCTGCTTGTCGAAGTCCa                  <  1:436725/150‑1 (MQ=255)
                                                                                                                             gCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGgc                                                      <  2:245832/103‑1 (MQ=255)
                                                                                                                             gCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGgc                                                      >  1:245832/1‑103 (MQ=255)
                                                                                                                              cTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAgctgc                               >  1:815707/1‑125 (MQ=255)
                                                                                                                              cTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAgctgc                               <  2:815707/125‑1 (MQ=255)
                                                                                                                               ttCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAGCTGCTTGTCGAAGTCCATGCCACCGACGaa     >  1:910772/1‑150 (MQ=255)
                                                                                                                                  gTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAGCTGCTTGTCGAAGTCCATGCCACCGACGAAGCt  <  2:1123606/150‑1 (MQ=255)
                                                                                                                                   tCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCa                                                           >  2:1111313/1‑92 (MQ=255)
                                                                                                                                   tCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCa                                                           <  1:1111313/92‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                               
CCAGGTTCATCACATCGTCGGTGAAGGTGGCGGAGAACAGCAGGGTCTGGCGCTCGCTTTTAGGCGGTGTCTGACGAATGATTTGGCGGACCTGGGGGATGAAGCCCATGTCGAGCATACGGTCGGCTTCGTCAAGCACCATCACCTCGACCATGTCCAGGTGCACTTCACCACGCTGGTTGAAGTCCAGCAGGCGGCCCGGCGTAGCCACCAGAATATCGCAGTGGCGCGCTTCCAGGGCTTTGAGCTGCTTGTCGAAGTCCATGCCACCGACGAAGCT  >  NC_002947/1480871‑1481150

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