Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,989,258	G→T	15.0%	R43R (CGG→AGG)	PP_2613 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,989,258	0	G	T	15.0%	101.5 / 8.8	40	R43R (CGG→AGG)	PP_2613	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base T (3/3);  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 = 4.89e-01

ACCAGCAAAGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAG  >  NC_002947/2989110‑2989399                                                                                                                                                     |                                                                                                                                              accaGCAAAGGCCCAAAATCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGa                                                                                                                                              >  2:747404/1‑150 (MQ=255)  ccaGCAAAGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGAt                                                                                                                                             >  1:582299/1‑150 (MQ=255)      cAAAGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGc                                                                                                                                         <  2:582299/150‑1 (MQ=255)       aaaGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGcc                                                                                                                                        >  1:468866/1‑150 (MQ=255)                      tGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCatc                                                                                                                                    <  1:715312/139‑3 (MQ=255)                      tGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCatc                                                                                                                                    >  2:715312/1‑137 (MQ=255)                          ggTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCat                                                                                                                                     <  1:615718/134‑2 (MQ=255)                          ggTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCat                                                                                                                                     >  2:615718/1‑133 (MQ=255)                              aTCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCAt                                                                                                                 <  1:370722/150‑1 (MQ=255)                                           ttcttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCAt                                                                                                                 <  2:440691/137‑1 (MQ=255)                                           ttcttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCAt                                                                                                                 >  1:440691/1‑137 (MQ=255)                                              ttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGt                                                                                                                    >  1:576542/1‑131 (MQ=255)                                              ttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGt                                                                                                                    <  2:576542/131‑1 (MQ=255)                                              ttcCATAACCATCCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGTGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCgg                                                                                                 >  2:371929/1‑150 (MQ=255)                                                ccATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCatc                                                                                                                                    <  2:455885/113‑3 (MQ=255)                                                ccATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCatc                                                                                                                                    >  1:455885/1‑111 (MQ=255)                                                     aCCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTtgt                                                                                          <  1:710340/150‑1 (MQ=255)                                                        aTGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAg                                                                                       >  2:554589/1‑150 (MQ=255)                                                         tGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGc                                                                                      >  2:98554/1‑150 (MQ=255)                                                            cTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGaa                                                                                   <  2:551336/150‑1 (MQ=255)                                                                 gAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGt                                                                                 <  1:537848/147‑1 (MQ=255)                                                                 gAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGt                                                                                 >  2:537848/1‑147 (MQ=255)                                                                   gTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCg                                                                            >  1:50800/1‑150 (MQ=255)                                                                                  gCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCGATTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGt                                                             >  1:50265/1‑150 (MQ=255)                                                                                    gTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtga                                                           >  1:749614/1‑150 (MQ=255)                                                                                    gTTTATAAGCTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtga                                                           <  1:747404/150‑1 (MQ=255)                                                                                     tttATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtgat                                                          <  2:749614/150‑1 (MQ=255)                                                                                     tttATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtgat                                                          <  2:50265/150‑1 (MQ=255)                                                                                             gTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAggt                                                  >  1:636184/1‑150 (MQ=255)                                                                                              tCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAggtg                                                 <  2:468866/150‑1 (MQ=255)                                                                                               cGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGTGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAggtgg                                                <  1:371929/150‑1 (MQ=255)                                                                                               cGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAggtgg                                                <  1:554589/150‑1 (MQ=255)                                                                                                     ccGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTg                                          >  2:456512/1‑150 (MQ=255)                                                                                                           gggCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTg                                    <  2:636184/150‑1 (MQ=255)                                                                                                                    cAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGAt                           >  1:595016/1‑150 (MQ=255)                                                                                                                        aaCCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACaa                       >  2:118317/1‑149 (MQ=255)                                                                                                                             gTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCa                  >  2:258266/1‑150 (MQ=255)                                                                                                                               tCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAAt                >  2:393909/1‑150 (MQ=255)                                                                                                                                       cTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAgg        <  2:595016/150‑1 (MQ=255)                                                                                                                                             gATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAg  <  1:393909/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              ACCAGCAAAGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAG  >  NC_002947/2989110‑2989399

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