Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	756,577	T→C	10.5%	S391S (TCT→TCC)	PP_0646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	756,577	0	T	C	10.5%	98.3 / 3.4	39	S391S (TCT→TCC)	PP_0646	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/21);  new base C (2/2);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GACAAAGTGGCGAAGCTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTATGTGGCATTGCGCTGGTCGG  >  NC_002947/756437‑756720                                                                                                                                              |                                                                                                                                                 gACAAAGTGGCGAAGCTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAg                                                                                                                                           <  7:221515/148‑1 (MQ=255)    aaaGTGGCGAAGCTTAGTAAGTACCTTCGGTATGGTGGGAATTTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGATGCAAAAAAA‑TTCGTATGTCTGAGGCAggatc                                                                                                                                       >  5:83337/1‑146 (MQ=255)      aGTGGCGAAGCTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCtt                                                                                                                                     <  3:45617/149‑1 (MQ=255)          gCGAAGCTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGc                                                                                                                                 >  6:356027/1‑149 (MQ=255)          gCGAAGCTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGc                                                                                                                                 >  4:123700/1‑149 (MQ=255)              aGCTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑ggg                                                                                                                            <  3:221336/149‑1 (MQ=255)                cTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑ggggg                                                                                                                          <  8:278055/149‑1 (MQ=255)                    gtatttaCCTTCGGTATGGTGGGCATGTCGGAATTGTCTTAGGTGTTACATATTCCTTTCTTAAAGTCCAAGAGGCCTGCCGTGCTGGGGAAACCGAAGGGTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTg                                                                                                                      <  6:83337/144‑1 (MQ=255)                          aCCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGGGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAAACGAAGGCTTCAAAAAAA‑TAAGTTTGTCTGAGGCAGGTGCTTTCGCT‑gggggacatcgc                                                                                                                   >  3:303027/1‑139 (MQ=255)                           ccTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCgggag                                                                                                               <  6:50347/149‑1 (MQ=255)                             ttCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATATTACTATCTGAAAGTCCAAGAGGAATGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGATGGCgggagg                                                                                                              >  3:362301/1‑148 (MQ=255)                                  tATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCgt                                                                                                        <  1:247089/149‑1 (MQ=255)                                     ggtggGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAATTTCGTTTGTCTGAGGCAGGTGCTTTCGTTGGGGGGGCTGGCGGGAGGGATCgt                                                                                                        <  1:94501/148‑1 (MQ=255)                                     ggtggGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCgtggg                                                                                                     >  3:33585/1‑149 (MQ=255)                                                   aaTTGGCTTAGGTGGTACATCTTCCTATCTGAACGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCCGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTACTg                                                                                       <  4:361539/149‑1 (MQ=255)                                                             ggtggtACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGTGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTCGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTgc                                                                             <  3:123700/149‑1 (MQ=255)                                                                 gtACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGcgctc                                                                         >  3:68780/1‑149 (MQ=255)                                                                 gtACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGcgctc                                                                         >  4:85008/1‑149 (MQ=255)                                                                 gtACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCG‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGcgctc                                                                         <  5:329753/149‑1 (MQ=255)                                                                     aTCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTg                                                                     <  4:339311/149‑1 (MQ=255)                                                                      tCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTg                                                                     >  5:301228/1‑148 (MQ=255)                                                                             aTCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTggg                                                             <  5:356027/149‑1 (MQ=255)                                                                                 gAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGgtgt                                                         >  2:325137/1‑149 (MQ=255)                                                                                     gTCCAACACGCCTGCCGTGCAGGTGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCAGAGGCAGGGGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTtctc                                                     <  6:317707/149‑1 (MQ=255)                                                                                        cAAGAGGCCTGCCGTGCAGGGGAATCCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATGGGGGGTGGGAAGATTGGTGGGTTCACTGGGGTCGCTGGATGTGGGGTGTTCTgtgc                                                  >  6:296282/1‑149 (MQ=255)                                                                                                   ccGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACcgc                                         <  4:68780/147‑1 (MQ=255)                                                                                                         aGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCg                                 >  8:78007/1‑149 (MQ=255)                                                                                                         aGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCg                                 <  4:56619/149‑1 (MQ=255)                                                                                                           gggAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCggg                               <  3:139443/149‑1 (MQ=255)                                                                                                             gAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACAGCGGGTTTCGGGgc                             >  6:89619/1‑149 (MQ=255)                                                                                                                 ccGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCg                         <  1:325137/149‑1 (MQ=255)                                                                                                                    aaGGCTGCAAAAAAA‑TTCGTTTGTCCGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTa                      >  5:146332/1‑149 (MQ=255)                                                                                                                    aaGGCTGCAAAAAAA‑TTCGTTTGTCCGAGGCAGGTGCTTTCGCT‑GGGGGGCGGGAGGGAGGTATCGTGGGTGGCAAGATTTCTGGGATCACTGCGCTCGCTTGATGTGGTGTGTTCTCTGCTGGTACCTCTGGTTTCGGTGCGCCGGTa                      >  5:140606/1‑149 (MQ=255)                                                                                                                       gCTGCAATAATA‑TTCGTTGGTCTGAGGCAGGTGCTTTCGCT‑GGGTGTCTGGCGGGATGGATCTTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTAtgt                   <  4:362301/149‑1 (MQ=255)                                                                                                                        cTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTAtgtg                  <  6:301228/149‑1 (MQ=255)                                                                                                                           cAAAAACA‑TTCGTTTGTCCGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTATGTGGCa               <  6:146332/149‑1 (MQ=255)                                                                                                                           cAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTATGTGGCa               <  3:85008/149‑1 (MQ=255)                                                                                                                           cAAAAAAA‑TTCGTTTGGCCGAGTCAGGTTCTTTCGCG‑GGTGGGCGTGCGGGAGGTAGCGGGGGTGGCAAGATTGCGGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGGACCGCGGGTTTCGGGGCGCCGGTATGTGGCa               <  6:140606/149‑1 (MQ=255)                                                                                                                                          tGTCTGAGGCAGGTGCTTTGGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTATGTGGCATTGCGCTGGTCgg  <  5:89619/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                 GACAAAGTGGCGAAGCTTAGTAAGTACCTTCGGTATGGTGGGCATGTCGGAATTGGCTTAGGTGGTACATCTTCCTATCTGAAAGTCCAAGAGGCCTGCCGTGCAGGGGAAACCGAAGGCTGCAAAAAAA‑TTCGTTTGTCTGAGGCAGGTGCTTTCGCT‑GGGGGGCTGGCGGGAGGGATCGTGGGTGGCAAGATTGCTGGGATCACTGCGCTCGCTGTATGTGGGGTGTTCTCTGCTGGTACCGCGGGTTTCGGGGCGCCGGTATGTGGCATTGCGCTGGTCGG  >  NC_002947/756437‑756720

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