Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,697,938	T→G	41.9%	R228R (CGT→CGG)	pyrC' →	dihydroorotase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,697,938	0	T	G	41.9%	9.4 / 5.0	26	R228R (CGT→CGG)	pyrC'	dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (4/11);  new base G (4/7);  total (8/18)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTC  >  NC_002947/5697808‑5698070                                                                                                                                   |                                                                                                                                     aCCTGGCCCAGGGCGGCCTGGCCCTCGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGgcgcc                                                                                                                                 >  2:506597/1‑135 (MQ=255)  ccTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCAGCCGGGCATCCCGGAGCGCGCCGAAACCGAGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCAc                                                                                                                                >  2:483632/1‑136 (MQ=255)     ggCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTc                                                                                                                             <  2:532688/136‑1 (MQ=255)                    ggCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAgc                                                                                                              >  2:29610/1‑136 (MQ=255)                      cccACGAAGGTGCCATGGCCAGCTTCCTCGGCCGGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCGGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAgcgc                                                                                                            <  1:128667/136‑1 (MQ=255)                         aCGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCAGCCCGATCACCAGCGcccc                                                                                                         >  2:507695/1‑135 (MQ=255)                                         cAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGGGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATc                                                                                         <  2:474862/136‑1 (MQ=255)                                         cAGCTGGCTGGGCCGGCGGGGCATCCGGGTGCGCGCCGAAACGGGGGCCCGGGGGCGCAACCTGGGGGTGGGGGAACAGAGCGGCGTGCGGGGGCAGTGCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATc                                                                                         <  1:161084/136‑1 (MQ=255)                                          cgCTTCGTCGGCCTGCGGGGCATCCCGTTGTGCGCCGAAACGGTGGCCCTGGCGCGAAACCTGCTGCTGGGGCAACAGCGGGGGGTGGGGGCGCACTCCCGCCAGATCTCCAGCTCCCGTGGCGCGCGGCTGATCg                                                                                        <  2:272844/135‑1 (MQ=255)                                               ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGTCGCCCTTCAGCCCGAACACCACCGCCCCGGGCGCCCCGCGGATCgccccc                                                                                   >  1:444458/1‑134 (MQ=255)                                               ccTCGGCCTGCCGGGCATCCCGCAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccag                                                                                   <  2:368467/136‑1 (MQ=255)                                                cTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                  <  1:399831/136‑1 (MQ=255)                                                           gggCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                       <  1:167600/136‑1 (MQ=255)                                                           gggCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                       <  1:170550/136‑1 (MQ=255)                                                                             cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACATAAGCCACATAACCAGCCCCCCGGGCGCGCGGCAGGTAGCCCCCGCCCCCGCCATGGGCCTGGCCGGAAccgcc                                                     >  1:294836/1‑136 (MQ=255)                                                                              gAAACCGTGGCCCGGGCGCGCAACCTGCTGCTGGTGGAACAGCGCGGCGTGCGGGCGCACTGCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                    <  2:121770/136‑1 (MQ=255)                                                                                   gggggCCCGGGCGCGCACCCGGCTGCGGGGGCAACAGCGCGGCGGGGGGGCGCACTCCCCCCAGATCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                               <  2:449320/133‑1 (MQ=255)                                                                                        cccTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCGGCCCGATCACCAGCGCCCGTGGCGCGCGGCGGCTCGCCCCGGCCCAGGCCGTGGGCCTGCCGGGAACCGCCGGTGTAgcgcg                                          >  2:338345/1‑135 (MQ=255)                                                                                                cgcAACCTGCGGCTGGTGCAACAGAGCGGCGTGCGTGCGCACTTCACCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGc                                  <  2:501998/136‑1 (MQ=255)                                                                                                      ctgctgctgGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATcc                            >  1:304957/1‑136 (MQ=255)                                                                                                           gctgGTGTAACAGTGGGGGGTGCGTGGGCACTTCTGCCAGATCACCCGCGCCCGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACg                       <  2:490128/136‑1 (MQ=255)                                                                                                               gtggAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACg                   <  2:134184/136‑1 (MQ=255)                                                                                                                         gcggcgTGGGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg         <  1:301194/136‑1 (MQ=255)                                                                                                                         gcggcgTGCGTGCGCACTTCAGCCAGATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg         <  1:525710/136‑1 (MQ=255)                                                                                                                         gcggcgTGCGGGCGCACTCCAGCCAGAACCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg         <  1:188083/136‑1 (MQ=255)                                                                                                                          gggcgTGCGTGCGCACTTCAGCCAGAGCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcgc        <  2:35290/135‑1 (MQ=255)                                                                                                                              gggcgtgcgAATTTCAGCCAGATCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt    <  1:377877/134‑1 (MQ=255)                                                                                                                                gggtgcgCCCATAAGCCAGATAACCCGCGCCCGTGGCGCGCGGCTGATAGCCCCGGCCCCGGGCGTGGGGCGGCCGGGAACCGCCCATGTAGCGCTGTTCCAGCAGATACTCACGGACGAGTCGCGGCGCGAGttg  >  2:453759/3‑135 (MQ=255)                                                                                                                                   |                                                                                                                                     ACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTC  >  NC_002947/5697808‑5698070

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