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%	146.7 / NA	46	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 (22/24);  total (22/24)

TGGCGATTTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA  >  NC_000913/1529803‑1530071                                                                                                                                        |                                                                                                                                      tgGCGATTTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAgg                                                                                                                                    >  1:987246/1‑139 (MQ=255)                              cATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCttt                                                                                                       >  1:354946/1‑139 (MQ=255)                              cATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCttt                                                                                                       >  1:45296/1‑139 (MQ=255)                                              taACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTa                                                                                       >  2:399059/1‑139 (MQ=255)                                                    tGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCTAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTgg                                                                                 >  2:787566/1‑139 (MQ=255)                                                    tGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                         >  1:1367935/1‑115 (MQ=255)                                                    tGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt                                                                                                         <  2:1367935/115‑1 (MQ=255)                                                      ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt                                                                                                               <  1:1302363/107‑1 (MQ=255)                                                      ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt                                                                                                               >  2:1302363/1‑107 (MQ=255)                                                      ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt                                                                               >  2:881319/1‑139 (MQ=255)                                                           ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTtgttgt                                                                                          >  2:476374/1‑123 (MQ=255)                                                           ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTtgttgt                                                                                          <  1:476374/123‑1 (MQ=255)                                                           ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTGGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGAt                                                                          <  1:841320/139‑1 (MQ=255)                                                           ctctTCGGGTTAGATGGTATAAAAGCAGACAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg                                                                                  <  1:372764/131‑1 (MQ=255)                                                           ctctTCGGGTTAGATGGTATAAAAGCAGACAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg                                                                                  >  2:372764/1‑131 (MQ=255)                                                            tctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATa                                                                         >  1:295522/1‑139 (MQ=255)                                                              ttCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCTAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAAc                                                                       <  1:787566/139‑1 (MQ=255)                                                              ttCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTAcc                                                                                     >  2:1320991/1‑125 (MQ=255)                                                              ttCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTAcc                                                                                     <  1:1320991/125‑1 (MQ=255)                                                                cGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCAtttt                                                                                                                     <  2:1018106/91‑1 (MQ=255)                                                                cGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCAtttt                                                                                                                     >  1:1018106/1‑91 (MQ=255)                                                                    ttAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCa                                                                 <  2:987246/139‑1 (MQ=255)                                                                              taAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACa                                                       <  1:399059/139‑1 (MQ=255)                                                                              taAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACa                                                       <  2:37501/139‑1 (MQ=255)                                                                                      gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACt                                                                                                                              <  2:636531/60‑1 (MQ=255)                                                                                      gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACt                                                                                                                              >  1:636531/1‑60 (MQ=255)                                                                                      gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATaa                                                                        >  2:1117732/1‑114 (MQ=255)                                                                                      gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATaa                                                                        <  1:1117732/114‑1 (MQ=255)                                                                                      gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTAttt                                               >  2:1065673/1‑139 (MQ=255)                                                                                      gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTAttt                                               >  1:246102/1‑139 (MQ=255)                                                                                       aTAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTAtttt                                              <  2:295522/139‑1 (MQ=255)                                                                                         aaTCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTAt                                            <  1:672353/139‑1 (MQ=255)                                                                                               aTAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCt                                                           <  2:574179/118‑1 (MQ=255)                                                                                               aTAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCt                                                           >  1:574179/1‑118 (MQ=255)                                                                                                tAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATa                                     <  2:354946/139‑1 (MQ=255)                                                                                                  aaagcaaagcAAGTGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaaa                                   <  1:745158/139‑1 (MQ=255)                                                                                                          gcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaaa                           <  1:881319/139‑1 (MQ=255)                                                                                                                 aGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGtt                    <  2:246102/139‑1 (MQ=255)                                                                                                                   ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa                                       >  1:741243/1‑118 (MQ=255)                                                                                                                   ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa                                       <  2:741243/118‑1 (MQ=255)                                                                                                                   ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTga                  >  2:1326630/1‑139 (MQ=255)                                                                                                                    tCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTgag                 <  1:1326630/139‑1 (MQ=255)                                                                                                                             tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACTCCAAAGAGCGTTGAGATATATCtt        >  2:762685/1‑139 (MQ=255)                                                                                                                             tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaaa                           <  1:479966/120‑1 (MQ=255)                                                                                                                             tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaaa                           >  2:479966/1‑120 (MQ=255)                                                                                                                                   tattaAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaa  <  1:1071291/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      TGGCGATTTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA  >  NC_000913/1529803‑1530071

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