Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,295,092	A→G	8.9%	T192A (ACG→GCG)	PP_4668 →	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,295,092	0	A	G	8.9%	121.0 / 3.3	45	T192A (ACG→GCG)	PP_4668	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (2/2);  total (23/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 = 6.48e-01

TCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTG  >  NC_002947/5294951‑5295235                                                                                                                                               |                                                                                                                                                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)                                                                                                                                               |                                                                                                                                                TCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCC‑GATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTG  >  NC_002947/5294951‑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: ‑