Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,774,941	T→C	12.9%	S57S (AGT→AGC)	PP_3333 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,774,941	0	T	C	12.9%	76.7 / 5.0	31	S57S (AGT→AGC)	PP_3333	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/13);  new base C (2/2);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

CTACGCCTGCGTGCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGCCCCACGAGCCTGCAAGCCCT  >  NC_002947/3774815‑3775071                                                                                                                               |                                                                                                                                   ctACGCCTGCGTGCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGAAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCaaa                                                                                                                           >  1:432123/1‑136 (MQ=255)     gCCTGCGTGCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGc                                                                                                                       >  2:548756/1‑136 (MQ=255)             gCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGcc                                                                                                               >  2:194448/1‑136 (MQ=255)             gCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGcc                                                                                                               >  1:363713/1‑136 (MQ=255)             gCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGCTTCGGCAAAGGGCTCGGCGcc                                                                                                               >  1:27681/1‑136 (MQ=255)             gCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGCTTCGGCAAAGGGCTCGGCGcc                                                                                                               >  1:30451/1‑136 (MQ=255)                 ttCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCt                                                                                                           <  2:286584/136‑1 (MQ=255)                                     ccTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATc                                                                                       >  2:14256/1‑136 (MQ=255)                                     ccTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATc                                                                                       >  1:510020/1‑136 (MQ=255)                                            cATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGcc                                                                                >  2:307766/1‑136 (MQ=255)                                                   aGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGTGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCtgtg                                                                         <  1:207469/136‑1 (MQ=255)                                                   aGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCtgtg                                                                         <  2:319540/136‑1 (MQ=255)                                                    ggCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGc                                                                        >  2:298336/1‑136 (MQ=255)                                                        ggcggcTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGcc                                                                    >  1:385761/1‑136 (MQ=255)                                                        ggcggcTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGcc                                                                    >  1:217595/1‑136 (MQ=255)                                                             cTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCg                                                               <  2:97694/136‑1 (MQ=255)                                                             cTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCg                                                               <  1:194448/136‑1 (MQ=255)                                                             cTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCg                                                               <  1:123831/136‑1 (MQ=255)                                                                  tGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCt                                                          <  1:532545/136‑1 (MQ=255)                                                                    cTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTgg                                                        <  1:260785/136‑1 (MQ=255)                                                                     tCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGc                                                       <  1:405745/136‑1 (MQ=255)                                                                          ggTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCa                                                  >  1:207847/1‑136 (MQ=255)                                                                            tGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCAcc                                                <  2:363713/136‑1 (MQ=255)                                                                                 gTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCa                                           >  2:356794/1‑136 (MQ=255)                                                                                         ttttGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTgg                                   >  2:43603/1‑136 (MQ=255)                                                                                                         cGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGcccc                   <  2:385761/136‑1 (MQ=255)                                                                                                               cGGCCAGCGCGGTCAGCTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGCCCCACGAGc             <  2:27681/136‑1 (MQ=255)                                                                                                               cGGCCAGCGCGGTCAGCTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGCCCCACGAGc             <  2:30451/136‑1 (MQ=255)                                                                                                                  ccAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGCCCCACGAGCCTg          <  1:69431/136‑1 (MQ=255)                                                                                                                  ccAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGCCCCACGAGCCTg          >  2:89969/1‑136 (MQ=255)                                                                                                                          gTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGCCCCACGAGCCTGCAAGCCCt  <  1:548756/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                   CTACGCCTGCGTGCTGTTCGCAGTGCTGCATTGCGGCCTGGGGCATGGTCAGGCCGGCGGCTTGATGCTCAACGGTGCAGGTGGCACTTTTTGTGGCCACTTAGCGACTGCGGCCAGCGCGGTCAGTTTCGGCAAAGGGCTCGGCGCCACCTTCGCTGGCGTATCCACCATCGATTGCCCGCTGTGCAGCCATGCGGGCCTGGCCATCACCCTGCAGTTGCTGGGCGTGGTATTTGCCCCACGAGCCTGCAAGCCCT  >  NC_002947/3774815‑3775071

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