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_000913247,4610GA100.0% 136.5 / NA 45*250* (TAG→TAA) yafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/23);  total (22/23)

TTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTC  >  NC_000913/247330‑247588
                                                                                                                                   |                                                                                                                               
ttGGGCGATGGGGAATTTTTCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGTCGGAAGAGACGATTTTGTAAGACgata                                                                                                                          <  2:916070/139‑4 (MQ=255)
               tttATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTAcgccgc                                                                                                          >  1:1350601/1‑139 (MQ=255)
                ttATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCa                                                                                                         <  2:688372/139‑1 (MQ=255)
                  aTCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATc                                                                                                       >  2:1202239/1‑139 (MQ=255)
                      aGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGCCGGAGGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCTGCATCCGGc                                                                                                   <  2:1016563/139‑1 (MQ=255)
                          gCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGtt                                                                                               >  1:173643/1‑139 (MQ=255)
                                     gCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTg                                                                                    >  1:1505928/1‑139 (MQ=255)
                                                cGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGAc                                                                         >  1:302538/1‑139 (MQ=255)
                                                        ccGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACg                                                                                          <  1:1065836/114‑1 (MQ=255)
                                                        ccGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACg                                                                                          >  2:1065836/1‑114 (MQ=255)
                                                         cGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGcgcg                                                                <  2:302538/139‑1 (MQ=255)
                                                             tAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCtt                                                            >  1:488332/1‑139 (MQ=255)
                                                                 cGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCa                                                        <  1:1202239/139‑1 (MQ=255)
                                                                  gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACg                                                                        <  2:1533894/122‑1 (MQ=255)
                                                                  gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACg                                                                        >  1:1533894/1‑122 (MQ=255)
                                                                         tCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACg                                                >  2:690793/1‑139 (MQ=255)
                                                                              cAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGccc                                           <  1:958723/139‑1 (MQ=255)
                                                                               aGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCcg                                          <  2:39137/139‑1 (MQ=255)
                                                                                ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATg                                                                             <  2:916659/103‑1 (MQ=255)
                                                                                ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATg                                                                             >  1:916659/1‑103 (MQ=255)
                                                                                ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCcgc                                         >  2:1103661/1‑139 (MQ=255)
                                                                                 gACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACgat                                                                                                                           <  2:785684/57‑3 (MQ=255)
                                                                                 gACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACgat                                                                                                                           >  1:785684/1‑55 (MQ=255)
                                                                                  aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTa                                                                                            <  1:746631/86‑1 (MQ=255)
                                                                                  aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                   <  2:100621/111‑1 (MQ=255)
                                                                                  aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                   >  1:100621/1‑111 (MQ=255)
                                                                                  aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGt                                                               >  2:474723/1‑115 (MQ=255)
                                                                                  aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGt                                                               <  1:474723/115‑1 (MQ=255)
                                                                                  aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATACGGCAGTTGTa                                                                                            >  2:746631/1‑86 (MQ=38)
                                                                                             ggCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAg                            >  1:1147776/1‑139 (MQ=255)
                                                                                             ggCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAg                            <  2:1505928/139‑1 (MQ=255)
                                                                                                 cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCa                                                                                                         >  2:1347676/1‑58 (MQ=37)
                                                                                                 cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCa                                                                                                         <  1:1347676/58‑1 (MQ=37)
                                                                                                 cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                 >  1:1488183/1‑66 (MQ=25)
                                                                                                 cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                 <  2:1488183/66‑1 (MQ=25)
                                                                                                 cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGc                        >  1:867488/1‑139 (MQ=255)
                                                                                                  gcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACg                                                                                          <  1:1084541/72‑1 (MQ=37)
                                                                                                  gcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGAGTTTACGCCGCATCCGGCAGTTGTACg                                                                                          >  2:1084541/1‑72 (MQ=37)
                                                                                                      aGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTTCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                   <  1:836221/139‑1 (MQ=37)
                                                                                                      aGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                   <  2:173643/139‑1 (MQ=255)
                                                                                                      aGGATTTTGACGGAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                   <  2:867488/139‑1 (MQ=255)
                                                                                                                  gAAGGGTCGATTTTGTAAGAC‑GATAAGGCGTTGAGGCCGCAACCGGCAGTTGTACGCAGCTGCCCAATGCGCCGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGGGAGCGGCTTTCCCGATATAAgaca       <  2:94582/139‑1 (MQ=21)
                                                                                                                  gAAGAGACGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAaaca       >  1:852207/1‑139 (MQ=255)
                                                                                                                       gaCGATTTTGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTc  <  2:488332/139‑1 (MQ=255)
                                                                                                                             tttGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGcgcg                                                                <  2:1288092/71‑1 (MQ=17)
                                                                                                                             tttGTAAGAC‑GATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGcgcg                                                                >  1:1288092/1‑71 (MQ=17)
                                                                                                                                   |                                                                                                                               
TTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTC  >  NC_000913/247330‑247588

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