Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,334,336	C→T	pseudogene (67/4605 nt) *208* (TAG→TAA)	yfaS ← yfaT ←	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein DUF1175 family protein, putative host defense protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,334,336	0	C	T	100.0%	190.9 / NA	60	pseudogene (67/4605 nt) *208* (TAG→TAA)	yfaS yfaT	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein DUF1175 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/31);  total (29/31)

CTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGG  >  NC_000913/2334209‑2334469                                                                                                                                |                                                                                                                                      ctgctAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATcgcg                                                                                                                                <  1:325923/135‑1 (MQ=255) ctgctAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATcgcg                                                                                                                                >  2:325923/1‑135 (MQ=255)      aaaaCTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAAt                                                                                                                       >  2:1335188/1‑139 (MQ=255)               gTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGAt                                                                                                              <  1:35028/139‑1 (MQ=255)                   gcaagcaaAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAAt                                                                                                          <  2:79181/139‑1 (MQ=255)                        caaAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                     <  1:856714/139‑1 (MQ=255)                        caaAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                     <  2:873487/139‑1 (MQ=255)                                  gTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTggg                                                                                           >  2:680316/1‑139 (MQ=255)                                      ccccGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTg                                                                                       >  1:241190/1‑139 (MQ=255)                                             aGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTGTCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAtt                                                                                              <  2:515082/125‑1 (MQ=255)                                             aGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTGTCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAtt                                                                                              >  1:515082/1‑125 (MQ=255)                                               gcggcgCATAGTTGCTGGAAGGAAGCGAATCATTAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTg                                                                              <  2:1452055/139‑1 (MQ=255)                                               gcggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTg                                                                              >  1:1452055/1‑139 (MQ=255)                                                  gcgcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTa                                                                                                                    >  2:266034/1‑98 (MQ=255)                                                  gcgcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTa                                                                                                                    <  1:266034/98‑1 (MQ=255)                                                    gcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTAt                                                                         <  2:694001/139‑1 (MQ=255)                                                             cTGGAAGGCAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGCATCCATCGGgt                                                                >  2:1048680/1‑139 (MQ=255)                                                             cTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                                >  2:383118/1‑139 (MQ=255)                                                              tGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgtg                                                               <  1:1335188/139‑1 (MQ=255)                                                                  aggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCATGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATTAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                                <  2:125987/134‑1 (MQ=255)                                                                  aggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                                >  1:125987/1‑134 (MQ=255)                                                                  aggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCt                                                           >  1:245903/1‑139 (MQ=255)                                                                  aggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCt                                                           >  2:795384/1‑139 (MQ=255)                                                                    gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTc                                                         <  2:241190/139‑1 (MQ=255)                                                                     aagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGa                                                                                          <  2:38615/105‑1 (MQ=255)                                                                     aagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGa                                                                                          >  1:38615/1‑105 (MQ=255)                                                                     aagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTcc                                                        <  1:383118/139‑1 (MQ=255)                                                                       gCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  <  2:191307/79‑1 (MQ=255)                                                                       gCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  >  1:191307/1‑79 (MQ=255)                                                                               tcaGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                              <  2:1126203/139‑1 (MQ=255)                                                                                     ttAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGttgttg                                        >  2:194630/1‑139 (MQ=255)                                                                                          aaGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAgac                                   <  1:680316/139‑1 (MQ=255)                                                                                           aGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAgact                                  <  1:795384/139‑1 (MQ=255)                                                                                                tGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGactgactg                             >  2:1536445/1‑139 (MQ=255)                                                                                                 gTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGaa                                                                                    <  1:477915/83‑1 (MQ=255)                                                                                                 gTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGaa                                                                                    >  2:477915/1‑83 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTg                                                                                       <  1:1600279/79‑1 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTg                                                                                       >  2:1600279/1‑79 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                            <  2:76693/90‑1 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                            >  1:76693/1‑90 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATcc                                                                       >  1:775950/1‑95 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATcc                                                                       <  2:775950/95‑1 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTgcg                           <  2:245903/139‑1 (MQ=255)                                                                                                  ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTgcg                           >  2:842209/1‑139 (MQ=255)                                                                                                       acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGttgttg                                        <  1:94140/121‑1 (MQ=255)                                                                                                       acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGttgttg                                        >  2:94140/1‑121 (MQ=255)                                                                                                        caAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCAtt                     >  2:138072/1‑139 (MQ=255)                                                                                                           gagaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATTAAATTGGGATTGGAGTCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCg                  <  1:194630/139‑1 (MQ=255)                                                                                                            agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGt                 >  1:1515482/1‑139 (MQ=255)                                                                                                            agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGt                 >  1:92235/1‑139 (MQ=255)                                                                                                              aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAAt                                                                                   >  2:970367/1‑71 (MQ=255)                                                                                                              aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAAt                                                                                   <  1:970367/71‑1 (MQ=255)                                                                                                                        gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAAt                                                                                                    <  2:130665/44‑1 (MQ=255)                                                                                                                        gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAAt                                                                                                    >  1:130665/1‑44 (MQ=255)                                                                                                                        gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTc                                                             >  2:404465/1‑83 (MQ=255)                                                                                                                        gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTc                                                             <  1:404465/83‑1 (MQ=255)                                                                                                                        gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                  >  2:164681/1‑94 (MQ=255)                                                                                                                        gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                  <  1:164681/94‑1 (MQ=255)                                                                                                                           cGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTgg  <  2:92235/139‑1 (MQ=255)                                                                                                                           cGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTgg  <  1:1297363/139‑1 (MQ=255)                                                                                                                             ctgttATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAt                                                                                               >  1:898053/1‑44 (MQ=255)                                                                                                                                |                                                                                                                                      CTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGG  >  NC_000913/2334209‑2334469

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