Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,998,329:1	+C	8.3%	coding (54/903 nt)	rfbD ←	dTDP‑4‑dehydrorhamnose reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,998,329	1	.	C	8.3%	178.7 / 10.5	48	coding (54/903 nt)	rfbD	dTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base . (25/19);  new base C (2/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGCAGCGTTGACGATCACGTCCGGCGCCAGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAGTACCAAGC  >  NC_002947/1998191‑1998474                                                                                                                                            |                                                                                                                                                  ggCAGCAGCGTTGACGATCACGTCCGGCGCCAGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCt                                                                                                                                         >  8:210286/1‑149 (MQ=255)   cagcagCGTTGACGATCACGTCCGGCGCCAGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTcc                                                                                                                                       <  7:188191/149‑1 (MQ=255)                       tCCGGCGCCAGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGc                                                                                                                                   >  3:305351/1‑133 (MQ=255)                       tCCGGCGCCAGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGc                                                                                                                                   <  4:305351/133‑1 (MQ=255)                           gcgcCAGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTa                                                                                                                                            <  3:257144/120‑1 (MQ=255)                                aGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAAcag                                                                                                          >  4:159187/1‑149 (MQ=255)                                         aTCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGCCCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGAttt                                                                                                  >  2:110730/1‑148 (MQ=255)                                                cGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCtt                                                                                          <  1:323777/149‑1 (MQ=255)                                                  gCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTa                                                                                        >  3:264588/1‑149 (MQ=255)                                                    cAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTAtt                                                                                      >  8:327833/1‑149 (MQ=255)                                                          gTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGcccc                                                                                >  2:358019/1‑149 (MQ=255)                                                                 ttGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTAttgt                                                                         >  7:48624/1‑149 (MQ=255)                                                                 ttGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTAttgt                                                                         >  5:257346/1‑149 (MQ=255)                                                                     aCAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTtgtg                                                                     >  8:164450/1‑149 (MQ=255)                                                                      cAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGc                                                                    <  7:210286/149‑1 (MQ=255)                                                                       aaGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGcc                                                                   >  4:168340/1‑149 (MQ=255)                                                                         gTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTc                                                                                                <  2:216560/118‑1 (MQ=255)                                                                         gTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTc                                                                                                >  1:216560/1‑118 (MQ=255)                                                                             ccACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCTGGATTTTCACGCTTTATTTGcccc                                                                                <  2:56367/130‑1 (MQ=255)                                                                             ccACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGcccc                                                                                >  1:56367/1‑130 (MQ=255)                                                                                caTAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCa                                                          <  4:264588/149‑1 (MQ=255)                                                                                 aTAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTa                                                                             <  6:219857/129‑1 (MQ=255)                                                                                 aTAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTa                                                                             >  5:219857/1‑129 (MQ=255)                                                                                      cccTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAg                                                    <  5:146372/149‑1 (MQ=255)                                                                                       ccTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATtgttgt                                                                      >  3:238903/1‑130 (MQ=255)                                                                                       ccTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATtgttgt                                                                      <  4:238903/130‑1 (MQ=255)                                                                                                ccTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTa                                                     >  6:232579/1‑138 (MQ=255)                                                                                                ccTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAg                                                    <  8:329457/139‑1 (MQ=255)                                                                                                ccTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAg                                                    >  7:329457/1‑139 (MQ=255)                                                                                                 cTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCAt                                         >  8:262172/1‑149 (MQ=255)                                                                                                   gACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGa                                       >  7:83777/1‑149 (MQ=255)                                                                                                            gcgcAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTa                              <  3:159187/149‑1 (MQ=255)                                                                                                             cgcAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTAc                             <  7:164450/149‑1 (MQ=255)                                                                                                                    aCCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCAttt                       <  7:327833/148‑1 (MQ=255)                                                                                                                      cTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTg                    >  8:269983/1‑149 (MQ=255)                                                                                                                      cTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTg                    <  1:358019/149‑1 (MQ=255)                                                                                                                           tacagcgcgcCTTGGCGCCGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCAt                         <  7:233432/134‑1 (MQ=255)                                                                                                                           tacagcgcgcCTTGGCGCCGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCAt                         >  8:233432/7‑140 (MQ=255)                                                                                                                           ccAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGCCCCCTTTTACCCATTTGTGGTTgg               >  3:71975/1‑149 (MQ=255)                                                                                                                             cagcgcgcCTTGGCGCCGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGa                                       <  2:130565/120‑1 (MQ=255)                                                                                                                             cagcgcgcCTTGGCGCCGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGa                                       >  1:130565/5‑124 (MQ=255)                                                                                                                             aGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGcc             <  3:168340/149‑1 (MQ=255)                                                                                                                                ggcgCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCGCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAg          >  1:113033/1‑149 (MQ=255)                                                                                                                                 gcgcCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAGt         <  2:113033/149‑1 (MQ=255)                                                                                                                                  cgcCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAGTa        >  2:245052/1‑149 (MQ=255)                                                                                                                                       aGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAGTACCAAg   >  6:91755/1‑149 (MQ=255)                                                                                                                                       aGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAGTACCAAg   <  1:110730/149‑1 (MQ=255)                                                                                                                                        ggCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAGTACCAAGc  >  1:182789/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                  GGCAGCAGCGTTGACGATCACGTCCGGCGCCAGCTGACGAATCGTAGCGGCCAAGCCGTCCAGGTTGGACAAGTCGCCACATAAGCCCTCGGCCCCCTGACGATCCAGCGCAATGACCTCACCCAGCGGCGCCAAGGCG‑CGCTGTAGCTCCCAGCCTACTTGCCCGTTTTTCCCCAACAGCAGGATTTTCACGCTTTATTTGCCCCGTATTGTTGTGCCACCCAGTCACGGTAGCTGCCGTCCATGACACCTTTTACCCATTTCTGGTTGGCCAAGTACCAAGC  >  NC_002947/1998191‑1998474

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