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%	155.9 / 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 (23/26);  total (23/26)

TGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGCC  >  NC_000913/247325‑247594                                                                                                                                         |                                                                                                                                      tgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGa                                                                                                                                     >  1:121971/1‑139 (MQ=255) tgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGa                                                                                                                                     >  1:2036095/1‑139 (MQ=255)  gtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAc                                                                                                                                    >  1:274505/1‑139 (MQ=255)     tttGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACgat                                                                                                                                 <  2:226099/139‑3 (MQ=255)     tttGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACgat                                                                                                                                 >  1:955082/1‑137 (MQ=255)             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCgat                                                                                                                        >  1:1208042/1‑137 (MQ=255)             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGtt                                                                                                                        >  2:1575536/1‑139 (MQ=255)             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGtt                                                                                                                        >  1:77400/1‑139 (MQ=255)             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGtt                                                                                                                        >  2:804022/1‑139 (MQ=255)                  cAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTAcgc                                                                                                                   <  2:716437/139‑1 (MQ=255)                      ttATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCa                                                                                                               <  2:2036095/139‑1 (MQ=255)                      ttATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCa                                                                                                               <  1:868457/139‑1 (MQ=255)                      ttATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCa                                                                                                               <  2:859485/139‑1 (MQ=255)                           gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCgg                                                                                                          <  1:804022/139‑1 (MQ=255)                             gTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCACATCCGGCa                                                                                                        <  1:1074943/139‑1 (MQ=255)                                   aCGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTCTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTa                                                                                                  <  1:1300968/139‑1 (MQ=255)                                        ccGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGA‑TTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                       >  2:86666/1‑128 (MQ=39)                                        ccGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGA‑TTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                       <  1:86666/128‑1 (MQ=39)                                          ggCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGt                                                                                           >  1:883309/1‑139 (MQ=255)                                           acGAAACCATACGGATAAGGCGATTAGCTGCACCTTTCTGGCATGACCCTTTCTTGGGCGCGCGCAGAAATATGACGGACGCGACGATCTTGTAAGA‑TGATATGGCGTTTACGCCGCATCCGGCAGTTGTACACAGGTg                                                                                          >  1:1976997/2‑139 (MQ=255)                                                  cATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATg                                                                                   <  1:1575536/139‑1 (MQ=255)                                                      cGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGAc                                                                               <  2:883309/139‑1 (MQ=255)                                                       ggATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACg                                                                              <  2:77400/139‑1 (MQ=255)                                                          tAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTg                                                                           >  2:588763/1‑139 (MQ=255)                                                           aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTgg                                                                          <  2:955082/139‑1 (MQ=255)                                                                        gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                             <  2:1854198/139‑1 (MQ=255)                                                                           ccTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                                                          >  1:864124/1‑139 (MQ=255)                                                                           ccTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                                                          >  2:1246030/1‑139 (MQ=255)                                                                              ttCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAg                                                                                             >  1:1318527/1‑101 (MQ=255)                                                                              ttCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAg                                                                                             <  2:1318527/101‑1 (MQ=255)                                                                                  ggCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGc                                                   <  2:274505/139‑1 (MQ=255)                                                                                        aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                       <  1:734991/81‑1 (MQ=38)                                                                                        aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                       >  2:734991/1‑81 (MQ=38)                                                                                        aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTgg                                                                          <  1:1428569/110‑1 (MQ=255)                                                                                        aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTgg                                                                          >  2:1428569/1‑110 (MQ=255)                                                                                        aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGa                                             >  2:1131982/1‑139 (MQ=255)                                                                                        aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGa                                             >  2:2000079/1‑139 (MQ=255)                                                                                                 tGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  2:1537412/102‑1 (MQ=255)                                                                                                 tGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         >  1:1537412/1‑102 (MQ=255)                                                                                                 tGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCg                                    <  1:1131982/139‑1 (MQ=255)                                                                                                            aGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                         <  2:864124/139‑1 (MQ=255)                                                                                                                        gAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGa                                                     >  2:36178/1‑99 (MQ=255)                                                                                                                        gAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGa                                                     <  1:36178/99‑1 (MQ=255)                                                                                                                        gAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCg                                    >  2:1296654/1‑116 (MQ=255)                                                                                                                        gAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCg                                    <  1:1296654/116‑1 (MQ=255)                                                                                                                        gAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAaaca             >  1:1248072/1‑139 (MQ=255)                                                                                                                          agagACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATc                                                                                                             <  2:338845/41‑1 (MQ=17)                                                                                                                          agagACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATc                                                                                                             >  1:338845/1‑41 (MQ=17)                                                                                                                            agaCGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCt         <  2:1248072/139‑1 (MQ=255)                                                                                                                                   tttGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGcc  <  1:1634262/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      TGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGCC  >  NC_000913/247325‑247594

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