Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,567,128:1 +C 7.8% coding (843/879 nt) PP_2250 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,567,1281.C 7.8% 197.7 / 11.2 51coding (843/879 nt)PP_2250LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (22/25);  new base C (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

ATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCAAGCC  >  NC_002947/2566993‑2567277
                                                                                                                                        |                                                                                                                                                     
aTGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCgg                                                                                                                                         >  2:442594/1‑150 (MQ=255)
 tctcgcgGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACGGACGTTGCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGc                                                                                                                                        <  2:1190051/147‑1 (MQ=255)
   aCGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCtg                                                                                                                                      <  1:893062/150‑1 (MQ=255)
     gcgGTGCGGGGCCTATCGACGGAGCTATCGGTTTGAGGGGTAAAAGGTTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGTAGCGGCtaga                                                                                                                                    >  1:1190051/1‑147 (MQ=255)
     gcgGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCtgtt                                                                                                                                    <  2:932533/150‑1 (MQ=255)
         tGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGttt                                                                                                                                <  1:1045559/150‑1 (MQ=255)
             gggCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTgcgc                                                                                                                            <  2:662333/150‑1 (MQ=255)
               gCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGACCGCGCCGCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCca                                                                                                                          >  1:73507/1‑150 (MQ=255)
                  tATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGTAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAg                                                                                                                       <  1:131288/150‑1 (MQ=255)
                  tATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTGGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAg                                                                                                                       >  2:808828/1‑150 (MQ=255)
                         cGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGc                                                                                                                <  1:660524/150‑1 (MQ=255)
                               atCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAgg                                                                                                                      >  2:795296/1‑138 (MQ=255)
                               atCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAgg                                                                                                                      <  1:795296/138‑1 (MQ=255)
                               atCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGc                                                                                                          >  2:52777/1‑150 (MQ=255)
                                 cGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAg                                                                                                        <  2:370769/150‑1 (MQ=255)
                                             tAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                        >  1:397522/1‑122 (MQ=255)
                                             tAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                        <  2:397522/122‑1 (MQ=255)
                                              aaaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTg                                                                                           <  1:52777/150‑1 (MQ=255)
                                               aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGc                                                                                          >  1:847318/1‑150 (MQ=255)
                                               aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGc                                                                                          >  1:827577/1‑150 (MQ=255)
                                                 aGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCgt                                                                                        <  2:576388/150‑1 (MQ=255)
                                                 aGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCgt                                                                                        <  2:643881/150‑1 (MQ=255)
                                                   aTTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGc                                                                                          >  2:485766/1‑146 (MQ=255)
                                                   aTTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGc                                                                                          <  1:485766/146‑1 (MQ=255)
                                                          tGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACgg                                                                               >  1:1096657/1‑150 (MQ=255)
                                                            cTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGaa                                                                             >  2:271193/1‑150 (MQ=255)
                                                                gggAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctgat                                                                                                                                             >  1:1050986/1‑80 (MQ=38)
                                                                gggAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctgat                                                                                                                                             <  2:1050986/82‑3 (MQ=38)
                                                                gggAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCAGACTGAAAAGCAGGCGTGCAACCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTcc                                                                         >  1:750399/1‑150 (MQ=255)
                                                                 ggAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCa                                                                        <  1:179306/150‑1 (MQ=255)
                                                                  gAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCTACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAg                                                                       <  2:496779/150‑1 (MQ=255)
                                                                  gAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTcc                                                                         <  1:468440/148‑1 (MQ=255)
                                                                  gAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTcc                                                                         >  2:468440/1‑148 (MQ=255)
                                                                            acacacGTTTCGACGGACGGCATTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGa                                                             <  1:170228/150‑1 (MQ=255)
                                                                              aGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctga                                                                                                                                              <  2:695742/67‑2 (MQ=37)
                                                                              aGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctga                                                                                                                                              >  1:695742/1‑66 (MQ=37)
                                                                                       gacggacgGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTgct                                                  >  1:71348/1‑148 (MQ=255)
                                                                                        acggacggCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGcc                                                 <  2:750399/150‑1 (MQ=255)
                                                                                          ggacggCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCcag                                               >  1:672077/1‑150 (MQ=255)
                                                                                            acggCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCcagca                                             >  2:1032473/1‑150 (MQ=255)
                                                                                                  gTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCgggg                                       <  1:442594/150‑1 (MQ=255)
                                                                                                   tCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGc                                      <  2:472160/150‑1 (MQ=255)
                                                                                                   tCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGc                                      >  1:472160/1‑150 (MQ=255)
                                                                                                   tCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGc                                      >  2:1014433/1‑150 (MQ=255)
                                                                                                       ggATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCca                                                <  1:579871/136‑1 (MQ=255)
                                                                                                       ggATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCca                                                >  2:579871/1‑136 (MQ=255)
                                                                                                                     ctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCAACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAg                    <  2:750752/150‑1 (MQ=255)
                                                                                                                               agGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCAAGCAGCCGGGCGACCCGCGACAACGGCAGATCTCCCCAg          >  2:1152845/1‑150 (MQ=255)
                                                                                                                                    tGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGAACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCCGGCGCACCCGACGCCCCCTGGCCAGCCCCGGGGCGACCCGCGACAACCGCAGATCTCCCCAGCCCaa     <  2:73507/150‑1 (MQ=255)
                                                                                                                                      cAGGCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCAAGc   >  2:301462/1‑150 (MQ=255)
                                                                                                                                      cAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCAAGc   <  1:808828/150‑1 (MQ=255)
                                                                                                                                       aGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCAAGcc  >  2:541439/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                                                                     
ATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCAAGCC  >  NC_002947/2566993‑2567277

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