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%	141.7 / NA	37	pseudogene (2561/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/17);  total (20/17)

TTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTACCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTC  >  NC_000913/1595293‑1595562                                                                                                                                       |                                                                                                                                        tttGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGAATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTg                                                                                                                                     >  1:600405/1‑139 (MQ=255)                 ttGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTc                                                                                                                     >  2:892131/1‑139 (MQ=255)                 ttGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTc                                                                                                                     <  2:426729/139‑1 (MQ=255)                           aCATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGt                                                                                                           >  2:740627/1‑139 (MQ=255)                                         gtgtTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACa                                                                                             >  2:488719/1‑139 (MQ=255)                                           gtTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAAc                                                                                           >  2:613057/1‑139 (MQ=255)                                                   aaaTTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTAt                                                                                   <  1:488719/139‑1 (MQ=255)                                                      ttGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGa                                                                                <  1:892131/139‑1 (MQ=255)                                                         aGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGc                                                                             <  2:600405/139‑1 (MQ=255)                                                              cgtcAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGata                                                                        >  2:343266/1‑139 (MQ=255)                                                                          aTTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCcgcg                                                            >  2:381140/1‑139 (MQ=255)                                                                                cATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGt                                                       >  1:558880/1‑138 (MQ=255)                                                                                cATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGt                                                       <  2:558880/138‑1 (MQ=255)                                                                                cATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTc                                                      >  1:420585/1‑139 (MQ=255)                                                                                  tGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCtt                                                                                                                       <  2:843276/72‑1 (MQ=255)                                                                                  tGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCtt                                                                                                                       >  1:843276/1‑72 (MQ=255)                                                                                     gATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCgcagc                                                 <  1:592219/139‑1 (MQ=255)                                                                                          ttCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGaa                                            >  1:453650/1‑139 (MQ=255)                                                                                                  tttAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTa                                    >  1:728663/1‑139 (MQ=255)                                                                                                     aGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGa                                                                                >  1:579718/1‑92 (MQ=255)                                                                                                     aGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGa                                                                                <  2:579718/92‑1 (MQ=255)                                                                                                       gCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGaa                               <  2:453650/139‑1 (MQ=255)                                                                                                        cACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGaaa                              <  2:728663/139‑1 (MQ=255)                                                                                                           gTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGc                                                                                                                          >  2:114490/1‑44 (MQ=255)                                                                                                           gTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGc                                                                                                                          <  1:114490/44‑1 (MQ=255)                                                                                                            tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTg                                                                           <  2:359183/90‑1 (MQ=255)                                                                                                            tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTg                                                                           >  1:359183/1‑90 (MQ=255)                                                                                                            tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGt                                                       <  2:492685/110‑1 (MQ=255)                                                                                                            tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGt                                                       >  1:492685/1‑110 (MQ=255)                                                                                                            tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAatga                          >  1:812691/1‑139 (MQ=255)                                                                                                                     aCGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATgg                       <  2:126035/133‑1 (MQ=255)                                                                                                                     aCGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATgg                       >  1:126035/1‑133 (MQ=255)                                                                                                                      cGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTcc                >  1:689366/1‑139 (MQ=255)                                                                                                                         cGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGc                                                                             >  2:323955/1‑75 (MQ=255)                                                                                                                         cGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGc                                                                             <  1:323955/75‑1 (MQ=255)                                                                                                                                    tttGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTc  <  2:689366/139‑1 (MQ=255)                                                                                                                                    tttGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTc  <  1:381140/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        TTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTACCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTC  >  NC_000913/1595293‑1595562

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