Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,610,048 C→A S173I (AGC→ATC)  PP_4929 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,610,0480CA96.1% 148.8 / ‑1.8 51S173I (AGC→ATC) PP_4929LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (27/22);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCT  >  NC_002947/5609902‑5610193
                                                                                                                                                  |                                                                                                                                                 
ccTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGc                                                                                                                                                <  2:1477647/150‑1 (MQ=255)
 cTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGcc                                                                                                                                               >  2:92512/1‑150 (MQ=255)
  tGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGccc                                                                                                                                              >  1:82171/1‑150 (MQ=255)
     gggCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAAc                                                                                                                                           >  2:540224/1‑150 (MQ=255)
           aTCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACgg                                                                                                                                     >  1:1191720/1‑150 (MQ=255)
                 tcgccGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTcc                                                                                                                               >  1:667371/1‑150 (MQ=255)
                          gggTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGcc                                                                                                                      >  2:305348/1‑150 (MQ=255)
                                   cTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCtt                                                                                                             <  2:667371/150‑1 (MQ=255)
                                          gtgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCACGGGTGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTgga                                                                                                      <  1:92512/150‑1 (MQ=255)
                                           tgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                     <  2:160819/150‑1 (MQ=255)
                                           tgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTGCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGGGACCATCTCCACGATGGCCAACGCACTGCGCTCCAAGGAAGCCGGGATGCTTTCGTgtag                                                                                                     >  2:305332/1‑150 (MQ=255)
                                             gCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGcc                                                                                                                                               >  1:211738/1‑106 (MQ=255)
                                             gCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGcc                                                                                                                                               <  2:211738/106‑1 (MQ=255)
                                             gCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggagga                                                                                                   >  2:1476636/1‑150 (MQ=255)
                                                 gATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGgaagaa                                                                                               >  1:1333173/1‑150 (MQ=255)
                                                 gATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGgaagaa                                                                                               <  1:1020744/150‑1 (MQ=255)
                                                   tCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTgga                                                                                                      >  2:1105984/1‑141 (MQ=255)
                                                   tCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTgga                                                                                                      <  1:1105984/141‑1 (MQ=255)
                                                    cAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAg                                                                                            <  2:490419/150‑1 (MQ=255)
                                                    cAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAg                                                                                            <  2:1142800/150‑1 (MQ=255)
                                                         aCTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCa                                                                                       <  2:939448/150‑1 (MQ=255)
                                                          cTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGCAGGGAGAACAGGCCCaa                                                                                      >  2:7740/1‑149 (MQ=255)
                                                           tGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGc                                                                                          <  2:1044762/145‑1 (MQ=255)
                                                           tGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGc                                                                                          >  1:1044762/1‑145 (MQ=255)
                                                                gTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGgc                                                                                <  1:305348/150‑1 (MQ=255)
                                                                    cGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGgcgcgc                                                                            <  2:82171/150‑1 (MQ=255)
                                                                     gggCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGGATGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGcgcc                                                                           >  1:28876/1‑150 (MQ=255)
                                                                        ccAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCt                                                                        <  2:1296327/150‑1 (MQ=255)
                                                                        ccAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCt                                                                        >  1:1236804/1‑150 (MQ=255)
                                                                        ccAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCt                                                                        >  2:1387188/1‑150 (MQ=255)
                                                                           aTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTcca                                                                     >  1:977120/1‑150 (MQ=255)
                                                                                 ggCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCt                                                               >  1:1337095/1‑150 (MQ=255)
                                                                                                      ccgccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTg                                                             <  2:972988/131‑1 (MQ=255)
                                                                                                      ccgccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTg                                                             <  2:246631/131‑1 (MQ=255)
                                                                                                      ccgccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTg                                                             >  1:246631/1‑131 (MQ=255)
                                                                                                      ccgccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTg                                                             >  1:972988/1‑131 (MQ=255)
                                                                                                        gccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGt                                                                                                         <  1:1104537/85‑1 (MQ=255)
                                                                                                        gccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGt                                                                                                         >  2:1104537/1‑85 (MQ=255)
                                                                                                          cAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTcc                                      <  1:540224/150‑1 (MQ=255)
                                                                                                                   gCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAg                             >  1:563074/1‑150 (MQ=255)
                                                                                                                      gACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTg                          <  2:1333173/150‑1 (MQ=255)
                                                                                                                            gcACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTc                    >  1:138325/1‑150 (MQ=255)
                                                                                                                              acacCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCGCCTCCa                  >  1:954052/1‑150 (MQ=255)
                                                                                                                              acacCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCa                  >  1:627013/1‑150 (MQ=255)
                                                                                                                                acCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTg                                                                                                        <  2:457417/62‑1 (MQ=255)
                                                                                                                                acCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTg                                                                                                        >  1:457417/1‑62 (MQ=255)
                                                                                                                                      aCCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACggg          <  2:28876/150‑1 (MQ=255)
                                                                                                                                            tcCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCgg    <  2:1337095/150‑1 (MQ=255)
                                                                                                                                              cACGATGCCCAACGCACGGCGCTCCAAGGAATCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCt  >  2:183521/1‑150 (MQ=255)
                                                                                                                                              cACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCg                                                                                  <  2:1096753/70‑1 (MQ=255)
                                                                                                                                              cACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCg                                                                                  >  1:1096753/1‑70 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCT  >  NC_002947/5609902‑5610193

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