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	97.1%	112.7 / ‑3.8	35	L21L (TTG→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (15/19);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

GGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAAC  >  NC_000913/4638001‑4638474                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           ggTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAAt                                                                                                                                                                                                                                 <  1:7136/251‑1 (MQ=255)               gTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAgg                                                                                                                                                                                                                   <  1:236134/251‑1 (MQ=255)                           ggTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc                                                                                                                                                                                                       <  1:333722/251‑1 (MQ=255)                              cGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg                                                                                                                                                                                                    <  1:53348/251‑1 (MQ=255)                                                 gAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGAtgt                                                                                                                                                                                 >  1:194425/1‑251 (MQ=255)                                                                     gCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCaaa                                                                                                                                                             <  1:95935/251‑1 (MQ=255)                                                                               ggCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCa                                                                                                                                                               >  1:7534/1‑239 (MQ=255)                                                                               ggCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCa                                                                                                                                                               <  2:7534/239‑1 (MQ=255)                                                                                            cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGa                                                                                                                                                                                                                                     >  1:88097/1‑156 (MQ=255)                                                                                            cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCGTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGa                                                                                                                                                                                                                                     <  2:88097/156‑1 (MQ=255)                                                                                              cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGt                                                                                                                                    >  1:4598/1‑251 (MQ=255)                                                                                                 cGACTTCACCGTCACTTCAAATAGCTTCAAATTCTTCCCACATCGTCTTCGTATCCTGCTACCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGCTTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTAACAGCGGGCCGCAAAAACTCGTTgg                                                                                                                                 >  2:404254/1‑251 (MQ=255)                                                                                                   gcTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAACTGACTAGCCTGTTTGGTGCAATATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGCACCAGCGGGCCGCAAAAACTCGTTGGGc                                                                                                                               >  2:150039/2‑251 (MQ=255)                                                                                                      tCACCGTCACTTCACATAGCTACAAATTCTCCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTATCCTGTTTTGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCTGGTGATTGGCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGcc                                                                                                                              <  1:373001/249‑1 (MQ=255)                                                                                                       cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCAATGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCt                                                                                                                           >  1:141302/1‑251 (MQ=255)                                                                                                                     aTAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCATTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGCTTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCtg                                                                                                             >  1:178436/1‑251 (MQ=255)                                                                                                                                       cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGATCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGTTGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTAACCGAGCTTTATACgg                                                                                           >  1:323740/1‑251 (MQ=255)                                                                                                                                          aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAACCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTTCTGTGACCGAGCTTTATACGGTGc                                                                                        <  1:343585/251‑1 (MQ=255)                                                                                                                                                           ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAcc                                                                       >  2:389285/1‑251 (MQ=255)                                                                                                                                                                 ttGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                                   >  2:5429/1‑249 (MQ=255)                                                                                                                                                                   gCAAAGGAGAAGAATATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTTGAGCg                                                               <  2:97232/251‑1 (MQ=255)                                                                                                                                                                         gagaAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCg                                                                                                                                                                  <  2:250150/146‑1 (MQ=255)                                                                                                                                                                         gagaAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCg                                                                                                                                                                  >  1:250150/1‑146 (MQ=255)                                                                                                                                                                                        aaaTCACCCATGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATACTGATGCTTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGCACCAACGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGAACGAGATTTCTACGGGGCAGGAACAGGATAAAACCGGGGAGCGGGAACGAAGTTTTATCCAttt                                          >  2:386300/1‑251 (MQ=255)                                                                                                                                                                                                  tGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAACCTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTTCAGGAAGAGGATAAAACCGTGGAGCGGAAACGaa                                                       <  1:62399/228‑1 (MQ=255)                                                                                                                                                                                                  tGTTATGGAAACTGACTTGCCTGTTTGGTGCCGGATTGCTGTTGTTAATTCCGCTAATGCTGCTTCGGCAGGGGATTGTCGAACGTGCTGATTACCGTGGCGATGTGGAAGATGCGATTCGCCAAAGTACAAGCGGGCCGCAAAAACTCGTTGGGCCGCTCACCGCTATTCCTGTGACCGAGATTTATACGGGGCAGGAAGAGGATAACACCGTGGAGCGGAAACGaaagatcggaagggagtaggggtgg                                >  2:62399/1‑228 (MQ=255)                                                                                                                                                                                                   gttttGGAAAATAACTAGCCTTTTTGGTGCATTATTGCTGTTGTTAATGCCGATAATGCTTATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAATATGCGATGCGCCAAAGGACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGTATAGGTTAAATCCGTGGCGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                               <  2:265459/247‑1 (MQ=255)                                                                                                                                                                                                         gAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGAGTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACGCGTTGGGTCGCTCATCGCTATTCCTGTGACCGAGCTTTAGACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAt                         <  1:258556/251‑1 (MQ=255)                                                                                                                                                                                                          aaaaTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAtt                        <  1:5429/251‑1 (MQ=255)                                                                                                                                                                                                                 cTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCgg                                                              <  1:262575/206‑1 (MQ=255)                                                                                                                                                                                                                   agcctgggtGGTGCAGTATTGCTGTTGTTAATGACGATAAGGCTGAGTAGGCAGGTGATTGTCGAACGTGCGGATTACCGTAGCGATGTGGAAGATGCGATTCGACAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGACGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGt                                                                     <  2:137533/189‑1 (MQ=255)                                                                                                                                                                                                                   aGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGt                                                                     >  1:137533/1‑197 (MQ=255)                                                                                                                                                                                                                   aGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttg                >  1:268620/1‑250 (MQ=255)                                                                                                                                                                                                                       ggTTTGGTGCAGTATTGCTGGTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGc           <  1:45423/250‑1 (MQ=255)                                                                                                                                                                                                                                cAGTACTGCTGTTGTTAAGTCCGATAATGCTGATTCGGAAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTAGACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAATTTTTAGCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAAc  <  1:404254/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           GGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAAC  >  NC_000913/4638001‑4638474

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