Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,549,041	G→A	12.4%	V125I (GTC→ATC)	PP_3133 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,549,041	0	G	A	12.4%	99.7 / 4.1	40	V125I (GTC→ATC)	PP_3133	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base A (3/2);  total (21/19)
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.93e-01

GGGTGGCACCTCCAATCACTGGGCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACTGGAGCAGACCCCAGGCCTGGA  >  NC_002947/3548899‑3549186                                                                                                                                                |                                                                                                                                                  gggTGGCACCTCCAATCACTGGGCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCAGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATc                                                                                                                                             >  6:96933/1‑149 (MQ=255)        aCCTCCAATCACTGGGCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CATCGATCTGGCGcc                                                                                                                                      >  2:130279/1‑149 (MQ=255)        aCCTCCAATCACTGGGCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CATCGATCTGGCGcc                                                                                                                                      >  2:125690/1‑149 (MQ=255)              aaTCACTGGGCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTc                                                                                                                                <  5:154182/149‑1 (MQ=255)                  aCTGGGCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGa                                                                                                                             <  1:107854/148‑1 (MQ=255)                      ggCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CATCGATCTGGCGCCTGGCTCGGATTTTc                                                                                                                        <  1:125690/149‑1 (MQ=255)                      ggCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CATCGATCTGGCGCCTGGCTCGGATTTTc                                                                                                                        <  1:130279/149‑1 (MQ=255)                         tGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTcgcg                                                                                                                     >  8:49896/1‑149 (MQ=255)                         tGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTcgcg                                                                                                                     >  6:301340/1‑149 (MQ=255)                          ggCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTcgcgc                                                                                                                    <  8:155093/149‑1 (MQ=255)                                       cATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGt                                                                                                       <  8:134405/149‑1 (MQ=255)                                         ttCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTc                                                                                                                        <  4:46434/130‑1 (MQ=255)                                         ttCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTc                                                                                                                        >  3:46434/1‑130 (MQ=255)                                             ccccTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCTGGTCTCAAt                                                                                                 >  5:61956/1‑149 (MQ=255)                                              cccTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCCGCAGCCAAGGAGATATAGCGAGCTGGCGCCTGGGTCGGATTTTCGCGCGTTGAATGCCGGGGTCAAt                                                                                                 >  1:253834/1‑149 (MQ=255)                                               ccTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATgg                                                                                               >  2:166401/1‑149 (MQ=255)                                                cTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTgg                                                                                                                                   <  3:35627/112‑1 (MQ=255)                                                cTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTgg                                                                                                                                   >  4:35627/1‑112 (MQ=255)                                                   cGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCg                                                                                             <  1:166401/147‑1 (MQ=255)                                                                    ccGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCg                                                                          >  2:288231/1‑149 (MQ=255)                                                                         aGGCGTCTTGCCGTGCTGTCCAATTCCCTTCTGGGAGATGGATGGCTTTCTGGCGGCTGCCATGGAGTT‑CGTCGATCTGGCGCCTGGGTCGGATTTTCGCGGGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTc                                                                     <  5:301340/149‑1 (MQ=255)                                                                         aGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTc                                                                     <  1:11545/149‑1 (MQ=255)                                                                              tCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGGCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTc                                                                <  5:96933/149‑1 (MQ=255)                                                                                       gCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTc                                                                                                                                >  2:305818/1‑76 (MQ=255)                                                                                       gCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTc                                                                                                                                <  1:305818/76‑1 (MQ=255)                                                                                             cAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTGTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGt                                                 >  1:252581/1‑149 (MQ=255)                                                                                                      aCTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTgcg                                           >  4:294258/1‑146 (MQ=255)                                                                                                        tCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAg                                        >  5:240140/1‑147 (MQ=255)                                                                                                          ggAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGaa                                      <  6:123630/147‑1 (MQ=255)                                                                                                                 gAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTa                              >  7:277833/1‑148 (MQ=255)                                                                                                                 gAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTa                              >  8:46453/1‑148 (MQ=255)                                                                                                                  aaGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTa                              >  6:140428/1‑147 (MQ=255)                                                                                                                         tCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCTGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACt                       <  6:61956/147‑1 (MQ=255)                                                                                                                            gCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACt                       <  6:234861/144‑1 (MQ=255)                                                                                                                            gCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACt                       >  5:234861/1‑144 (MQ=255)                                                                                                                             ccGGCATCAATGGAGAT‑CGTCGTTTTGGCGCCTGGCTCGTATTTTCCCGCTATTAATGTCGGTTTTAATGGCGGTGACTTCTATGTGTACCGTTTTTTGCTTATCACGCCGAAACTTTTTGTGTAGAAGTACGCTACAGCACTGGa                    <  7:46453/146‑1 (MQ=255)                                                                                                                             ccGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACTGGa                    <  6:240140/146‑1 (MQ=255)                                                                                                                             ccGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACTGGAg                   >  1:113431/1‑147 (MQ=255)                                                                                                                                gCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACTGGAGCAGAc              >  8:48653/1‑149 (MQ=255)                                                                                                                                            aT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACTGGAGCAGACCCCAGGCCTGGa  <  3:294258/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  GGGTGGCACCTCCAATCACTGGGCTGGCCGGTGCCGACCATTCACCCCTTCGGATTTTGCCGTCGCCCCGCCAGGCGTCTTGCCGGGCTGGCCAATTCCCTACTCGGAGATCGAAGGCTTTCTGCCGGCAGCCATGGAGAT‑CGTCGATCTGGCGCCTGGCTCGGATTTTCGCGCGATGAATGCCGGTCTCAATGGCGGTGACTTCGAGGCGGACCGTTTTCTGCTCAGCACGCCGACACGTTTTGCGCAGAAGTACGCTACAGCACTGGAGCAGACCCCAGGCCTGGA  >  NC_002947/3548899‑3549186

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