Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,484,088	(C)6→7	coding (887/1023 nt)	yheT →	putative hydrolase YheT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,484,082	1	.	C	86.8%	114.0 / 10.6	38	coding (881/1023 nt)	yheT	putative hydrolase YheT
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (16/17);  total (19/19)
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

CGCCGACGCTATCGACTATTATCGTCAGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑C‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTGACTGGTTAACAACGTATCTGGAGG  >  NC_000913/3483948‑3484210                                                                                                                                         |                                                                                                                                    cGCCGACGCTATCGACTATTATCGTCAGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑c‑cc‑cc                                                                                                                              >  1:466518/1‑140 (MQ=255)         cTATCGACTATTATCGTCAGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGt                                                                                                                       <  1:65341/140‑1 (MQ=255)           aTCGACTATTATCGTCAGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTgg                                                                                                                     >  2:485893/1‑140 (MQ=255)                    tatCGTCAGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAac                                                                                                            <  1:18147/140‑1 (MQ=255)                       cGTCAGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAactga                                                                                                         <  1:660543/140‑1 (MQ=255)                           aGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAActgactga                                                                                                     >  1:417536/1‑140 (MQ=255)                             tgtAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCTGAAACCGGAAAGGC‑T‑CACC‑CCCGCAAGTGGGGTATCAACTGACTGAAc                                                                                                   >  2:205047/1‑140 (MQ=255)                             tgtAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGAATCAACTGACTGAAc                                                                                                   >  2:117305/1‑140 (MQ=255)                               tAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGTAACCGGAAGGTCGC‑C‑CC‑CCCGCAGGGGGAGTATCAACTGACTGAACAt                                                                                                 >  2:263116/1‑140 (MQ=255)                                 gcgcCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCAC‑C‑CC‑CCCACAGATGAAGTATCGACTGACTGAACATgg                                                                                               >  1:173568/1‑140 (MQ=255)                                   gcCATGCCGATGCTGAACCGGAGCGACAACCCGACGCCGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCg                                                                                             <  1:205047/140‑1 (MQ=255)                                    ccATGCCGATGCTGAACCGGATTGCCAAACCTACGTTGACTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCgg                                                                                            <  2:173568/140‑1 (MQ=255)                                       tGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCa                                                                                         >  2:69465/1‑140 (MQ=255)                                       tGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCa                                                                                         >  1:649255/1‑140 (MQ=255)                                       tGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCa                                                                                         >  1:612684/1‑140 (MQ=255)                                       tGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCa                                                                                         >  1:511697/1‑140 (MQ=255)                                       tGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCa                                                                                         >  1:425009/1‑140 (MQ=255)                                       tGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCa                                                                                         <  1:26415/140‑1 (MQ=255)                                            aTGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTg                                                                                    >  2:644416/1‑140 (MQ=255)                                                 gAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCttt                                                                               <  2:511697/140‑1 (MQ=255)                                                   aCCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTAt                                                                             >  1:434311/1‑140 (MQ=255)                                                      ggATCGCCAAACCGACGATGATCATTCACGCCAAACACGCTCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTgg                                                                          <  2:174343/140‑1 (MQ=255)                                                       gATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGc                                                                         >  2:109215/1‑140 (MQ=255)                                                        aTCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑g                                                                       <  1:109215/140‑1 (MQ=255)                                                        aTCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑g                                                                       <  2:612684/140‑1 (MQ=255)                                                                     cGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCa                                                          <  2:649255/140‑1 (MQ=255)                                                                                    ccAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCGTTATTGGA‑GGTACAATACTTCTTCCCCAAATGTGGCt                                           >  1:301661/1‑140 (MQ=255)                                                                                      aaaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGGACATTACTTCATCCGCCAAAgtgggtgg                                         >  2:373565/1‑140 (MQ=255)                                                                                          aCGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTc                                     <  2:109845/140‑1 (MQ=255)                                                                                            gATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCAc                                   <  2:466518/140‑1 (MQ=255)                                                                                                    tATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCt                           <  2:425009/140‑1 (MQ=255)                                                                                                     aTGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTg                          <  2:417536/140‑1 (MQ=255)                                                                                                      tGGATCATCAGGTGATTCCGAACCCGGAAAGCC‑C‑CTCC‑CCCGCAGGTGCAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTGa                         <  2:301661/140‑1 (MQ=255)                                                                                                            atcaGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTGACTGGtt                   >  1:99863/1‑140 (MQ=255)                                                                                                            atcaGGTGATCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTGACTGGtt                   <  1:644416/140‑1 (MQ=255)                                                                                                              caGGTGATCCCGAAACCGGAAATTC‑A‑C‑CCACCCGCAGTTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTGACTGGTTaa                 <  1:320740/140‑1 (MQ=255)                                                                                                                     tCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGGTAATTTATTTTTTCCCCCAA‑GTGGGGGGGGCTACACATTCCTGACTGGTTaacaaagga          >  2:479900/1‑136 (MQ=255)                                                                                                                     tCCCGAAACCGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGGACATTTCTTCATCCCCCAAAGGGGCTGGGGGCAGCCATTTCTGGCTGGGTTAAAACgga          >  1:305482/1‑138 (MQ=255)                                                                                                                             ccGGAAAGTC‑TCC‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTGACTGGTTAACAACGTATCTGGAgg  <  2:434311/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    CGCCGACGCTATCGACTATTATCGTCAGTGTAGCGCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑C‑CC‑CCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGC‑GGTACATTACTTCATCCGCAAATGTGGCTGGAGTCACGCATTCCTGACTGGTTAACAACGTATCTGGAGG  >  NC_000913/3483948‑3484210

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