Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,938	G→A	pseudogene (2037/2037 nt) R6K (AGA→AAA)	rhsE → ydcD →	pseudogene, Rhs family putative immunity protein for RhsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,529,938	0	G	A	100.0%	218.3 / NA	67	pseudogene (2037/2037 nt) R6K (AGA→AAA)	rhsE ydcD	pseudogene, Rhs family putative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/37);  total (30/37)

TTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT  >  NC_000913/1529812‑1530068                                                                                                                               |                                                                                                                                   ttCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGt                                                                                                                        <  1:88020/139‑1 (MQ=255)              aGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAtt                                                                                                           <  2:375772/139‑1 (MQ=255)                cGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAc                                                                                                         >  2:1457179/1‑139 (MQ=255)                 gATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAcc                                                                                                        <  1:1457179/139‑1 (MQ=255)                    gCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCtt                                                                                                     <  1:1337146/139‑1 (MQ=255)                       tGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGAATAAAAGCAGATAATCCAATAAAAGCAAGGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGt                                                                                                  >  1:333707/1‑139 (MQ=255)                               aaaaTATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt                                                                                                            >  1:1636608/1‑121 (MQ=255)                               aaaaTATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt                                                                                                            <  2:1636608/121‑1 (MQ=255)                                   tataACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGttt                                                                                                >  1:110367/1‑129 (MQ=255)                                   tataACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGttt                                                                                                <  2:110367/129‑1 (MQ=255)                                       acacTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTAcc                                                                                  <  1:1608507/139‑1 (MQ=255)                                       acacTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTAcc                                                                                  <  1:55193/139‑1 (MQ=255)                                         acTGGGCAACTCTTCGGGTTAGATGGAATAAAAGCAGATAATCCAATAAAAGCAAGGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCtt                                                                                <  2:333707/139‑1 (MQ=255)                                          cTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg                                                                               <  2:1552957/139‑1 (MQ=255)                                             ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt                                                                                      >  2:1035060/1‑129 (MQ=255)                                             ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt                                                                                      <  1:1035060/129‑1 (MQ=255)                                             ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                            >  1:14087/1‑139 (MQ=255)                                             ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                            >  1:166466/1‑139 (MQ=255)                                             ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                            >  2:469836/1‑139 (MQ=255)                                             ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                            >  1:257951/1‑139 (MQ=255)                                               cAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGttt                                                                          <  2:14087/139‑1 (MQ=255)                                                  ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGAt                                                                       >  2:200249/1‑139 (MQ=255)                                                         ggTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAAtt                                                                <  2:1683022/139‑1 (MQ=255)                                                             aGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCaaa                                                            >  1:1243140/1‑139 (MQ=255)                                                               aTGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCCAGGGAGTTCAGATGCTGCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAcc                                                                                                        <  2:1145243/93‑1 (MQ=255)                                                               aTGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAcc                                                                                                        >  1:1145243/1‑93 (MQ=255)                                                               aTGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAAtt                                                          >  2:782337/1‑139 (MQ=255)                                                                    ataAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAtt                                                                                             >  1:1153003/1‑99 (MQ=255)                                                                    ataAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAtt                                                                                             <  2:1153003/99‑1 (MQ=255)                                                                    ataAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGAc                                                     >  2:64722/1‑139 (MQ=255)                                                                     taAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACa                                                    <  1:456893/139‑1 (MQ=255)                                                                             gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                      >  1:1582893/1‑81 (MQ=255)                                                                             gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                      <  2:1582893/81‑1 (MQ=255)                                                                             gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTtgttgt                                                                                       >  2:1687183/1‑96 (MQ=255)                                                                             gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTtgttgt                                                                                       <  1:1687183/96‑1 (MQ=255)                                                                               tAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTa                                          <  1:782337/139‑1 (MQ=255)                                                                                  tCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTg                                                                                                   <  2:237899/79‑1 (MQ=255)                                                                                  tCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTg                                                                                                   >  1:237899/1‑79 (MQ=255)                                                                                  tCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATgg                                       <  1:64722/139‑1 (MQ=255)                                                                                            gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCa                                                              >  1:1471992/1‑106 (MQ=255)                                                                                            gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCa                                                              <  2:1471992/106‑1 (MQ=255)                                                                                            gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaaaa                               <  2:1207746/137‑1 (MQ=255)                                                                                            gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaaaa                               >  1:1207746/1‑137 (MQ=255)                                                                                                agcaaGGGAGTTCATATGCTTCAAATTATTAAAGGCAAACTTATCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa                                    <  2:137793/128‑1 (MQ=255)                                                                                                agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAtt                                                                                             <  1:386042/71‑1 (MQ=255)                                                                                                agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAtt                                                                                             >  2:386042/1‑71 (MQ=255)                                                                                                agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa                                    <  2:153105/128‑1 (MQ=255)                                                                                                agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa                                    >  1:137793/1‑128 (MQ=255)                                                                                                agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa                                    >  1:153105/1‑128 (MQ=255)                                                                                                 gcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaaa                        <  2:166466/139‑1 (MQ=255)                                                                                                      ggAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTGTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCg                   <  2:1243140/139‑1 (MQ=255)                                                                                                       gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt                  <  2:257951/139‑1 (MQ=255)                                                                                                        aGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGtt                 <  1:469836/139‑1 (MQ=255)                                                                                                          ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt                                                                                 <  1:8638/73‑1 (MQ=255)                                                                                                          ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt                                                                                 >  2:8638/1‑73 (MQ=255)                                                                                                          ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTAtttt                                           <  2:853368/111‑1 (MQ=255)                                                                                                          ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTAtttt                                           >  1:853368/1‑111 (MQ=255)                                                                                                               aTGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGatat          >  1:1427217/1‑139 (MQ=255)                                                                                                                  cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTaa                                                                                                             >  2:1623854/1‑37 (MQ=255)                                                                                                                  cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTaa                                                                                                             <  1:1623854/37‑1 (MQ=255)                                                                                                                  cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATc                                                  >  2:96681/1‑96 (MQ=255)                                                                                                                  cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATc                                                  <  1:96681/96‑1 (MQ=255)                                                                                                                    tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGa                                                      <  1:958986/90‑1 (MQ=255)                                                                                                                    tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGa                                                      >  2:958986/1‑90 (MQ=255)                                                                                                                       aaTTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCt  <  1:200249/139‑1 (MQ=255)                                                                                                                         ttattaAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaaaa                               <  2:557965/108‑1 (MQ=255)                                                                                                                         ttattaAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaaaa                               >  1:557965/1‑108 (MQ=255)                                                                                                                               |                                                                                                                                   TTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT  >  NC_000913/1529812‑1530068

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