Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,595,427	A→G	pseudogene (2561/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,595,427	0	A	G	100.0%	97.2 / NA	27	pseudogene (2561/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/14);  total (13/14)

TTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTACCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGC  >  NC_000913/1595294‑1595557                                                                                                                                      |                                                                                                                                   ttGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGa                                                                                                                               <  1:305973/139‑1 (MQ=255) ttGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGACGATAAGCACGTCGCCCGTCCCGTTTGCCTGa                                                                                                                               <  2:272147/139‑1 (MQ=255)    cgcCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAg                                                                                                                            >  1:444821/1‑139 (MQ=255)           aCCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGc                                                                                                                     <  1:218908/139‑1 (MQ=255)                ttGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTc                                                                                                                <  1:536672/139‑1 (MQ=255)                 tGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCTACGGTTTAGGCACGTCGATAAGCACGTAGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTcc                                                                                                               >  2:371298/1‑139 (MQ=255)                       ccaaCATGGCTGTCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGgtgt                                                                                                         >  1:357416/4‑139 (MQ=255)                             tGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGATTATTCGCCGTGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCAc                                                                                                   >  1:379463/1‑139 (MQ=255)                                   tcgtcgTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCgcgc                                                                                             >  1:254677/1‑139 (MQ=255)                                       cgtGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGAc                                                                                         <  2:357416/139‑1 (MQ=255)                                          gtTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAAc                                                                                      <  2:444821/139‑1 (MQ=255)                                                  aaaTTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTAt                                                                              <  2:254677/139‑1 (MQ=255)                                                        aGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCATGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGc                                                                        <  2:445573/139‑1 (MQ=255)                                                              gtcAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACg                                                                                                                       <  1:263401/86‑1 (MQ=255)                                                              gtcAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACg                                                                                                                       >  2:263401/1‑86 (MQ=255)                                                                         aTTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCcgcg                                                       >  1:395326/1‑139 (MQ=255)                                                                              ccATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGt                                                  <  2:395326/139‑1 (MQ=255)                                                                               cATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACg                                                                         >  1:473782/1‑115 (MQ=255)                                                                               cATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACg                                                                         <  2:473782/115‑1 (MQ=255)                                                                                 tGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCgc                                               >  1:14993/1‑139 (MQ=255)                                                                                                    aGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCg                            >  2:322900/1‑139 (MQ=255)                                                                                                           tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAatga                     >  2:297145/1‑139 (MQ=255)                                                                                                              aTAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATgg                  <  1:311512/139‑1 (MQ=255)                                                                                                                  gCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTAcc              >  2:106500/1‑139 (MQ=255)                                                                                                                    aCGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTc            <  1:297145/139‑1 (MQ=255)                                                                                                                      gTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCa          >  1:232157/1‑139 (MQ=255)                                                                                                                              tcccgtTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGc  <  2:14993/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   TTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTACCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGC  >  NC_000913/1595294‑1595557

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