Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,989,261	A→C	15.0%	Y42D (TAT→GAT)	PP_2613 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,989,261	0	A	C	15.0%	85.6 / 11.7	40	Y42D (TAT→GAT)	PP_2613	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/17);  new base C (3/3);  total (20/20)
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.68e-01

CAAAGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCATCA  >  NC_002947/2989115‑2989408                                                                                                                                                   |                                                                                                                                                    cAAAGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGc                                                                                                                                                  <  2:582299/150‑1 (MQ=255)  aaaGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGcc                                                                                                                                                 >  1:468866/1‑150 (MQ=255)                 tGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCatc                                                                                                                                             >  2:715312/1‑137 (MQ=255)                 tGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCatc                                                                                                                                             <  1:715312/139‑3 (MQ=255)                     ggTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCat                                                                                                                                              >  2:615718/1‑133 (MQ=255)                     ggTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCat                                                                                                                                              <  1:615718/134‑2 (MQ=255)                         aTCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCAt                                                                                                                          <  1:370722/150‑1 (MQ=255)                                      ttcttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCAt                                                                                                                          >  1:440691/1‑137 (MQ=255)                                      ttcttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCAt                                                                                                                          <  2:440691/137‑1 (MQ=255)                                         ttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGt                                                                                                                             <  2:576542/131‑1 (MQ=255)                                         ttcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGt                                                                                                                             >  1:576542/1‑131 (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                                                                                          >  2:537848/1‑147 (MQ=255)                                                            gAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGt                                                                                          <  1:537848/147‑1 (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:50265/150‑1 (MQ=255)                                                                                tttATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtgat                                                                   <  2:749614/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)                                                                                                                                                 aTAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCATCa  <  1:456512/150‑1 (MQ=255)                                                                                                                                                 aTAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCATCa  >  1:288192/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAAAGGCCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCATCA  >  NC_002947/2989115‑2989408

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