Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,082,281 T→A 7.3% *318L (TAA→TTA)  PP_1861 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,082,2810TA 7.3% 151.1 / 3.1 55*318L (TAA→TTA) PP_1861LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (22/29);  new base A (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

TCAGCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCT  >  NC_002947/2082137‑2082421
                                                                                                                                                |                                                                                                                                            
tcaGCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATc                                                                                                                                          <  6:157951/149‑1 (MQ=255)
      aCGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAATTTCTGGCACGCTGATCATCAGTGCGCGATTAATTCGCGCac                                                                                                                                    <  4:312553/149‑2 (MQ=255)
      aCGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGATTAATTCGCGCac                                                                                                                                    <  5:346802/149‑2 (MQ=255)
              gAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTgc                                                                                                                            <  5:376340/149‑1 (MQ=255)
              gAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTgc                                                                                                                            >  5:138216/1‑149 (MQ=255)
                   cAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCaa                                                                                                                        <  5:335705/148‑1 (MQ=255)
                      ggCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGc                                                                                                                     >  2:32395/1‑148 (MQ=255)
                          aaaCAGACCTTGTTGTCCAGCAGCAGCTCGTCACAGGAGTCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGt                                                                                                                <  7:44086/149‑1 (MQ=255)
                           aaCAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTc                                                                                                               >  8:323191/1‑149 (MQ=255)
                           aaCAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTc                                                                                                               >  7:149218/1‑149 (MQ=255)
                                aCCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAgg                                                                                                          >  5:57818/1‑149 (MQ=255)
                                aCCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAgg                                                                                                          <  7:312358/149‑1 (MQ=255)
                                  cTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGTGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGt                                                                                                        <  2:332566/149‑1 (MQ=255)
                                  cTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTATTc                                                                                                                                          <  5:195472/115‑1 (MQ=255)
                                  cTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAGTc                                                                                                                                          >  6:195472/1‑115 (MQ=255)
                                  cTGGGTGTCCAGCAGCAGCTAGTCACAGGAGGCTTGGGTTTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCCCGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGt                                                                                                        <  4:64067/149‑1 (MQ=255)
                                   tGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTg                                                                                                       >  6:174521/1‑149 (MQ=255)
                                     gTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTg                                                                                                       <  8:149218/147‑1 (MQ=255)
                                      ttGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGaa                                                                                                    <  3:330651/149‑1 (MQ=255)
                                            agcagcagcTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGc                                                                                              >  3:150674/1‑149 (MQ=255)
                                                    cTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCgg                                                                                      <  6:138216/149‑1 (MQ=255)
                                                       gTCACAGGAGGCTTGGGCGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGctc                                                                                   <  7:44636/149‑1 (MQ=255)
                                                        tCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCAGGCTGATCATAAGTGCGCGATTTAATCGCGCAATTCACTGCGCAAGGACAGGTCCCAGGGTGGAATCGGGCTGAAGCGGctct                                                                                  <  1:107494/149‑1 (MQ=255)
                                                           caGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGa                                                                               <  8:11251/149‑1 (MQ=255)
                                                                   ttGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTc                                                                        >  6:318953/1‑148 (MQ=255)
                                                                   ttGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTc                                                                        >  8:278406/1‑148 (MQ=255)
                                                                     gggTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGTATCGGGCTGAAGCGGCTCTTGAGGAATTCg                                                                       <  1:32395/147‑1 (MQ=255)
                                                                      ggTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGagc                                                                    <  6:255060/149‑1 (MQ=255)
                                                                       gtgtCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGagca                                                                   <  5:318953/149‑1 (MQ=255)
                                                                            ggCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCaaa                                                              >  4:254909/1‑149 (MQ=255)
                                                                                  cATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTg                                                        >  6:164365/1‑149 (MQ=255)
                                                                                   aTGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGa                                                                               <  1:43000/125‑1 (MQ=255)
                                                                                   aTGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGa                                                                               >  2:43000/1‑125 (MQ=255)
                                                                                   aTGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGc                                                       <  6:57818/149‑1 (MQ=255)
                                                                                            cTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTc                                               <  5:164365/148‑1 (MQ=255)
                                                                                             ttGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGt                                             >  5:118183/1‑149 (MQ=255)
                                                                                             ttGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGt                                             >  3:270882/1‑149 (MQ=255)
                                                                                                      atcaATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGc                                    <  3:26458/149‑1 (MQ=255)
                                                                                                         aaTGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAg                                 <  4:150674/149‑1 (MQ=255)
                                                                                                          aTGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAgg                                >  3:337826/1‑149 (MQ=255)
                                                                                                           tGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGc                               >  5:283647/1‑149 (MQ=255)
                                                                                                              aaaaTTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTa                            <  6:283647/149‑1 (MQ=255)
                                                                                                                  ttCTTGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAgcgc                        <  4:216005/149‑1 (MQ=255)
                                                                                                                  ttCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAgcgc                        >  7:198051/1‑149 (MQ=255)
                                                                                                                        cACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAAt                  <  4:337826/149‑1 (MQ=255)
                                                                                                                               atcatcAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCt            <  7:323191/148‑1 (MQ=255)
                                                                                                                                  atcaGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAGACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGt         <  5:174521/148‑1 (MQ=255)
                                                                                                                                   tcaGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTGCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGtct       >  5:1261/1‑149 (MQ=255)
                                                                                                                                     attGCGCGATTAATTCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGtct       >  7:151789/3‑147 (MQ=255)
                                                                                                                                     attGCGCGATTAATTCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGACAGGAGTAGCGCGTAAATGCCGCTGGtct       >  3:253334/3‑147 (MQ=255)
                                                                                                                                      gTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCgaga                                                   <  4:144786/102‑1 (MQ=255)
                                                                                                                                      gTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCgaga                                                   >  3:144786/1‑102 (MQ=255)
                                                                                                                                        gcgcGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCt  >  8:352584/1‑149 (MQ=255)
                                                                                                                                            gAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGtt                                            <  6:353013/103‑1 (MQ=255)
                                                                                                                                            gAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGtt                                            >  5:353013/1‑103 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TCAGCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCT  >  NC_002947/2082137‑2082421

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