Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,575,536	A→G	pseudogene (1402/2721 nt)	yjiV →	putative uncharacterized protein YjiV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,575,536	0	A	G	100.0%	173.1 / NA	48	pseudogene (1402/2721 nt)	yjiV	putative uncharacterized protein YjiV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/20);  total (28/20)

TCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCA  >  NC_000913/4575406‑4575671                                                                                                                                   |                                                                                                                                        tCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGa                                                                                                                                >  2:382840/1‑140 (MQ=255)          tAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCAt                                                                                                                       <  1:382840/140‑1 (MQ=255)          tAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCAt                                                                                                                       <  1:659384/140‑1 (MQ=255)           aaTCAGAACAGACGCTGGCCTTAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGGAAGAAGCATg                                                                                                                      >  1:806653/1‑140 (MQ=255)           aaTCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATg                                                                                                                      >  1:264059/1‑140 (MQ=255)             tCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATggg                                                                                                                    <  2:264059/140‑1 (MQ=255)                  aCAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCaa                                                                                                               >  2:197552/1‑140 (MQ=255)                      aCGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGcc                                                                                                                  <  2:713638/133‑1 (MQ=255)                      aCGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGcc                                                                                                                  >  1:713638/1‑133 (MQ=255)                      aCGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTa                                                                                                           >  2:169360/1‑140 (MQ=255)                      aCGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTa                                                                                                           >  2:580073/1‑140 (MQ=255)                      aCGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTa                                                                                                           >  2:606108/1‑140 (MQ=255)                      aCGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTa                                                                                                           >  2:609771/1‑140 (MQ=255)                           ggCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGa                                                                                                      >  1:560791/1‑140 (MQ=255)                             ccTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGAcc                                                                                                    <  1:580073/140‑1 (MQ=255)                               tCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGc                                                                                                  <  1:527470/140‑1 (MQ=255)                                cAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCt                                                                                                 <  2:662672/140‑1 (MQ=255)                                                       ccATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCa                                                                          >  2:44678/1‑140 (MQ=255)                                                       ccATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCa                                                                          >  2:426799/1‑140 (MQ=255)                                                           cGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCg                                                                      >  1:23284/1‑140 (MQ=255)                                                            gTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAAc                                                                        >  1:556892/1‑137 (MQ=255)                                                            gTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAAc                                                                        <  2:556892/137‑1 (MQ=255)                                                             tttAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGtt                                                                    >  2:310852/1‑140 (MQ=255)                                                               tAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGtttt                                                                  >  1:731979/1‑140 (MQ=255)                                                                 gCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTttta                                                                <  1:477545/140‑1 (MQ=255)                                                                      gTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTAc                                                           <  1:126092/140‑1 (MQ=255)                                                                       ttGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTAcc                                                          >  1:100827/1‑140 (MQ=255)                                                                               gCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTtttattt                                                             >  2:736456/1‑129 (MQ=255)                                                                               gCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTtttattt                                                             <  1:736456/129‑1 (MQ=255)                                                                                 tGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCat                                                <  1:426799/140‑1 (MQ=255)                                                                                                    cGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGGTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCa                             <  2:275819/140‑1 (MQ=255)                                                                                                      tCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACg                           >  1:791602/1‑140 (MQ=255)                                                                                                      tCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACg                           <  2:731979/140‑1 (MQ=255)                                                                                                      tCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACg                           <  1:606108/140‑1 (MQ=255)                                                                                                      tCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACg                           <  1:609771/140‑1 (MQ=255)                                                                                                      tCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACg                           <  2:23284/140‑1 (MQ=255)                                                                                                                  cTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt               <  1:15915/140‑1 (MQ=255)                                                                                                                   tATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAttt              <  2:806653/140‑1 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGc                      >  1:96491/1‑122 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGc                      <  2:96491/122‑1 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAc    >  1:657477/1‑140 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAc    >  1:548159/1‑140 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAc    >  1:671021/1‑140 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAc    >  1:230962/1‑140 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAc    >  2:740369/1‑140 (MQ=255)                                                                                                                             atCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAc    >  2:853561/1‑140 (MQ=255)                                                                                                                               cAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCa  >  2:655393/1‑140 (MQ=255)                                                                                                                               cAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCa  >  1:675965/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                        TCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCA  >  NC_000913/4575406‑4575671

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