Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	890,350	T→C	22.7%	S88P (TCA→CCA)	ybjM →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	890,350	0	T	C	22.7%	43.8 / 7.3	29	S88P (TCA→CCA)	ybjM	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (10/12);  new base C (4/3);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCACAGAAAAAATCAGTAAAGCCCTCAACGCGA  >  NC_000913/890209‑890480                                                                                                                                              |                                                                                                                                   gtgttGCTATTTTTCAGTCTTCCTGGATCAGTCTCCAGCTTCTTGTCACAGCTTAGATAAGTCCTGAAACCTCTGTTGGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCtt                                                                                                                             <  2:103796/146‑2 (MQ=255) ggCTTGCTATTTTTCATTCTTCCAGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCTGCCCCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGCCGCGCTCCTTCTg                                                                                                                             <  2:492198/149‑1 (MQ=255)  gCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCCCATTCTgg                                                                                                                            >  2:700935/1‑149 (MQ=255)  gCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTAGGCGGCTGTTTTTTGCACCGACGGGCCCATTCTgg                                                                                                                            >  2:743541/1‑149 (MQ=255)  gCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCAGTTTGGCCCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCCCAATcgtg                                                                                                                            >  2:277488/1‑146 (MQ=255)        tATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAagag                                                                                                                      >  1:706809/1‑149 (MQ=255)           ttttCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACCGCGTAGAGAAGTCCTGAAACCTCTGTGTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTg                                                                                                                   <  2:9008/149‑1 (MQ=255)           tttgCATTCTGCCTGGGGCAGTCGCCAGCTTCGTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTg                                                                                                                   <  2:9003/149‑1 (MQ=255)                           aGCAGTCGCCAGCTTTTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAgcg                                                                                                   <  2:731665/149‑1 (MQ=255)                           aGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAgcg                                                                                                   <  1:127451/149‑1 (MQ=255)                                            tttCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCCCATTCTGGCAAGAGCTGGCCTGGTTACTAAGCGCGGTGTTCCGGTGcgcg                                                                                  >  1:212836/1‑149 (MQ=255)                                                   gCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTggggg                                                                           >  1:108457/1‑149 (MQ=255)                                                   gCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAAGGCGGGCGGCACCCGGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGGCGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTggggg                                                                           <  1:387153/149‑1 (MQ=255)                                                               ccggAAACCCCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCt                                                               <  2:661055/146‑1 (MQ=255)                                                                  gAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCttat                                                            >  2:598088/1‑149 (MQ=255)                                                                  gAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGAGGTGTTCTGATGTGCGCTGGGGGCACTGTGTTTCttat                                                            >  1:70886/1‑149 (MQ=255)                                                                        tctGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCa                                                      >  1:421511/1‑149 (MQ=255)                                                                        tctGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCCCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCa                                                      <  1:700935/149‑1 (MQ=255)                                                                        tctGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTGTTTTTCACCGACGCGCCCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCa                                                      <  1:277488/149‑1 (MQ=255)                                                                        ccggTTTGGCGCAATGCGGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCCCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCa                                                      <  1:743541/146‑1 (MQ=255)                                                                            tttGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACGGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCagtag                                                  >  1:416636/1‑149 (MQ=255)                                                                                   cAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAgtttgttt                                           <  2:663029/149‑1 (MQ=255)                                                                                        cTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTTGCATGGTTACTAAGAGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAAcc                                      >  2:729832/1‑149 (MQ=255)                                                                                        cTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAAcc                                      >  2:251042/1‑149 (MQ=255)                                                                                                   cTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCACAGaa                           >  1:328795/1‑149 (MQ=255)                                                                                                   cTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCACAGaa                           >  1:244582/1‑149 (MQ=255)                                                                                                           tGCTCCTTATGCGGCTGTTTTTTTCACCGACGCGCTCAGTCTGGCAAGGGCGGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCACAGAAAAAATCAg                   <  2:15836/149‑1 (MQ=255)                                                                                                                          tGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCACAGAAAAAATCAGTAAAGCCCTCAAcgc    <  1:598088/149‑1 (MQ=255)                                                                                                                            tttttttCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCACAGAAAAAATCAGTAAAGCCCTCAACGCGa  <  1:214180/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                   GGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCACAGAAAAAATCAGTAAAGCCCTCAACGCGA  >  NC_000913/890209‑890480

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