Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,402,360	Δ1 bp	22.9%	coding (50/330 nt)	PP_2103 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,402,360	0	T	.	22.9%	86.2 / 23.9	35	coding (50/330 nt)	PP_2103	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/12);  new base . (4/4);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTA  >  NC_002947/2402216‑2402507                                                                                                                                                  |                                                                                                                                                     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:404807/1‑150 (MQ=255)                                                     tACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAg                                                                                            >  2:404354/1‑150 (MQ=255)                                                     tACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAg                                                                                            <  1:204871/150‑1 (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:203568/150‑1 (MQ=255)                                                                       gCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                                          <  1:404354/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                 >  1:325331/2‑148 (MQ=255)                                                                                                                                  aaggcgctgCTGCGCC‑TGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg                 <  2:325331/147‑1 (MQ=255)                                                                                                                                   aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGgc               <  2:373074/150‑1 (MQ=255)                                                                                                                                            tgcgcC‑TGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCa     <  1:279473/149‑1 (MQ=255)                                                                                                                                                ccTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTa  <  2:217158/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGT‑TACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTA  >  NC_002947/2402216‑2402507

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