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%	181.2 / NA	57	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 (27/30);  total (27/30)

GGCGATTTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT  >  NC_000913/1529804‑1530062                                                                                                                                       |                                                                                                                             ggCGATTTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGc                                                                                                                          >  1:1604060/1‑139 (MQ=255)                 aaTGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCAtttttt                                                                                                          <  1:375388/139‑1 (MQ=255)                        cGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAc                                                                                                   >  1:615521/1‑139 (MQ=255)                           tGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                <  2:1641543/139‑1 (MQ=255)                           tGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                <  1:1265597/139‑1 (MQ=255)                             cATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCttt                                                                                              >  2:125082/1‑139 (MQ=255)                             cATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                >  1:631578/1‑137 (MQ=255)                             cATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                <  2:631578/137‑1 (MQ=255)                               tGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTaa                                                                                                       <  1:706700/128‑1 (MQ=255)                               tGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTaa                                                                                                       >  2:706700/1‑128 (MQ=255)                               tGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGt                                                                                            >  2:956725/1‑139 (MQ=255)                                    tGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAtt                                                                                       <  2:28447/139‑1 (MQ=255)                                     gAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAttt                                                                                      <  1:125082/139‑1 (MQ=255)                                           tataACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt                                                                                <  2:1604060/139‑1 (MQ=255)                                                     ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt                                                                           <  2:261210/134‑1 (MQ=255)                                                     ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt                                                                           >  1:261210/1‑134 (MQ=255)                                                     ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                      <  2:1317325/139‑1 (MQ=255)                                                     ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                      >  2:429958/1‑139 (MQ=255)                                                     ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                      >  2:900076/1‑139 (MQ=255)                                                     ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                      >  1:1317325/1‑139 (MQ=255)                                                     ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                      >  2:607828/1‑139 (MQ=255)                                                       cAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAt                                                                                        <  1:794223/119‑1 (MQ=255)                                                       cAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAt                                                                                        >  2:794223/1‑119 (MQ=255)                                                          ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAcc                                                                                                  <  1:546223/106‑1 (MQ=255)                                                          ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAcc                                                                                                  >  2:546223/1‑106 (MQ=255)                                                               cGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTCGTTTATTTGTTGTTTACCTTGGGTTTGATAACaa                                                            >  2:375327/1‑139 (MQ=255)                                                                  gTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTc                                                         <  1:900076/139‑1 (MQ=255)                                                                  gTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTc                                                         <  1:1508955/139‑1 (MQ=255)                                                                   ttAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCa                                                        <  1:375327/139‑1 (MQ=255)                                                                       aTGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg                                                                         >  1:608156/1‑118 (MQ=255)                                                                       aTGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg                                                                         <  2:608156/118‑1 (MQ=255)                                                                           tataAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAttt                                                                                      >  2:295845/1‑101 (MQ=255)                                                                           tataAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAttt                                                                                      <  1:295845/101‑1 (MQ=255)                                                                             taAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACa                                              <  1:956725/139‑1 (MQ=255)                                                                                        aaTCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTAt                                   <  2:1229531/139‑1 (MQ=255)                                                                                          tCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATgg                                 >  1:1249315/1‑139 (MQ=255)                                                                                                 aaagcaaagcAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaaa                          <  2:615521/139‑1 (MQ=255)                                                                                                       aagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATgg                                 <  2:1045234/126‑1 (MQ=255)                                                                                                       aagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATgg                                 >  1:1045234/1‑126 (MQ=255)                                                                                                        agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTaa                                                                                                       <  2:1068138/55‑1 (MQ=255)                                                                                                        agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTaa                                                                                                       >  1:1068138/1‑55 (MQ=255)                                                                                                        agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaa                   >  2:1081648/1‑139 (MQ=255)                                                                                                        agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaa                   >  1:1499305/1‑139 (MQ=255)                                                                                                         gcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaaa                  <  2:1499305/139‑1 (MQ=255)                                                                                                         gcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaaa                  <  2:1249315/139‑1 (MQ=255)                                                                                                              ggAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCg             <  2:1265555/139‑1 (MQ=255)                                                                                                               gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt            <  1:429958/139‑1 (MQ=255)                                                                                                                    cAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTgaga       >  1:44276/1‑139 (MQ=255)                                                                                                                         gCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGatatat  >  2:903811/1‑139 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                <  2:357070/44‑1 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                <  2:874144/44‑1 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                >  1:357070/1‑44 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                >  1:874144/1‑44 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt            >  1:539121/1‑128 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt            <  2:539121/128‑1 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCA‑TTTTTTAATTAc                                                                                                   >  1:429877/1‑40 (MQ=255)                                                                                                                          cTTCAAATTATTAAAGGCAAACTTGTCA‑TTTTTTAATTAc                                                                                                   <  2:429877/40‑1 (MQ=255)                                                                                                                                       |                                                                                                                             GGCGATTTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT  >  NC_000913/1529804‑1530062

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