Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,295,089 G→C 8.7% V191L (GTG→CTG)  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,0890GC 8.7% 128.2 / 3.9 46V191L (GTG→CTG) PP_4668LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base C (2/2);  total (23/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 = 6.88e-01

ATCAGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTG  >  NC_002947/5294942‑5295235
                                                                                                                                                    |                                                                                                                                                  
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)
                                                                                                                                                 cAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATcccc                                                           <  2:90702/93‑1 (MQ=255)
                                                                                                                                                 cAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATcccc                                                           >  1:90702/1‑93 (MQ=255)
                                                                                                                                                  aGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTg  >  6:293597/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATCAGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTG  >  NC_002947/5294942‑5295235

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: ‑