Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,544,384	C→T	pseudogene (381/381 nt)	yddK ←	pseudogene, leucine‑rich protein, putative glycoportein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,544,384	0	C	T	100.0%	231.5 / NA	71	pseudogene (381/381 nt)	yddK	pseudogene, leucine‑rich protein, putative glycoportein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/34);  total (37/34)

AACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTTTTAT  >  NC_000913/1544261‑1544518                                                                                                                            |                                                                                                                                       aaCGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGaa                                                                                                                         >  2:437202/1‑139 (MQ=255)              tttGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACg                                                                                                            <  1:437202/139‑1 (MQ=255)                  aaCGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCt                                                                                                        >  2:444360/1‑139 (MQ=255)                                ttCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTa                                                                                          >  2:937828/1‑139 (MQ=255)                                         tCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 >  1:1164804/1‑123 (MQ=255)                                         tCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 <  2:1164804/123‑1 (MQ=255)                                               cGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     <  2:1281677/129‑1 (MQ=255)                                               cGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     >  1:1281677/1‑129 (MQ=255)                                                   aaTCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                                       <  2:1146544/139‑1 (MQ=255)                                                                aCATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCa                                                                                                                             <  1:942636/72‑1 (MQ=255)                                                                aCATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCa                                                                                                                             >  2:942636/1‑72 (MQ=255)                                                                aCATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCt                                                                     >  1:1283750/1‑128 (MQ=255)                                                                aCATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCt                                                                     <  2:1283750/128‑1 (MQ=255)                                                                      tAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     >  1:1397442/1‑106 (MQ=255)                                                                      tAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     <  2:1397442/106‑1 (MQ=255)                                                                      tAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGtta                                                    <  2:693476/139‑1 (MQ=255)                                                                           atatTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACa                                               >  2:1367709/1‑139 (MQ=255)                                                                            tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTaaaa                                                                                                                                   >  2:258802/1‑54 (MQ=255)                                                                            tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTaaaa                                                                                                                                   <  1:258802/54‑1 (MQ=255)                                                                             attatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                                          >  1:843769/1‑78 (MQ=255)                                                                             attatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                                          <  2:843769/78‑1 (MQ=255)                                                                               tatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCgg                                                                                      <  2:679690/96‑1 (MQ=255)                                                                               tatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCgg                                                                                      >  1:679690/1‑96 (MQ=255)                                                                                atCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 <  1:1403306/84‑1 (MQ=255)                                                                                atCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 >  2:1403306/1‑84 (MQ=255)                                                                                      aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                                        >  1:244826/1‑135 (MQ=255)                                                                                      aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                                        <  2:244826/135‑1 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                 >  2:731068/1‑92 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                 <  2:742897/92‑1 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                 <  1:731068/92‑1 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                 >  1:742897/1‑92 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         >  1:1045869/1‑116 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         <  2:1045869/116‑1 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttatt                                                  >  2:224637/1‑123 (MQ=255)                                                                                        gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttatt                                                  <  1:224637/123‑1 (MQ=255)                                                                                          ccAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAt                                <  1:1367709/139‑1 (MQ=255)                                                                                           cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                 >  1:140265/1‑89 (MQ=255)                                                                                           cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                 <  2:140265/89‑1 (MQ=255)                                                                                           cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                                   <  2:1228477/119‑1 (MQ=255)                                                                                           cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                                   >  1:1228477/1‑119 (MQ=255)                                                                                             gACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 <  2:687929/71‑1 (MQ=255)                                                                                             gACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 >  1:687929/1‑71 (MQ=255)                                                                                                   tataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAt                                <  2:1310346/130‑1 (MQ=255)                                                                                                   tataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAt                                >  1:1310346/1‑130 (MQ=255)                                                                                                    ataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     <  2:1411069/76‑1 (MQ=255)                                                                                                    ataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     >  1:1411069/1‑76 (MQ=255)                                                                                                       gAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                           >  1:369083/1‑99 (MQ=255)                                                                                                       gAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                           <  2:369083/99‑1 (MQ=255)                                                                                                          ttAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 >  2:509794/1‑58 (MQ=255)                                                                                                          ttAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                 <  1:509794/58‑1 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                                          <  1:435052/47‑1 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                                          >  2:435052/1‑47 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                           <  2:1244671/94‑1 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                           >  1:1244671/1‑94 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         <  1:1426371/96‑1 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         <  1:1156638/96‑1 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         >  2:1426371/1‑96 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         >  2:1156638/1‑96 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                                        >  1:996452/1‑113 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                                        <  2:996452/113‑1 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtt              >  2:820269/1‑139 (MQ=255)                                                                                                            aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtt              >  1:1430793/1‑139 (MQ=255)                                                                                                              tgatgaTTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                                       <  2:463985/80‑1 (MQ=255)                                                                                                              tgatgaTTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                                       >  1:463985/1‑80 (MQ=255)                                                                                                              tgatgaTTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                             >  2:684728/1‑106 (MQ=255)                                                                                                              tgatgaTTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                             <  1:684728/106‑1 (MQ=255)                                                                                                               gatgatTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     <  2:685790/65‑1 (MQ=255)                                                                                                               gatgatTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                     >  1:685790/1‑65 (MQ=255)                                                                                                                  gatTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         <  2:1060580/90‑1 (MQ=255)                                                                                                                  gatTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                         >  1:1060580/1‑90 (MQ=255)                                                                                                                        ttCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTTTtat  >  1:313285/1‑139 (MQ=255)                                                                                                                            |                                                                                                                                       AACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTTTTAT  >  NC_000913/1544261‑1544518

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