Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	7,939	G→A	F7F (TTC→TTT)	yaaJ ←	putative transporter YaaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	7,939	0	G	A	97.0%	101.2 / ‑6.6	33	F7F (TTC→TTT)	yaaJ	putative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/16);  minor base T (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

GCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCC‑TGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCAT  >  NC_000913/7712‑8184                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        cccagccaaaCGCCGGCAAGGTTGCCGCTCCCCACGCGCGCCGCAAGACTGGTCCACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACGGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGgcatgca                                                                                                                                                                                                                                  <  2:443645/243‑1 (MQ=255)                                      gcgcCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtgtg                                                                                                                                                                                             <  1:125236/251‑1 (MQ=255)                                             aaGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGgca                                                                                                                                                                                                                                      >  2:269606/1‑203 (MQ=255)                                             aaGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGgca                                                                                                                                                                                                                                      <  1:269606/203‑1 (MQ=255)                                                 cTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTAcc                                                                                                                                                                                  >  1:30257/1‑251 (MQ=255)                                                   ggTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                                                                                                                                                                 <  2:30257/250‑1 (MQ=255)                                                                 gAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGCAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGCTGACTTTCCTGGGTAAAAGTTACCA‑CTTGTTTaagacgt                                                                                                                                                                  >  2:481563/1‑247 (MQ=255)                                                                  aaaTGAGGTTAAACCGCCTGCCTGTGGATGAATGCTATTTATAAGACTTTTGCCAAACTGGCGGATGTAGCGAGACTGCACAAATACGGTGCGAAAAGTGGACCAACAACCTGCGCCGAAGAGCAGGTAACTCATTACCGCTCCCCCAAGGACACTGTTAATTAGGGAGACAAAATCAGGCAAGCACATCCCTCTTATTGCCGGGTGCGACGACTTGCCTGTGGAAACGTTACCACATT‑TTTAAGACGtat                                                                                                                                                                >  2:455268/1‑251 (MQ=255)                                                                       aGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGCGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCCTGCATATCCCTCTTATTGCCGGTCGCGCTGACTTTCCGGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATAcc                                                                                                                                                           >  2:221045/1‑250 (MQ=255)                                                                             aaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACg                                                                                                                                                      >  2:134037/1‑250 (MQ=255)                                                                              aCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAg                                                                                                                                                    <  1:349499/251‑1 (MQ=255)                                                                                  ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGAGACCCACAAAACTGCGACGAAGAGTAGGAAAATCATTACCGATAACCAAAGGACGCTGCTAATAAAGGAGAAAAAATCTGGCATGACTATCCCTCTTATTGCCGGTCGCGATGACTTTCATGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTCCGAGGTAc                                                                                                                                                >  1:65866/1‑251 (MQ=255)                                                                                           ggATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAAACTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATAAACGCTACGAGGTACTTGATaaat                                                                                                                                       >  1:217296/1‑249 (MQ=255)                                                                                                      aTTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGc                                                                                                                             <  2:185128/250‑1 (MQ=255)                                                                                                           tAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTTCTTGATAACTTCTGCGTAGCATACa                                                                                                                        <  2:149829/250‑1 (MQ=255)                                                                                                                          aCTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATaaa                                                                                                        <  1:190682/251‑1 (MQ=255)                                                                                                                                   tGTAGCGAAACTGCACAAATCCGGTGCGAAAAGCGACCCAACAACCTGCGCCGAAGAGCAGGGAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATggcggg                                                                                               >  1:147383/1‑251 (MQ=255)                                                                                                                                         gAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATAGCGGGCGatat                                                                                         >  1:481537/1‑251 (MQ=255)                                                                                                                                          aaaCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATc                                                                                        <  2:271435/251‑1 (MQ=255)                                                                                                                                          aaaCAGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATc                                                                                        <  1:211637/251‑1 (MQ=255)                                                                                                                                               gCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAAc                                                                                     <  1:179556/249‑1 (MQ=255)                                                                                                                                                  caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAgg                                                                                >  1:278542/1‑250 (MQ=255)                                                                                                                                                                  agtaaaCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAAGCCGAAAc                                                                <  2:347794/247‑1 (MQ=255)                                                                                                                                                                    tGAACCAACAACCTGCGCCGCAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATCAGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGGGGCAGAAATACGAAAc                                                                >  2:399147/1‑249 (MQ=255)                                                                                                                                                                                cTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCC‑TGGCg                                                 >  1:331967/1‑251 (MQ=255)                                                                                                                                                                                     ccGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCC‑TGGCGATAAc                                            <  2:278542/251‑1 (MQ=255)                                                                                                                                                                                             cAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTa                                                                                                                               >  2:105180/1‑161 (MQ=255)                                                                                                                                                                                             cAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCGGCGTa                                                                                                                               <  1:105180/161‑1 (MQ=255)                                                                                                                                                                                              aGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCC‑TGGCGATAACCGTCTTGTc                                   <  1:403092/251‑1 (MQ=255)                                                                                                                                                                                                   actCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGGTATCAACGCAGTGGCAGCAATACGAAACAGTCTCGCC‑GGGCGATAACCCGCTTGTcgacgg                              >  2:372635/3‑251 (MQ=255)                                                                                                                                                                                                      cATTACCGATCCCCAAAGGACGCTGTTAATAACGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGAAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGGGATATACACGCATGGAAAGAAATACGGAACAATCTCGCC‑GGGCGGTAACCAGCTTTTCGGCGGtgg                           >  2:350621/1‑249 (MQ=255)                                                                                                                                                                                                          aCCGAACCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGGCTCCCC‑TGGCGATAACCGTCTTGTCGGCGGTTGCGCt                       <  1:369741/251‑1 (MQ=255)                                                                                                                                                                                                                  ccAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCC‑TGGCGATAACCGTCTTGTCGGCGGGTTCGCTgtcggtgt               >  1:90874/1‑244 (MQ=255)                                                                                                                                                                                                                                taataaAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGGTTAAGAAGAATATACGCTACGAGGAACTTGGTAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTAAGGAATCCGAAAAAAGCGCGCAATGACG‑TAACCGTCTTGTCGGCGGTGGGGCTGAAGTTGCgtggtggagagaa  >  1:290860/1‑243 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        GCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA‑ATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCC‑TGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCAT  >  NC_000913/7712‑8184

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