Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,329,835	G→A	pseudogene (4568/4605 nt)	yfaS ←	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,329,835	0	G	A	100.0%	204.7 / NA	61	pseudogene (4568/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/28);  total (33/28)

CCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAG  >  NC_000913/2329740‑2329931                                                                                                |                                                                                                 ccGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAac                                                                                              <  1:944983/100‑1 (MQ=255)  cGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCAAATCACACACTGCCacaaca                                                                                             >  2:814540/1‑100 (MQ=255)    cAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACt                                                                                           >  1:1309135/1‑100 (MQ=255)      aGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACttt                                                                                         <  1:972458/100‑1 (MQ=255)      aGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACttt                                                                                         <  1:1739453/100‑1 (MQ=255)        tgtTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTtgc                                                                                       >  2:361206/1‑100 (MQ=255)            ggCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTgcg                                                                                   >  1:275020/1‑100 (MQ=255)               agagTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAg                                                                                >  1:199447/1‑100 (MQ=255)                gagTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAgg                                                                               <  1:947381/100‑1 (MQ=255)                  gTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTg                                                                             <  2:624425/100‑1 (MQ=255)                   tAACCGACGTGAGTAGCCAGTCGGTTCTGCGCGAGTTCACTTATTTCATCTGCTTCCCGGTCCATTCACTCCCTGCCACAATACTTTGCTGCGCAGGTg                                                                             <  1:1667618/99‑1 (MQ=255)                    aaCCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCa                                                                           >  2:1214389/1‑100 (MQ=255)                    aaCCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCa                                                                           >  1:18768/1‑100 (MQ=255)                           gCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAAc                                                                    >  1:131337/1‑100 (MQ=255)                            ccAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACg                                                                   <  1:678145/100‑1 (MQ=255)                             cAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGc                                                                  >  1:1746461/1‑100 (MQ=255)                              aaTAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGca                                                                 <  1:757581/100‑1 (MQ=255)                                     aGACGATTCTGCGCCAGTTCTCTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACg                                                          <  2:1543878/100‑1 (MQ=255)                                     aGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACg                                                          >  2:287930/1‑100 (MQ=255)                                      gACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGa                                                         >  1:450145/1‑100 (MQ=255)                                      gACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTACATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGa                                                         <  1:1274877/100‑1 (MQ=255)                                        cGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGc                                                       >  1:1349049/1‑100 (MQ=255)                                           ttCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCagg                                                    <  1:1465571/100‑1 (MQ=255)                                             cTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCaggag                                                  <  1:288233/100‑1 (MQ=255)                                              tGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCaggagg                                                 >  1:868992/1‑100 (MQ=255)                                                cgcCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCa                                               >  2:1729429/1‑100 (MQ=255)                                                  ccAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGa                                             >  2:1507881/1‑100 (MQ=255)                                                    aGTTCACTTATTTCTCCTGCATTCTGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAAc                                           <  2:63340/100‑1 (MQ=255)                                                    aGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAAc                                           >  2:14596/1‑100 (MQ=255)                                                     gTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACg                                          >  2:633964/1‑100 (MQ=255)                                                          cTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTg                                     <  1:288243/100‑1 (MQ=255)                                                            tATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGcc                                   >  1:1418411/1‑100 (MQ=255)                                                            tATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGcc                                   >  2:1471128/1‑100 (MQ=255)                                                            tATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCCCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACAAATTGc                                    >  2:1715004/1‑99 (MQ=255)                                                             aTTTCACCTGCATCACGATCCATACACTCCCTGCCACAACACTTTGCTGCGCATGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGcc                                   >  1:261735/1‑99 (MQ=255)                                                              tttCACCTGCATCCCGCTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGcc                                   <  1:1293225/98‑1 (MQ=255)                                                               ttCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCtt                                >  2:708357/1‑100 (MQ=255)                                                                tCACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCttt                               >  2:668086/1‑100 (MQ=255)                                                                 cACCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCACtttt                              >  1:766079/1‑100 (MQ=255)                                                                  aCCTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCttttt                             <  2:1746474/100‑1 (MQ=255)                                                                   ccTGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTg                            >  1:10390/1‑100 (MQ=255)                                                                     tGCATCCGTCTCCATTCACTCCCTTCCACAACACTTTGCTGCGCAGGTGCATAGGAATGCACATAACGAGCAGGAGGTTGAACGAATTGCCCTTTTTGc                           <  1:1677427/99‑1 (MQ=255)                                                                     tGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCg                          <  2:131337/100‑1 (MQ=255)                                                                     tGCATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCg                          <  1:749508/100‑1 (MQ=255)                                                                       cATCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCgag                        <  2:1309145/100‑1 (MQ=255)                                                                        aTCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCgaga                       <  1:460954/100‑1 (MQ=255)                                                                         tctctgTCCAGTCTCGCCCTGCGACAACACTTGTCGGGGCAGTGGCATAGTAACGCACATAACGATCATGAGGCAGAACGGATTTCCCTTTTTGCgaga                       <  1:1727096/94‑1 (MQ=255)                                                                         tCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGaa                      >  2:1508870/1‑100 (MQ=255)                                                                         tCCCGGTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGaa                      >  1:1264155/1‑100 (MQ=255)                                                                           ccGGTCCATTCACTCCCTACCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAgc                    >  1:451981/1‑100 (MQ=255)                                                                            cgttcttttCACTCCCTGTCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAgc                    <  2:1632663/92‑1 (MQ=255)                                                                             ggTCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAgcgc                  <  1:1412862/100‑1 (MQ=255)                                                                               tCCATTCACTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCgcag                <  1:1212055/100‑1 (MQ=255)                                                                                      aCTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGAGCAACaaaa           >  2:761840/1‑97 (MQ=255)                                                                                       cTCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGc        <  2:622467/100‑1 (MQ=255)                                                                                        tCCCTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCg       <  2:275021/100‑1 (MQ=255)                                                                                         cccTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCgg      <  1:645842/100‑1 (MQ=255)                                                                                         cccTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCgg      >  1:1797473/1‑100 (MQ=255)                                                                                          ccTGCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAACAAATGGCGGa     >  2:1485173/1‑100 (MQ=255)                                                                                             gCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAg  >  1:359847/1‑100 (MQ=255)                                                                                             gCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAg  >  2:988947/1‑100 (MQ=255)                                                                                             gCCACAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAAGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAg  >  2:15908/1‑100 (MQ=255)                                                                                                |                                                                                                 CCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAG  >  NC_000913/2329740‑2329931

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