Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,038,094	C→A	15.4%	R155L (CGA→CTA)	PP_5284 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,038,094	0	C	A	15.4%	79.9 / 3.5	26	R155L (CGA→CTA)	PP_5284	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base A (2/2);  total (14/12)
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.46e-01

AACGCCAGCCCCCGCTTCAACGACAGCGTCAAATGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAACAGACACCCGAACACCTCGCTGGCCTCATGCC  >  NC_002947/6037962‑6038205                                                                                                                                     |                                                                                                                aaCGCCAGCCCCCGCTTCAACGACAGCGTCAAATGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTcgcg                                                                                                              >  1:551619/1‑136 (MQ=255)           cccGCTTCAACGACAGCGTCAAATGCACATCGTCCTGACTCGGTTCTCTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACt                                                                                                    <  2:507561/136‑1 (MQ=255)                ttCAACGACAGCGTCAAATGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTAGCGGATAGACACTCGccc                                                                                               >  2:261593/1‑136 (MQ=255)                ttCAACGACAGCGTCAAATGCACATCATCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTAGCGGATAGACACTCGccc                                                                                               >  2:630899/1‑136 (MQ=255)                      gACAGCGTCAAATGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTa                                                                                         <  2:625560/136‑1 (MQ=255)                             tCAAATGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGcc                                                                                  >  1:196478/1‑136 (MQ=255)                                aaTGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCt                                                                               <  1:426622/136‑1 (MQ=255)                                aaTGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCt                                                                               <  1:494289/136‑1 (MQ=255)                                          gTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGAtc                                                                     >  2:77620/1‑136 (MQ=255)                                                    ggTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCa                                                           >  1:371510/1‑136 (MQ=255)                                                      ttCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGt                                                         <  1:77620/136‑1 (MQ=255)                                                      ttCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTAGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGt                                                         <  1:261593/136‑1 (MQ=255)                                                      ttCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTAGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGt                                                         <  1:630899/136‑1 (MQ=255)                                                        cACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTgg                                                       >  1:335494/1‑136 (MQ=255)                                                        cACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTgg                                                       >  2:42518/1‑136 (MQ=255)                                                                   gAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTgg                                            >  2:635981/1‑136 (MQ=255)                                                                   gAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTgg                                            >  2:8090/1‑136 (MQ=255)                                                                            ttGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGaa                                   >  2:613679/1‑136 (MQ=255)                                                                            ttGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGaa                                   <  2:165879/136‑1 (MQ=255)                                                                            ttGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGaa                                   <  1:613679/136‑1 (MQ=255)                                                                              gtgGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAAc                                 >  1:406909/1‑136 (MQ=255)                                                                               tgGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCACCTGAACAGGATCTCGAATTCCAGTGGCAGGTGCTTGGTATCGAGAAACa                                >  1:479529/1‑136 (MQ=255)                                                                                                    gCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAACAGACACCCGAACACCTCGCTgg           <  2:514517/136‑1 (MQ=255)                                                                                                    gCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAACAGACACCCGAACACCTCGCTgg           <  2:551619/136‑1 (MQ=255)                                                                                                          tgCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAACAGACACCCGAACACCTCGCTGGcc         <  2:317636/132‑1 (MQ=255)                                                                                                          tgCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAACAGACACCCGAACACCTCGCTGGcc         >  1:317636/1‑132 (MQ=255)                                                                                                             aCCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAACAGACACCCGAACACCTCGCTGGCCTCATGcc  >  2:419151/1‑136 (MQ=255)                                                                                                                                     |                                                                                                                AACGCCAGCCCCCGCTTCAACGACAGCGTCAAATGCACATCGTCCTGACTCGGTTCACTGATCCCTGAAGGATGGTTGTGGCACAGGATCAATGCCGCAGCGTTGTGCACCAGCGCCCGCCGCACCACCTCTCGCGGATAGACACTCGCCCGATCTATCGTGCCCCTGAACAGGATCTCGAATGCCAGTGGCCGGTGCTTGGTATCGAGAAACAGACACCCGAACACCTCGCTGGCCTCATGCC  >  NC_002947/6037962‑6038205

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