Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,295,085 C→T 9.5% R189R (CGC→CGT)  PP_4668 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,295,0850CT 9.5% 114.3 / 2.8 43R189R (CGC→CGT) PP_4668LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (18/20);  new base T (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCA  >  NC_002947/5294942‑5295229
                                                                                                                                                |                                                                                                                                                
atcaGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCCAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGt                                                                                                                                             <  2:163720/149‑1 (MQ=255)
         tCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCt                                                                                                                                    <  1:42626/149‑1 (MQ=255)
           cAAGCGTGAGGCGTACATCGTCATCGCGCTGGAGCGGACCTAGCATGGTCCTTATTTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTTGACAGCCACGCGCC‑GATCCGCAAGGTGACGGATCTgg                                                                                                                                  <  7:307535/149‑1 (MQ=255)
            aaGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCAGATCCGTCACCTGGCGGATcggc                                                                                                                                  <  4:82378/149‑4 (MQ=255)
                                       cTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGa                                                                                                      >  4:28126/1‑149 (MQ=255)
                                        tGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGTCACCTGGCGGATcggcgc                                                                                                                                <  7:291114/122‑6 (MQ=255)
                                        tGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGTCACCTGGCGGATcggcgc                                                                                                                                >  8:291114/1‑117 (MQ=255)
                                        tGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGAt                                                                                                                                      <  5:244885/116‑1 (MQ=255)
                                        tGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGAt                                                                                                                                      >  6:244885/1‑116 (MQ=255)
                                            gCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGAcc                                                                                                 >  6:198519/1‑149 (MQ=255)
                                                        tGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCt                                                                                     <  8:285245/149‑1 (MQ=255)
                                                        tGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCt                                                                                     >  3:32006/1‑149 (MQ=255)
                                                                   tgtgCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTGCGAGctgctg                                                                          >  3:287205/1‑149 (MQ=255)
                                                                   tgtgCTGCAAGCTGTGCGACTACCGCCTGAGGCTGTATGCCAAACAGGGCTACCTGGACAGCCACGCGCC‑GATCCGACAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACATGGCGTTAAGCTCCGAGctgctg                                                                          >  4:330024/1‑149 (MQ=255)
                                                                    gtgCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGt                                                                         <  8:235340/149‑1 (MQ=255)
                                                                     tgctgcAAGCTGTGTGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTtcatca                                                                                                                >  2:312615/1‑109 (MQ=255)
                                                                     tgctgcAAGCTGTGTGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTtcatca                                                                                                                <  1:312615/109‑1 (MQ=255)
                                                                          cAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTgg                                                                   <  1:344874/149‑1 (MQ=255)
                                                                            aGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTg                                                                                                        >  1:186357/1‑110 (MQ=255)
                                                                            aGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTg                                                                                                        <  2:186357/110‑1 (MQ=255)
                                                                            aGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGc                                                                  >  6:61895/1‑148 (MQ=255)
                                                                                gtgCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGc                                                                  >  5:112688/1‑144 (MQ=255)
                                                                                gtgCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGc                                                                  <  6:112688/144‑1 (MQ=255)
                                                                                 tgCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAAcc                                                             <  5:198519/148‑1 (MQ=255)
                                                                                   cgacacccgcccgcGGCTGCACGCCACCCAGGCCTACCTGGACAGCCACGCGCC‑CATCCGCCAGGTGACCCATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGa                                                          <  4:287205/143‑1 (MQ=255)
                                                                                     aCTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATc                                                        <  4:32006/149‑1 (MQ=255)
                                                                                     aCTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATc                                                        <  6:230410/149‑1 (MQ=255)
                                                                                      cTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATcc                                                       >  3:230838/1‑149 (MQ=255)
                                                                                         ccGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTg                                                           <  3:356241/142‑1 (MQ=255)
                                                                                         ccGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTg                                                           >  4:356241/1‑142 (MQ=255)
                                                                                         ccGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATccc                                                      <  4:230838/147‑1 (MQ=255)
                                                                                          cGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCaa                                                   <  3:28126/149‑1 (MQ=255)
                                                                                                 ggCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAgc                                            >  6:264923/1‑149 (MQ=255)
                                                                                                        tGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACAt                                      <  1:25885/148‑1 (MQ=255)
                                                                                                             cccAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTgcg                                 >  4:240936/1‑148 (MQ=255)
                                                                                                              ccAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTgcg                                 >  5:158069/1‑147 (MQ=255)
                                                                                                                      aCCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGc                       >  3:329952/1‑149 (MQ=255)
                                                                                                                                    gcgcC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTa         >  6:21106/1‑149 (MQ=255)
                                                                                                                                    gcgcC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTCATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTa         >  8:308753/1‑149 (MQ=255)
                                                                                                                                     cgcC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGAt                                                         >  4:315368/1‑100 (MQ=255)
                                                                                                                                     cgcC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGAt                                                         <  3:315368/100‑1 (MQ=255)
                                                                                                                                       nnagATCCGTCACCTGGCGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTaca       >  2:317048/4‑149 (MQ=255)
                                                                                                                                            tCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCa  <  5:264923/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATCAGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCA  >  NC_002947/5294942‑5295229

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