Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,175,797	(C)6→7	coding (210/582 nt)	yqiA ←	esterase YqiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,175,791	1	.	C	97.2%	121.8 / ‑3.8	36	coding (216/582 nt)	yqiA	esterase YqiA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (19/16);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑T‑CC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCA  >  NC_000913/3175663‑3175908                                                                                                                                   |                                                                                                                       tatTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑T‑cc‑cccccagtg                                                                                                            >  2:211008/1‑135 (MQ=255)  atTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACGGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGa                                                                                                            >  2:661872/1‑140 (MQ=255)     gctgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGa                                                                                                         <  2:590005/140‑1 (MQ=255)      ctgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGTACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGaa                                                                                                        >  1:408529/1‑140 (MQ=255)      ctgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTTCCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑T‑CCNCCCCAGTGACGaa                                                                                                        <  1:675003/140‑1 (MQ=255)      ctgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGaa                                                                                                        >  2:367337/1‑140 (MQ=255)      ctgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGaa                                                                                                        <  2:125923/140‑1 (MQ=255)       tgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAAc                                                                                                       <  2:462970/140‑1 (MQ=255)       tgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAAc                                                                                                       <  1:100631/140‑1 (MQ=255)                tAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGC‑GTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGa                                                                                                         <  1:125902/128‑1 (MQ=255)                tAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGC‑GTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGa                                                                                                         >  2:125902/1‑128 (MQ=255)                        ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAAc                                                                                                       >  2:466775/1‑123 (MQ=255)                        ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAAc                                                                                                       <  1:466775/123‑1 (MQ=255)                        ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGaa                                                                                      >  1:417116/1‑140 (MQ=255)                        ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGaa                                                                                      >  2:143454/1‑140 (MQ=255)                                gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGcc                                                                              >  1:177844/1‑140 (MQ=255)                                gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCAACGcc                                                                              >  1:586932/1‑140 (MQ=255)                                 gccGAGATGGTCCGCCAGCATTTCAAAGGGGCCCACCGCGGGTTTTCCCCCCACTGCGGGAAGCATAAAACATTGTGACACCCAGGTGGCGTAATA‑TCCC‑CCCCAGGGACGACCCGACAATACCCAGCGAATCACCGCCa                                                                             <  2:184803/139‑1 (MQ=255)                                 ccggggaTAGCCCGTCAGCAGTTAAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATACC‑CC‑CCCCAGTGACGACCCGACAATCCCCAGCGAATCACCGCCa                                                                             <  2:177844/135‑1 (MQ=255)                                 aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCGGCGAATCACCGCCa                                                                             <  2:394829/140‑1 (MQ=255)                                 aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCa                                                                             <  1:211008/140‑1 (MQ=255)                                             cGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGAcc                                                                 >  1:184803/1‑139 (MQ=255)                                              gTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACaa                                                                <  1:74119/140‑1 (MQ=255)                                                     gTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTc                                                         >  2:647151/1‑140 (MQ=255)                                                                      gccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGt                                        <  2:408529/140‑1 (MQ=255)                                                                                      cTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggcgg                        <  1:143454/140‑1 (MQ=255)                                                                                         cGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCaa                     >  2:271912/1‑140 (MQ=255)                                                                                           ggCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACt                   >  2:560649/1‑140 (MQ=255)                                                                                              agcaTAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggcgg                        >  1:275412/1‑132 (MQ=255)                                                                                              agcaTAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggcgg                        <  2:275412/132‑1 (MQ=255)                                                                                              agcaTAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                >  1:507374/1‑140 (MQ=255)                                                                                               gcaTAAAACATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCgg               >  2:156643/1‑140 (MQ=255)                                                                                                     aaCATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGat         <  2:374284/140‑1 (MQ=255)                                                                                                      aCATTGTGACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCAACGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatc        <  2:586932/140‑1 (MQ=255)                                                                                                          tgtgACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCaa                     <  2:289014/123‑1 (MQ=255)                                                                                                          tgtgACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCaa                     >  1:289014/1‑123 (MQ=255)                                                                                                          tgtgACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct    >  2:89215/1‑140 (MQ=255)                                                                                                            tgACAACCAGGTGGCGTAATA‑TCCC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCa  >  2:14473/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                       TATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATA‑T‑CC‑CCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCA  >  NC_000913/3175663‑3175908

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