Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	116,569	A→G	G177G (GGT→GGC)	hofB ←	T2SSE family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	116,569	0	A	G	97.0%	104.2 / NA	33	G177G (GGT→GGC)	hofB	T2SSE family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/16);  minor base C (1/0);  total (17/16)
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

CAACTGGATCTTCGACGCTACAAATATTAATGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCACCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAGGCATTGTCCGCTGGTTCGATATGAATATCAGACGCGCG  >  NC_000913/116340‑116814                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cAACTGGATCTTCGACGCTACAAATATTAATGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTc                                                                                                                                                                                                                                  <  1:309147/251‑1 (MQ=255)      ggATCTTCGACGCTACAAATATTAATGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCaa                                                                                                                                                                                                                              <  1:184882/250‑1 (MQ=255)        aTCTTCGACGCTACAAATATTAATGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATAc                                                                                                                                                                                                                           >  1:130888/1‑251 (MQ=255)            tCGACGCTACAAATATTAATGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATa                                                                                                                                                                                                                            <  1:101612/246‑1 (MQ=255)            tCGAAGCTACAAATATTAATGTCAGCGGTATTCAGCTTTTGCAGGGCACTATACAGCGTGACCGTTTTGCCGCTGCCGGTAGGGACAGTTACCAGCACCCGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATa                                                                                                                                                                                                                            <  2:193190/246‑1 (MQ=255)                           ttAATGTCAGCGGTCTTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCACGGCATGAGCAACGTCCGCCAGTTGTAACGGCTGCATTACAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGCGACGGCGTTTcc                                                                                                                                                                                                        >  2:344338/1‑251 (MQ=255)                           ttAATGTCAGCGGTCTTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTCACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCAC‑TGGTGTAACCACCCTAATTCCAACCTTTTCCCCCCCCCGCCTGGGTAAGGTCGCATTACCAAATGAGACGCCGttccc                                                                                                                                                                                                        >  2:344334/1‑248 (MQ=255)                             aaTGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCt                                                                                                                                                                                                       <  1:114141/250‑1 (MQ=255)                             aaTGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTg                                                                                                                                                                                                      >  2:117152/1‑251 (MQ=255)                                  cAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCCCCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGt                                                                                                                                                                                                 <  1:32331/251‑1 (MQ=255)                                            cAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTg                                                                                                                                                                                        <  2:299233/250‑1 (MQ=255)                                               cTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAAtt                                                                                                                                                                                    >  2:104739/1‑251 (MQ=255)                                                ttttGCAGGGAACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTACCGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTg                                                                                                                                                                                   >  2:160280/1‑251 (MQ=255)                                                    gCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGccc                                                                                                                                                                                <  2:130888/250‑1 (MQ=255)                                                           aCTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCg                                                                                                                                                                        >  1:286234/1‑251 (MQ=255)                                                              ataAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGc                                                                                                                                                                      >  1:158490/1‑250 (MQ=255)                                                                         ccGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTcc                                                                                                                                                          <  2:97881/251‑1 (MQ=255)                                                                               tGCCGCTGCCGGTAGGGCCAGTTACCAGAACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGCCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCGTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGAGAGGTTCCGCAat                                                                                                                                                     <  2:284736/250‑1 (MQ=255)                                                                                      gCCGGTTGGGCCAGTTACCCGCACCTGTCCTTGTGGTTGTTGCATGTCATGATCAAAGTCCGCCAGTTGTATCTGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCTCCTGCTGTAACATCCTTAAGACCT‑CTTTTTCGCCCCCCCGACATGGTAAGGTCGCAATACGAAATGAGTCGGCGTTTCCTGCCATTACGACAGTGAATTGGCCGTCCTGCTGCTGGCGATGTTCCGCAATATCCAGGt                                                                                                                                             <  1:360682/251‑1 (MQ=255)                                                                                                ccAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCCGCAc                                                                                                                                   >  2:107325/1‑251 (MQ=255)                                                                                                       ccAGCACGAGGCCCTGTGGTTGTTGCAAGGCAGGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTGCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAAt                                                                                                                            <  2:286234/251‑1 (MQ=255)                                                                                                                 tCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGtt                                                                                                                   >  2:95038/1‑250 (MQ=255)                                                                                                                     tgtgGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGc                                                                                                              >  1:325856/1‑251 (MQ=255)                                                                                                                                       tGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTg                                                                                             >  1:17651/1‑250 (MQ=255)                                                                                                                                                            aaCGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCaa                                                                       <  1:183086/251‑1 (MQ=255)                                                                                                                                                               ggCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATAc                                                                    >  1:33278/1‑251 (MQ=255)                                                                                                                                                                                    aCATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGAt                                               <  1:197062/251‑1 (MQ=255)                                                                                                                                                                                             gCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTACGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAggt                                      >  1:99108/1‑250 (MQ=255)                                                                                                                                                                                               cTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAAGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAGGCAt                                    >  1:246810/1‑251 (MQ=255)                                                                                                                                                                                                                acaacCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAGGCATTGTCCGCTGGTTCGata                   <  1:130763/251‑1 (MQ=255)                                                                                                                                                                                                                  aacCTTAATACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAGGCATTGTCCGCTGGTTCGATATg                 <  1:16967/251‑1 (MQ=255)                                                                                                                                                                                                                        aaTACCA‑CCTTTTCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAGGCATTGTCCGCTGGTTCGATATGAATATc           >  2:234134/1‑251 (MQ=255)                                                                                                                                                                                                                                  ttttCGCCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAGGCATTGTCCGCTGGTTCGATATGAATATCAGAcgcgcg  >  2:66079/1‑251 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CAACTGGATCTTCGACGCTACAAATATTAATGTCAGCGGTATTCAGCTTTTGCAGGGCACTATAAAGCGTGACCGTTTTGCCGCTGCCGGTAGGGCCAGTTACCAGCACCAGTCCCTGTGGTTGTTGCAAGGCATGAGCAAAGTCCGCCAGTTGTAACGGCTGCATTCCAAGCGTGTTGACATCCAGTGCCTGACCCACCTGCTGTAACAACCTTAATACCA‑CCTTTTCACCACCCCGACATGGTAAGGTCGCAATACGAAATGAGACGGCGTTTCCTGCCAGTTCGACAGTGAATTGCCCGTCCTGCGGCAGGCGATGTTCCGCAATATCCAGGTTTCCCAGCACTTTTAATCTGGCGGTTAATGCGACTCCGGCATCCGGTGAAACATCCGGTAAAGGATGCAATACGCCGTCGATACGCAAGCGGATGCGGTAGGCATTGTCCGCTGGTTCGATATGAATATCAGACGCGCG  >  NC_000913/116340‑116814

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