Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	247,461	G→A	*250* (TAG→TAA)	yafL →	putative lipoprotein and C40 family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	247,461	0	G	A	100.0%	151.5 / NA	49	*250* (TAG→TAA)	yafL	putative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/28);  total (21/28)

TATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGC  >  NC_000913/247327‑247593                                                                                                                                       |                                                                                                                                     tATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACg                                                                                                                                  >  2:493594/1‑139 (MQ=255)         cGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACgat                                                                                                                                >  1:1169112/1‑131 (MQ=255)         cGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACgat                                                                                                                                <  2:1169112/133‑3 (MQ=255)           aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGtt                                                                                                                       >  1:930484/1‑139 (MQ=255)               gCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTAcg                                                                                                                   <  1:523847/139‑1 (MQ=255)                         gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCgg                                                                                                         >  1:1118756/1‑139 (MQ=255)                         gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCgg                                                                                                         <  1:1195380/139‑1 (MQ=255)                         gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCgg                                                                                                         <  1:275853/139‑1 (MQ=255)                         gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCgg                                                                                                         <  1:405727/139‑1 (MQ=255)                            tCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                      >  1:795322/1‑139 (MQ=255)                                 aCGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTa                                                                                                 <  1:900578/139‑1 (MQ=255)                                   gTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGTGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACg                                                                                               <  2:930484/139‑1 (MQ=255)                                              aCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTAcgc                                                                                                                  <  1:499854/109‑1 (MQ=255)                                              aCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTAcgc                                                                                                                  >  2:499854/1‑109 (MQ=255)                                              aCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGa                                                                                    >  1:959042/1‑139 (MQ=255)                                                cATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATg                                                                                  <  2:795322/139‑1 (MQ=255)                                                cATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATg                                                                                  <  2:374431/139‑1 (MQ=255)                                                          aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                        <  1:640423/139‑1 (MQ=255)                                                             cGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGcgcg                                                                     <  2:959042/139‑1 (MQ=255)                                                                     cGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCgg                                                                                                         <  1:367356/95‑1 (MQ=255)                                                                     cGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCgg                                                                                                         >  2:367356/1‑95 (MQ=255)                                                                      gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTAcg                                                                                                                   <  1:685009/84‑1 (MQ=255)                                                                      gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTAcg                                                                                                                   >  2:685009/1‑84 (MQ=255)                                                                             tCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACCCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACg                                                     >  2:929644/1‑139 (MQ=255)                                                                                    ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGt                                                                    <  1:885769/117‑1 (MQ=255)                                                                                    ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGt                                                                    >  2:885769/1‑117 (MQ=255)                                                                                      aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTgg                                                                         >  2:995856/1‑110 (MQ=255)                                                                                      aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTgg                                                                         <  1:995856/110‑1 (MQ=255)                                                                                       ccATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGaa                                           <  2:164465/139‑1 (MQ=255)                                                                                               tGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCg                                   >  2:609177/1‑139 (MQ=255)                                                                                                   cgcgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                            >  2:443619/1‑110 (MQ=255)                                                                                                   cgcgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                            <  1:443619/110‑1 (MQ=255)                                                                                                    gcgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCg                              >  2:858282/1‑139 (MQ=255)                                                                                                    gcgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCg                              >  1:860674/1‑139 (MQ=255)                                                                                                     cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGACgataa                                                                                                                              <  2:1242106/43‑5 (MQ=255)                                                                                                     cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGACgataa                                                                                                                              >  1:1242106/1‑39 (MQ=255)                                                                                                       cgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAg                                                                                            <  2:341015/74‑1 (MQ=37)                                                                                                       cgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAg                                                                                            >  1:341015/1‑74 (MQ=37)                                                                                                          aGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCCTCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                        <  1:858282/139‑1 (MQ=255)                                                                                                          aGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGTGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                        <  1:1071151/139‑1 (MQ=255)                                                                                                              ttttGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGat                    <  2:860674/139‑1 (MQ=255)                                                                                                                    cGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATaaaa              <  2:565450/139‑1 (MQ=255)                                                                                                                        agagACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGcgcg                                                                     <  1:357410/80‑1 (MQ=21)                                                                                                                        agagACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGcgcg                                                                     >  2:357410/1‑80 (MQ=21)                                                                                                                        agagACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                        >  1:1173159/1‑125 (MQ=255)                                                                                                                        agagACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGGTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                        <  2:1173159/125‑1 (MQ=255)                                                                                                                          agaCGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCt        <  1:493594/139‑1 (MQ=255)                                                                                                                            aCGATTTTGTAAAAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCa      >  2:39842/1‑139 (MQ=255)                                                                                                                                ttttGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGc  <  2:433679/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     TATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGC  >  NC_000913/247327‑247593

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