Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,325,026	G→T	13.3%	P36P (CCC→CCA)	PP_2924 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,325,026	0	G	T	13.3%	97.3 / 4.8	37	P36P (CCC→CCA)	PP_2924	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/13);  new base T (3/2);  total (22/15)
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.79e-01

GTACAGGCAGAAGTTGGTCAGCAGTACCTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCAAGTGTTTGTGC  >  NC_002947/3324900‑3325172                                                                                                                               |                                                                                                                                                   gTACAGGCAGAAGTTGGTCAGCAGTACCTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCACGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCag                                                                                                                               >  7:272508/1‑148 (MQ=255)                 gTCAGCAGTACCTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTggc                                                                                                              <  5:240585/149‑1 (MQ=255)                  tCAGCAGTACCTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCag                                                                                                                               <  4:13360/131‑1 (MQ=255)                  tCAGCAGTACCTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCag                                                                                                                               >  3:13360/1‑131 (MQ=255)                      cagTACCTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGt                                                                                                         >  1:135095/1‑149 (MQ=255)                            cTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCaa                                                                                                   <  8:272508/149‑1 (MQ=255)                                 gCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATg                                                                                              >  3:148550/1‑149 (MQ=255)                                     tGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGcc                                                                                            >  7:256582/1‑147 (MQ=255)                                     tGAGACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGa                                                                                          <  8:50236/149‑1 (MQ=255)                                           gCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCTCGCGAATTGCGGTTTGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTTGCGGAGTCGTCAAGGATGCCGAGTTTgg                                                                                    >  1:57469/1‑149 (MQ=255)                                           gCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTgg                                                                                    >  8:136473/1‑149 (MQ=255)                                                    gTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGtt                                                                           <  7:136473/149‑1 (MQ=255)                                                            tataGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCAtggtg                                                                   >  2:172369/1‑149 (MQ=255)                                                                      aCAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGc                                                         >  1:209428/1‑149 (MQ=255)                                                                       cAGCTACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGc                                                         >  1:202724/1‑148 (MQ=255)                                                                        aGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGa                                                       <  2:202724/149‑1 (MQ=255)                                                                          aTACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATAGTGGCTGACCTGCGATa                                                     >  3:25543/1‑149 (MQ=255)                                                                           tACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATAc                                                    >  3:43343/1‑149 (MQ=255)                                                                              ccATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGa                                                 >  5:254855/1‑149 (MQ=255)                                                                                   ccGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTAtt                                            >  5:29459/1‑149 (MQ=255)                                                                                   ccGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTAt                                             >  3:21135/1‑148 (MQ=255)                                                                                    cGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTAtt                                            >  4:88949/1‑148 (MQ=255)                                                                                      ttGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCAt                                         <  8:256582/149‑1 (MQ=255)                                                                                         gtggCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGAGTCGGCAAGGATGCAGAGTTTGGAGATGAGTTGCATGGTGGCTGACCTGCGATAATGAGTATTCATa                                        >  2:192174/1‑147 (MQ=255)                                                                                            gCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTTGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTa                                   >  8:32039/1‑149 (MQ=255)                                                                                            gCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTTGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTa                                   >  8:36826/1‑149 (MQ=255)                                                                                              aCCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGa                                                 <  1:59547/133‑1 (MQ=255)                                                                                              aCCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGa                                                 >  2:59547/1‑133 (MQ=255)                                                                                              aCCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGGACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTAtt                                 >  2:250147/1‑149 (MQ=255)                                                                                                 cAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTg                              <  4:21135/149‑1 (MQ=255)                                                                                                   ggCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCCATACTGAGTATTCATACAGTATTTTgag                            <  1:172369/149‑1 (MQ=255)                                                                                                     ccGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGAt                          <  4:148550/149‑1 (MQ=255)                                                                                                                   cgAATTGCGCTTTGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCaa            <  7:32039/149‑1 (MQ=255)                                                                                                                   cgAATTGCGCTTTGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCaa            <  7:36826/149‑1 (MQ=255)                                                                                                                    gAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCAAg           >  4:160693/1‑149 (MQ=255)                                                                                                                     aaTTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCAAgt          <  3:88949/149‑1 (MQ=255)                                                                                                                          cgcTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCAAGTGTTtg     <  2:135095/149‑1 (MQ=255)                                                                                                                               gggCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCAAGTGTTTGTGc  <  3:160693/147‑1 (MQ=255)                                                                                                                               |                                                                                                                                                   GTACAGGCAGAAGTTGGTCAGCAGTACCTTGAGCAGTGATACGCAGCGGCCGTCCGGGGTATAGCTGTGACAGATACCCATGCCGTTGGTGGCACCCAGGCCGTCAGCGCCGCGCGAATTGCGCTTGGGCGCACCACTGCTGGCGCAGGAGGCGTCGTACTTGGCGGCGTCGGCAAGGATGCCGAGTTTGGCGATGAGTTGCATGGTGGCTGACCTGCGATACTGAGTATTCATACAGTATTTTGAGATGTGGCCGATTGCAAGTGTTTGTGC  >  NC_002947/3324900‑3325172

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