Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,420,750	G→A	8.5%	A51V (GCC→GTC)	PP_5554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,420,750	0	G	A	8.5%	133.3 / 2.2	47	A51V (GCC→GTC)	PP_5554	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base A (2/2);  total (23/24)
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.65e-01

GTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCAGTAGCACCACGATGCTGCCGGTTC  >  NC_002947/3420605‑3420873                                                                                                                                                    |                                                                                                                            gtgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                                                                                            >  3:49502/1‑147 (MQ=255) gtgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGa                                                                                                                          <  3:277846/149‑1 (MQ=255)  tgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAg                                                                                                                         <  4:49502/149‑1 (MQ=255)    cacaCCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAgtg                                                                                                                       >  4:283377/1‑149 (MQ=255)         cacGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGGCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAaca                                                                                                                  >  6:25044/1‑149 (MQ=255)         cacGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAaca                                                                                                                  >  6:312522/1‑149 (MQ=255)                        cTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGc                                                                                                    <  6:312286/148‑1 (MQ=255)                         ttGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAgtgt                                                                                                                      <  3:37359/129‑1 (MQ=255)                         ttGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAgtgt                                                                                                                      >  4:37359/1‑129 (MQ=255)                                 gggggTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                                                                                            <  4:221456/115‑1 (MQ=255)                                 gggggTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                                                                                            >  3:221456/1‑115 (MQ=255)                                  ggggTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCa                                                                                         <  3:283377/149‑1 (MQ=255)                                        cAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCtt                                                                                     >  3:12071/1‑147 (MQ=255)                                        cAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTaa                                                                                   <  6:120114/149‑1 (MQ=255)                                           gCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGt                                                                                 >  5:312128/1‑148 (MQ=255)                                               aaGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACg                                                                             >  3:101717/1‑148 (MQ=255)                                                     cTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGG‑AGACATCACTGGTAATAGCTCTTAACGAGCCATGCGACGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGc                                                                       >  1:86489/1‑149 (MQ=255)                                                     cTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGGTCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGACGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATg                                                                        >  1:90438/1‑148 (MQ=255)                                                       gggTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCAt                                                                    >  5:18743/1‑149 (MQ=255)                                                              gccTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCtcgtcg                                                                                                     <  2:59914/109‑1 (MQ=255)                                                              gccTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCtcgtcg                                                                                                     >  1:59914/1‑109 (MQ=255)                                                                        accacAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGc                                                                       >  2:316777/1‑129 (MQ=255)                                                                        accacAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGc                                                                       <  1:316777/129‑1 (MQ=255)                                                                        accacAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGc                                                                <  6:237835/136‑1 (MQ=255)                                                                        accacAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGc                                                                >  5:237835/1‑136 (MQ=255)                                                                          cacaGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAgaa                                                           <  8:118408/139‑1 (MQ=255)                                                                          cacaGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAgaa                                                           >  7:118408/1‑139 (MQ=255)                                                                          cacaGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGa                                                 <  4:124786/149‑1 (MQ=255)                                                                                   ggcggcGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTc                                        >  7:100323/1‑149 (MQ=255)                                                                                    gcggcgGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATACCATGTGCCAGAAGGGAAGGCGAGTGATCGTCt                                       <  5:312522/149‑1 (MQ=255)                                                                                       gcggGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGGTCGTCTTcc                                    <  6:18743/149‑1 (MQ=255)                                                                                         ggggTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAgaagg                                                         <  1:185144/126‑1 (MQ=255)                                                                                         ggggTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAgaagg                                                         >  2:185144/1‑126 (MQ=255)                                                                                          gggTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCaaa                                 <  4:101717/149‑1 (MQ=255)                                                                                           ggTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCGCCGCTCGTCGCTGCGACTGGCAGCGTAAGTGACGTAAGGCCATGTGCCAGAAGGGAAGGGGAGTGATCGTCTTCCacac                                >  7:200434/1‑149 (MQ=255)                                                                                               ccAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGt                                         <  8:117041/136‑1 (MQ=255)                                                                                               ccAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGt                                         >  7:117041/1‑136 (MQ=255)                                                                                                 aGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCa                          <  5:25044/149‑1 (MQ=255)                                                                                                     aTAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGc                                                                >  1:133670/1‑107 (MQ=255)                                                                                                     aTAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGc                                                                <  2:133670/107‑1 (MQ=255)                                                                                                              ggagACATAACTGGTAATATATCATAACGAGCCATGCGGCGAGTGAAACACAAACGCTCGTCGCTGGGACTTTAAGCTTAAGTTAAGTAATGCAATGTTCCAGAAGGGAAGGAGAGATATCGACTTCCAAACCCGGCAGTAGcgcc                  >  2:250575/4‑146 (MQ=255)                                                                                                              g‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATc                                           <  3:215994/120‑1 (MQ=255)                                                                                                              g‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATc                                           >  4:215994/1‑120 (MQ=255)                                                                                                                       aCTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCAGTAGCACCACGAtgctgc        >  6:80614/1‑147 (MQ=255)                                                                                                                       aCTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCAGTAGCACCACGAtgctgc        <  5:80614/147‑1 (MQ=255)                                                                                                                        cTGGTAATAGCTCTTAACGAGCCATGCGACGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCAGTAGCACCACGATGCTGCCgg     <  2:86489/149‑1 (MQ=255)                                                                                                                        cTGGTAATAGCTCTTAACGAGCCATGCGACGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCAGTAGCACCACGATGCTGCCgg     <  2:90438/149‑1 (MQ=255)                                                                                                                           gTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCAGTAGCACCACGATGCTGCCGGTTc  <  8:100323/149‑1 (MQ=255)                                                                                                                                            gCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATc                                           >  3:197286/1‑91 (MQ=255)                                                                                                                                            gCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATc                                           <  4:197286/91‑1 (MQ=255)                                                                                                                                                    |                                                                                                                            GTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTC‑GGG‑CGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGACGTAATGCCATGTGCCAGAAGGGAAGGCGAGTGATCGTCTTCCAAACCCGGCAGTAGCACCACGATGCTGCCGGTTC  >  NC_002947/3420605‑3420873

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