Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,240	G→A	L21L (TTG→TTA)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,240	0	G	A	93.9%	96.1 / ‑2.5	33	L21L (TTG→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (12/19);  total (13/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 = 9.65e-01

GCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGA  >  NC_000913/4637996‑4638479                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                gcaccggtccgtgggtGGCGTGCGACAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCAAGCTGCGCAACGTGCAGGCAGGTGGCTTGGTGGCCGCGCTGCGCCTTCACAGTCACTTCAAATCGCTACAAAGTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCTTGTTCTGGAAAATGAATAGCATGTTCGGTGCAGTATTGCTGTTGTTAAtttcg                                                                                                                                                                                                                                           <  2:386085/236‑4 (MQ=255)                 ggCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCgg                                                                                                                                                                                                                            <  1:395924/250‑1 (MQ=255)                   cGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAg                                                                                                                                                                                                                         <  1:442144/251‑1 (MQ=255)                    gTTCGTCAGTGAGGTCGCCCGGTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGAATTGCTGTTGTTAATTCAGATAATGCTGATTCGGCAgg                                                                                                                                                                                                                        >  2:342371/1‑251 (MQ=255)                                                       aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGCGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGAtt                                                                                                                                                                                                                               >  2:271792/1‑209 (MQ=255)                                                       aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGAtt                                                                                                                                                                                                                               <  1:271792/209‑1 (MQ=255)                                                           cACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCGGATTCGGCAGGTGATTGTCGAACGTGCTGATTTCCGTAGCGATGTGGAAg                                                                                                                                                                                 <  1:363451/251‑1 (MQ=255)                                                                                   gggCGGGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCTATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCg                                                                                                                                                         <  2:76264/250‑1 (MQ=255)                                                                                    ggCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCgg                                                                                                                                                        >  1:48498/1‑251 (MQ=255)                                                                                                        aCTTCACCGTAACTTCACATAGCTTCAAATTCTTACCACATAGTCCTCGTATCCTGCTACCATTGCAAAGAAGAAGACTATGTTGAAATCCCCCCTGTTCTGGACAATGACTCCCCTGTTTGGTGCAGAATTGCTGTTGTTAATTCAGATATTGCTGATTCGGCAGGTGATTGTCGCAATTGCTGATTACCGTAGCGATGTGGAAGATTGGATTCGCCAAAGTAACCGCGGGACGGAAAAACGCGTTgggg                                                                                                                                    >  2:215124/1‑250 (MQ=255)                                                                                                                   ccctCACATAGCTACAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGTCTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGTGCCGCTCATCGCt                                                                                                                         <  2:121736/248‑1 (MQ=255)                                                                                                                      tCACATAGCTTCAAATTCTTCCCACATAGTATTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCACCTGTTCTGGAAAATGACTAGCATGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAACGTACCAGCGGGCCGCAAAAACTCGTTGGGCGGCTCCTCGCTAtt                                                                                                                      >  1:197310/1‑251 (MQ=255)                                                                                                                                     ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAACGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCtt                                                                                                       >  1:330225/1‑251 (MQ=255)                                                                                                                                                      tCGTAGCCTGCTGCCATTGCAAAGGATAATACTAGGTTGAAATCCCCCCTGTTTTGGAAAATGACTAGCCTGTTTTGTGCAGTATTGCTGTGGTTAAGTACGATAATGCTGATTCGGCAGGTGAGTGTCGAACGTTCTGATTCCCGTAGCGATGTGGAAGATGCGAGTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTTGGCCGCTCATCGCTATTCCTGTGACATAGCTTTATACGGTgtcggaag                                                                                       <  2:173488/250‑8 (MQ=255)                                                                                                                                                       cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATGCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAaga                                                                                      <  1:198564/250‑1 (MQ=255)                                                                                                                                                       cGTATCCTGCTGCAATTGCAAAGGAGAAGACAATGTTGAAATCCCCCCTTTTCTGGAAAATGACTAGCATGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGaa                                                                                        <  2:395042/248‑1 (MQ=255)                                                                                                                                                        gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAgg                                                                                    >  1:170804/1‑251 (MQ=255)                                                                                                                                                                  gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCACAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAcc                                                                            <  1:267605/249‑1 (MQ=255)                                                                                                                                                                           aaGGAGAAGATTATGTTGAAATCCCCCCGTTTCTGGAAAATGGCTTGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTTATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGAGCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGaa                                                                 <  2:22946/251‑1 (MQ=255)                                                                                                                                                                                 aagaCTATGTTGAAATCCCCCCTGTTCGGGAAAATGACTAGCCTGTTTGGTGCAGTATGTCTGTTGTTAATGCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGTAGCGGAAACGAAg                                                           <  2:317876/251‑1 (MQ=255)                                                                                                                                                                                    aCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGttt                                                        <  1:97737/251‑1 (MQ=255)                                                                                                                                                                                    aCTATGGTGAAATCCCCCCTGTTCTGGAAAATGACTAGTCTGTCTGGTGCAGGTTTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGGTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGGGCAGGAAGAGGATAGAACCGTGGAGCGGAAACGAAGttt                                                        <  1:156718/251‑1 (MQ=255)                                                                                                                                                                                              aaTCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACg                                                                                                 >  1:457114/1‑200 (MQ=255)                                                                                                                                                                                                  cccccTGTTCTGGAAAATGACTAGCCTGGTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGt                                          <  1:347421/251‑1 (MQ=255)                                                                                                                                                                                                            gggAAAATGCCTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTTATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTc                                <  2:361402/250‑1 (MQ=255)                                                                                                                                                                                                                 aaTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCCAAGTACCAGCGGGCCGCCAAAACTCGTTGGGCCGCTAATCGCTATTCCTGGGACCGGGTTTTATAGGGTGACGGAAGAGGATAAAACCGTGGAGAGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGCCAttgg                           >  2:187391/1‑250 (MQ=255)                                                                                                                                                                                                                  aTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTc                                >  1:350582/1‑245 (MQ=255)                                                                                                                                                                                                                  aTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTc                                <  2:350582/245‑1 (MQ=255)                                                                                                                                                                                                                          ccTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttgatg                  <  1:342371/251‑1 (MQ=255)                                                                                                                                                                                                                            tgttgtgTGCTGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGATCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttgatg                  <  2:68120/243‑1 (MQ=255)                                                                                                                                                                                                                             gTTTGGTGCAGTATTGCTGATGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGTTTACCGTAGTGATGTGGAAGATGCGATTCGACAAAGTACCAGCGGACCGCAAAAACGCGGTGGGACGCTGATCGCTATTACTGTGACCGGGACTTATCCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTATGGTTACCTGAgaccattatggtgaatg                  >  2:361067/1‑232 (MQ=255)                                                                                                                                                                                                                              tttGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCa               >  1:137343/1‑250 (MQ=255)                                                                                                                                                                                                                                           tGCTGATGTTGATTCCGATAATGCTGAGTCGGCAAGAGATTGTCGTACGTGCCGATTACCGTAGCGATGTGGACGATGCGATTAGACAAAGTAACAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTCTTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGCTAAAACCGTGGAGCGGAACCGAAGTTTTATCCTTTTTTGGTTACCTGAGTCATTGATGGTTGATGGCATTCAGaaaggtga  >  2:236018/1‑244 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                GCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGA  >  NC_000913/4637996‑4638479

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