Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,059,725	C→A	15.4%	L436I (CTC→ATC)	PP_3573 →	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,059,725	0	C	A	15.4%	76.3 / 2.7	26	L436I (CTC→ATC)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (11/11);  new base A (2/2);  total (13/13)
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.99e-01

GTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGC  >  NC_002947/4059599‑4059852                                                                                                                               |                                                                                                                                gtgtGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGccc                                                                                                                        <  2:245204/136‑1 (MQ=255)     ggCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTg                                                                                                                    >  2:91475/1‑136 (MQ=255)     ggCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTg                                                                                                                    >  2:413515/1‑136 (MQ=255)     ggCACGACGCCAAGCACGTCGCCGGCCACAACGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGAACAGCAATGTCAGCACCCTCGGAGCCCACtt                                                                                                                    >  2:175344/1‑135 (MQ=255)      gCACGACGCCAAGCACGTCGCCGGCCACATGGCAACACAACGCACATACACCGCCTACCGCGACCTCGAGCAACGCGCAGCGGATGAGCATCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGa                                                                                                                   <  1:457849/136‑1 (MQ=255)                         gCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCAt                                                                                                >  2:308792/1‑136 (MQ=255)                          ccGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATc                                                                                               <  1:413515/136‑1 (MQ=255)                          ccGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATc                                                                                               <  1:91475/136‑1 (MQ=255)                          ccGGCCACATCGCAACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCTCTGATGCCTACCCATACTCGCATc                                                                                               <  1:175344/136‑1 (MQ=255)                                    cGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGTTCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCa                                                                                     <  1:88362/136‑1 (MQ=255)                                                 aCCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCATCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGa                                                                           <  2:351632/133‑1 (MQ=255)                                                 aCCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCATCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGa                                                                           >  1:351632/1‑133 (MQ=255)                                                 aCCTACACCGCCCACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCATCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGa                                                                           <  2:366412/133‑1 (MQ=255)                                                 aCCTACACCGCCCACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCATCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGa                                                                           >  1:366412/1‑133 (MQ=255)                                                         cGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACACATACTCGCATCCGCATGTACAACACAAAGGAAACCTACCCCa                                                                >  1:63226/1‑136 (MQ=255)                                                              accgcgaccgcgAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAg                                                           >  1:167589/1‑136 (MQ=255)                                                              aacgcgacagctaGCAACGCGCAGCGGATGAGCAGCAAACCGCCAACAACAGCACTGTCAGCACCCTCGGAGCCCACTGATCCCTACCCAAACACGCATCCGCATGAACAACACAAAGGAAACATACCCCAACCaa                                                           >  1:271648/3‑135 (MQ=255)                                                                cgcgaccgcgaGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGAc                                                         >  1:229325/1‑136 (MQ=255)                                                                          aGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAAc                                               <  1:448392/136‑1 (MQ=255)                                                                          aGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAAc                                               <  1:150791/136‑1 (MQ=255)                                                                            cAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCATCAATGTCAGCACCCTCTGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACAAGACCCTGGACAACGa                                             <  2:303089/136‑1 (MQ=255)                                                                                         aTGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAACGACTTGTGTCAGGcc                                <  1:135352/136‑1 (MQ=255)                                                                                                 caAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAACGACTTGTGTCAGGCCGTgcgggc                        >  2:500568/1‑136 (MQ=255)                                                                                                                   aaTGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGt      >  1:414696/1‑136 (MQ=255)                                                                                                                     tGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTgc    >  1:12386/1‑136 (MQ=255)                                                                                                                       tCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCATTCCCTGGACATCGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTgcgc  <  2:271648/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                GTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCAGCACCCTCGGAGCCCACTGATGCCTACCCATACTCGCATCCGCATGTTCAACACAAAGGAAACCTACCCCAACCAGACCCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGC  >  NC_002947/4059599‑4059852

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