Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,402,351	A→T	15.7%	L20Q (CTG→CAG)	PP_2103 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,402,351	0	A	T	15.7%	73.9 / 9.2	32	L20Q (CTG→CAG)	PP_2103	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/12);  new base T (3/2);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

GGCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGC  >  NC_002947/2402206‑2402494                                                                                                                                                   |                                                                                                                                                 ggCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGcagcag                                                                                                                                              >  2:275921/1‑150 (MQ=255)    aCCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAgcgc                                                                                                                                           <  1:158395/150‑1 (MQ=255)           ccccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCTGCTGCGCC‑TGCGCa                                                                                                                                   <  2:247990/150‑1 (MQ=255)             ccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCat                                                                                                                                  >  2:185451/1‑149 (MQ=255)              cacaGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGt                                                                                                                                 >  2:204871/1‑150 (MQ=255)                 aGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcg                                                                                                                              >  1:100454/1‑150 (MQ=255)                   tGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatggaaat                                                                                                                            <  1:185451/150‑8 (MQ=255)                                     cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTgg                                                                                                         >  1:25661/1‑150 (MQ=255)                                      aGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGc                                                                                                        >  2:203568/1‑150 (MQ=255)                                                   gCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcgc                                                                                                                                    <  1:155936/110‑1 (MQ=255)                                                   gCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcgc                                                                                                                                    >  2:155936/1‑110 (MQ=255)                                                           cccATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATgg                                                                                     <  2:31055/148‑1 (MQ=255)                                                           cccATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCCCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATgg                                                                                     >  1:31055/1‑148 (MQ=255)                                                               tACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAg                                                                               >  2:404354/1‑150 (MQ=255)                                                               tACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAg                                                                               <  1:204871/150‑1 (MQ=255)                                                               tACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAg                                                                               >  2:404807/1‑150 (MQ=255)                                                                aCCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAgcgc                                                                                                                                           <  1:151309/90‑1 (MQ=255)                                                                aCCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAgcgc                                                                                                                                           >  2:151309/1‑90 (MQ=255)                                                                      cACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTcc                                                                        <  1:353714/150‑1 (MQ=255)                                                                             cgcAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCt                                                                 <  2:100454/150‑1 (MQ=255)                                                                               cAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTg                                                               <  1:275921/150‑1 (MQ=255)                                                                                aGCGGGTCAAGCAACTGACCGTACTCATGACGGATTTCCTGTCGGTGCATTACCATGCGC‑AGGCGGAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTCGGCATTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTCGCCTTgg                                                              >  2:188126/1‑150 (MQ=255)                                                                                 gCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                             <  1:404354/150‑1 (MQ=255)                                                                                 gCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                             <  1:203568/150‑1 (MQ=255)                                                                                     ttCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGt                                                         >  1:373074/1‑150 (MQ=255)                                                                                         aGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgt                                                     >  1:217158/1‑150 (MQ=255)                                                                                          gCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgtc                                                    >  2:222276/1‑150 (MQ=255)                                                                                                   ctgACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCGGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGACTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgc                                           >  1:399989/3‑150 (MQ=255)                                                                                                               ggctttCGTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCATTGTG‑TTCTGGAGCGGGTGCACGCCCGCATGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCga                               <  1:188126/147‑1 (MQ=255)                                                                                                                           ggTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGaa                   <  1:404807/150‑1 (MQ=255)                                                                                                                              gCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAacac                >  2:183291/1‑150 (MQ=255)                                                                                                                                    aactGCGCAAGGCGCTGCTGCGCC‑TGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACcgccgc          >  2:279473/4‑150 (MQ=255)                                                                                                                                      ctGCGCAAGGCGCTGCTGCGCC‑TGTGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCt        >  1:158514/2‑150 (MQ=255)                                                                                                                                        gcgc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTc      >  2:136192/1‑150 (MQ=255)                                                                                                                                            aaggcgctgCTGCGCC‑TGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg    <  2:325331/147‑1 (MQ=255)                                                                                                                                            aaggcgctgCTGCGCC‑TGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg    >  1:325331/2‑148 (MQ=255)                                                                                                                                             aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGgc  <  2:373074/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGC  >  NC_002947/2402206‑2402494

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