Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,071,172	C→A	pseudogene (387/426 nt)	yeeP →	CP4‑44 prophage; putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,071,172	0	C	A	100.0%	159.0 / NA	46	pseudogene (387/426 nt)	yeeP	CP4‑44 prophage; putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/25);  total (21/25)

GAATCTGTCCGGTCTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTCTCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAG  >  NC_000913/2071044‑2071297                                                                                                                                 |                                                                                                                              gAATCTGTCCGGTCTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGCATCTGTTGCCCg                                                                                                                    >  2:129573/1‑140 (MQ=255)   aTCTGTCCGGTTTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGGCCGGATATTTGACGCGGCGGAGGGCGTCTGTATTGCCTTTGTTGCGCGCACGGTTCTGCGTGCCGTCTGTGACACGGTGGTATCTGTTGCCCGTg                                                                                                                  <  2:285282/140‑1 (MQ=255)   aTCTGTCCGGTCTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTg                                                                                                                  <  2:140241/140‑1 (MQ=255)   aTCTGTCCGGTCTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTg                                                                                                                  <  2:86633/140‑1 (MQ=255)     cTGTCCGGTCTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCg                                                                                                                <  2:331518/140‑1 (MQ=255)            gTCTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGa                                                                                                         >  1:251985/1‑140 (MQ=255)                    cccGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGAtt                                                                                                 <  2:318775/140‑1 (MQ=255)                               tttACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAAc                                                                                      >  2:93093/1‑140 (MQ=255)                               tttACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAAc                                                                                      <  1:474429/140‑1 (MQ=255)                               tttACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAAc                                                                                      <  1:284355/140‑1 (MQ=255)                                      gtgcggtgGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGt                                                                               <  2:23267/140‑1 (MQ=255)                                       tgcggtgGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTg                                                                              >  2:89704/1‑140 (MQ=255)                                         cggtgGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGAc                                                                            >  2:615877/1‑140 (MQ=255)                                         cggtgGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGAc                                                                            >  2:226984/1‑140 (MQ=255)                                         cggtgGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGGATGGAACTGGATTTTCTTCTGAACCTGTCGTGAc                                                                            >  1:123157/1‑140 (MQ=255)                                             ggACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGAt                                                                        <  1:318627/140‑1 (MQ=255)                                                   gATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGccctc                                                                  >  1:560493/1‑140 (MQ=255)                                                   gATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGccctc                                                                  >  1:252046/1‑140 (MQ=255)                                                     tatTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGccctccc                                                                <  1:679567/140‑1 (MQ=255)                                                      atTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGccctccct                                                               <  2:567733/140‑1 (MQ=255)                                                      atTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGccctccct                                                               <  2:87689/140‑1 (MQ=255)                                                                gcggAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGt                                                     <  2:251985/140‑1 (MQ=255)                                                                     gagCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCt                                                <  1:447539/140‑1 (MQ=255)                                                                         gTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTcaca                                            >  2:514384/1‑140 (MQ=255)                                                                          tCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACaa                                           >  2:224266/1‑140 (MQ=255)                                                                           cTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACaaa                                          >  1:350246/1‑140 (MQ=255)                                                                             gTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAg                                        <  2:560493/140‑1 (MQ=255)                                                                               aTTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCa                                      <  2:421494/140‑1 (MQ=255)                                                                                        gTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCg                             >  1:613308/1‑140 (MQ=255)                                                                                        gTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCg                             <  1:226984/140‑1 (MQ=255)                                                                                        gTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCg                             >  2:624526/1‑140 (MQ=255)                                                                                          tGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGtt                           <  1:423848/140‑1 (MQ=255)                                                                                              cgcacgGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTg                       <  1:583261/140‑1 (MQ=255)                                                                                              cgcacgGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTg                       <  2:548040/140‑1 (MQ=255)                                                                                               gcacgGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGt                      <  1:430751/140‑1 (MQ=255)                                                                                                 acgGTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGtct                    >  2:180295/1‑140 (MQ=255)                                                                                                   ggTTCTGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGtctct                  <  1:224266/140‑1 (MQ=255)                                                                                                        tGCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGt             <  2:435128/140‑1 (MQ=255)                                                                                                         gCGTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTc            >  2:80485/1‑140 (MQ=255)                                                                                                           gTGCCGTCCGTGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCg          <  1:514384/140‑1 (MQ=255)                                                                                                               cgtccgtGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCa      >  1:93343/1‑140 (MQ=255)                                                                                                               cgtccgtGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCa      >  2:631203/1‑140 (MQ=255)                                                                                                               cgtccgtGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCa      >  2:320988/1‑140 (MQ=255)                                                                                                                 tccgtGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAAt    >  2:656809/1‑140 (MQ=255)                                                                                                                 tccgtGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAAt    >  2:436369/1‑140 (MQ=255)                                                                                                                 tccgtGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAAt    <  1:129573/140‑1 (MQ=255)                                                                                                                   cgtGACACGGTGGTATCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAg  <  1:615877/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                              GAATCTGTCCGGTCTCAGGCCCGTGAACAGTTTACCGGTGCGGTGGACCGGATATTTGACACGGCGGAGAGCGTCTGTATTGCCTCTGTTGCACGCACGGTTCTGCGTGCCGTCCGTGACACGGTGGTCTCTGTTGCCCGTGCGGTATGGAACTGGATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAG  >  NC_000913/2071044‑2071297

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