Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	856,914	A→G	N318D (AAC→GAC)	ybiT →	putative ATP‑binding protein YbiT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	856,914	0	A	G	96.9%	93.5 / NA	32	N318D (AAC→GAC)	ybiT	putative ATP‑binding protein YbiT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (15/16);  minor base T (0/1);  total (15/17)
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

TACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAA  >  NC_000913/856671‑857150                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             tACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCt                                                                                                                                                                                                                                       >  1:33821/1‑251 (MQ=255)            cATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCtggaagtggaag                                                                                                                                                                                                                            <  2:239061/251‑1 (MQ=255)               gACGGCGGCGACCCAGGCGCGTGAACGTCGGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTc                                                                                                                                                                                                                         <  2:33821/251‑1 (MQ=255)                aCGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACTCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCt                                                                                                                                                                                                                        <  1:152727/251‑1 (MQ=255)                        gaccctggcgctggaccGTCTGCTGGCCGATAACGCAAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCGCGCCAGGCAACTTCGCGCGCGCGCCAGGTTGATAAAATCAAACTGGAAGAGTTGAAATCCTCCAGCCGTCAGAACCCGTTCTTCCGTTTGGAACAGGATAAGAAACTGTTCCGTGACGCGCGGTAAGAGTAAGGTCTGACCAAAg                                                                                                                                                                                                                <  2:197488/235‑1 (MQ=255)                        gACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCGAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAg                                                                                                                                                                                                                <  1:265994/251‑1 (MQ=255)                         aCCCAGGCGCGTTAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGATTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCAGCCGTTGTTAACAGGATAAGAAACTGTTCCGTGACGCGCTGTAAGTTGAAGGTCTGACCAAAgg                                                                                                                                                                                                               <  2:208350/251‑1 (MQ=255)                              cgacgcgtttccgatcttgGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTTCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTGCATCCGGTTTGAACAGGATAAGAAACTGTTCCGGGACGCGCTGGAAGTGGAAGGTCTTACCAAAGGGTTTg                                                                                                                                                                                                          <  2:313199/234‑1 (MQ=255)                                             gctgGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTGCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTg                                                                                                                                                                                            <  1:83743/250‑1 (MQ=255)                                               tgGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGCGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTg                                                                                                                                                                                                          >  1:313199/1‑234 (MQ=255)                                                           ccAAGATGAAAGCGCAGATTGCTGCGTTGCAATCTTTCGTTAGCCGCTTTAGAGCCAACGCCTCGCAATGTCGCCAGGCAAATTCACGCGCGCGCCAGATTGATAAAATCAAACTGGATGAGGTGAAAGCCACCAGCCGTCAGAAACCGTTCATCCGTTTTGAACAGGAGAAGAAACTGTTCTGTGCCGCGCTGCAAGTGGAAGGTATGACCAAAGGGTGTGATAACGGTCCGCTGTTTAAAAATCTCAAc                                                                                                                                                                             >  2:148297/1‑251 (MQ=255)                                                             aagaagAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCACAACct                                                                                                                                                                           <  1:85871/251‑1 (MQ=255)                                                                     aGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTgg                                                                                                                                                                   >  2:229717/1‑251 (MQ=255)                                                                                cTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGa                                                                                                                                                         >  1:131134/1‑250 (MQ=255)                                                                                  gatttgcatgcttccgttgtcctcTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACGACTAGCGCGCGCCAGATGGATAAAAGCAAACTGGAAGAGGTGAAAGCCGCCAACCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTTCGTGACTCGCTGGATGTGGAAGGTCTGAACAAAGGGTTTGATAACGTTCCGCTTTTTAAAAATCTCACCCTACTGCTAGAAGTGGGTaaaa                                                                                                                                                       <  2:226891/227‑1 (MQ=255)                                                                                                cGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAACATCAAACTGGAAGAGGAGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACCGGATCAGAAACTGTTCCGTGACGAGCTGGAAGGGGAAGGACTGCCCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATATCAACATGCTGCTGGAAGTGGGTGACAACCTGGCGGTAcg                                                                                                                                          >  2:294188/1‑248 (MQ=255)                                                                                                cGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAACTCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTgg                                                                                                                                        >  2:93132/1‑251 (MQ=255)                                                                                                                gcCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCgg                                                                                                                        >  1:250281/1‑251 (MQ=255)                                                                                                                           gAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAACGCCTCCAGCCGTCAGAAACCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGAAGCGATGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGACGTGGGTGAAGAACTGTCGGTCCTGGGTATCAACGGCGTCGGTAAATCAACgc                                                                                                             >  2:111999/1‑251 (MQ=255)                                                                                                                                            cTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTAAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGACGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGAGGGTGAAAAACTGGCGGTACTGGGTAACAACGGCGTCGGTAAATCAACGCTGATGAAACCGCtggtg                                                                                            >  1:238701/1‑251 (MQ=255)                                                                                                                                                                gATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAg                                                                        <  1:98023/251‑1 (MQ=255)                                                                                                                                                                     aaTCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGc                                                                    >  1:167008/1‑250 (MQ=255)                                                                                                                                                                              ggAAGAGGTGAAATCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTTACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGTTACCAGCGGCGTCGGTAAATCAACGCTGCGGAAAACGCTGGTGGGCGAGCTGCAACCTGACAGCGGCACCGTaaa                                                            <  2:131134/249‑1 (MQ=255)                                                                                                                                                                              ggAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGCACCGTaaaa                                                           <  1:69210/250‑1 (MQ=255)                                                                                                                                                                                 agagGTGAACGACTCCAACCGCCAGAAACCCTTCATCCGTTTTGAACAGGCTAACAAAATGTTACGTGACGCGCAGGAAGTGAAAGGACTGACCAAAAGGTTTGATCACGGTCCGATGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGCGACAAACTGGCGGTACTGGATACAAACGGCATCGGTACATCAACGCTGCTGAAAACGCTGGTGCGCAATCTGAAACCGGACAGTGGCACCGCAAAATgtg                                                       >  2:127906/1‑249 (MQ=255)                                                                                                                                                                                       gAAAGCCTCCAGCCGACAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCAGGAAGTGGACGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAACTCTCAACCTGCTGCTGGAAGTGGGTGAAAACCTGGAGGGACTGGGTACCAAAGGCGTCGGTAAACCAACGCTGCTGAAAACGCTGGGGGGCGAGCAGCAACCGGACCGCGGCACCGTAAACTGGtacgagg                                                 >  2:186715/1‑245 (MQ=255)                                                                                                                                                                                                cAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACACGCGCa                                        <  1:93132/251‑1 (MQ=255)                                                                                                                                                                                                 aGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGGTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCAt                                       <  2:102185/251‑1 (MQ=255)                                                                                                                                                                                                   ccGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTg                                     <  1:88752/251‑1 (MQ=255)                                                                                                                                                                                                              cGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGGGGGCGATCTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCt                          >  1:50599/1‑251 (MQ=255)                                                                                                                                                                                                                         tttGAACAGGATAAGAAACTGTTCCGTGACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACg                <  1:106647/250‑1 (MQ=255)                                                                                                                                                                                                                                       gAAACTGTTCCGTGACGCGATGGACGTGGCAGGTCTGACCAAAGAGTTTGATAACGGTCCGCTGTTTAAAAATCTAAACCTGCTGCTGGAAATGGGTGAAAAACTGGCGGTACTGGGTACCCACGGCGTCGGTAAATCAAAGCTGCTGAAAACGATGGGGGGCGATCTGCAAACGGCCAGCGGCACCGTAAAATGGTCTGAGAACGCGGGCATTGGTTTCTATGCTCAGGAACACGAATATGAGTTTGaa  >  2:313177/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             TACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGTCTGACCAAAGGGTTTGATAACGGTCCGCTGTTTAAAAATCTCAACCTGCTGCTGGAAGTGGGTGAAAAACTGGCGGTACTGGGTACCAACGGCGTCGGTAAATCAACGCTGCTGAAAACGCTGGTGGGCGATCTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAA  >  NC_000913/856671‑857150

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