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

ACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTACCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTC  >  NC_000913/1595319‑1595562                                                                                                             |                                                                                                                                        aCATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGt                                                                                                           >  2:386775/1‑139 (MQ=255) aCATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGt                                                                                                           >  1:785604/1‑139 (MQ=255)     ggCTATCGCCGTGTTCCGGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCAcc                                                                                                       <  2:678483/139‑1 (MQ=255)                tgtTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACaa                                                                                            >  1:66016/1‑139 (MQ=255)                       gCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGa                                                                                                                                    <  2:787172/92‑1 (MQ=255)                       gCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGa                                                                                                                                    >  1:787172/1‑92 (MQ=255)                          aaTTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATc                                                                                  <  1:490510/139‑1 (MQ=255)                            ttGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAg                                                                                        >  2:768876/1‑131 (MQ=255)                            ttGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAg                                                                                        <  1:768876/131‑1 (MQ=255)                            ttGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGa                                                                                <  2:642330/139‑1 (MQ=255)                              gAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAATTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACg                                                                              <  2:692692/139‑1 (MQ=255)                               aGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGc                                                                             <  1:343053/139‑1 (MQ=255)                                    cgtcAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCttt                                                                                                                      <  1:441038/93‑1 (MQ=255)                                    cgtcAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCttt                                                                                                                      >  2:441038/1‑93 (MQ=255)                                    cgtcAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGc                                                                             <  2:729352/134‑1 (MQ=255)                                    cgtcAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGc                                                                             >  1:729352/1‑134 (MQ=255)                                       cAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGatatat                                                                     <  2:66016/139‑1 (MQ=255)                                          agGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTAt                                                                                   <  1:937679/122‑1 (MQ=255)                                          agGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTAt                                                                                   >  2:937679/1‑122 (MQ=255)                                                aTTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCcgcg                                                            >  2:1023613/1‑139 (MQ=255)                                                 ttCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGa                                                           >  2:748237/1‑139 (MQ=255)                                                       aTGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCg                                                     <  2:903922/139‑1 (MQ=255)                                                        tGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCgc                                                    >  2:178060/1‑139 (MQ=255)                                                          ggATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCgcag                                                  <  2:542867/139‑1 (MQ=255)                                                           gaacATTCCACCGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCgcagc                                                 >  2:651968/4‑139 (MQ=255)                                                           gATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCgcagc                                                 <  1:178060/139‑1 (MQ=255)                                                                     cGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGa                                                                                >  2:866092/1‑98 (MQ=255)                                                                     cGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGa                                                                                <  1:866092/98‑1 (MQ=255)                                                                     cGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATcc                                       >  1:937503/1‑139 (MQ=255)                                                                          tAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATc                                  <  1:790340/139‑1 (MQ=255)                                                                                 gTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGcc                                                                                                                         <  2:491918/45‑1 (MQ=255)                                                                                 gTCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGcc                                                                                                                         >  1:491918/1‑45 (MQ=255)                                                                                  tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCcg                                                                                                     <  1:1028892/64‑1 (MQ=255)                                                                                  tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCcg                                                                                                     >  2:1028892/1‑64 (MQ=255)                                                                                  tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGt                                                       <  2:1018006/110‑1 (MQ=255)                                                                                  tCGATAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGt                                                       >  1:1018006/1‑110 (MQ=255)                                                                                      tAAGCACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGt                      >  2:79444/1‑139 (MQ=255)                                                                                          cACGTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCt                  <  1:1023613/139‑1 (MQ=255)                                                                                             gTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCa               >  1:189173/1‑139 (MQ=255)                                                                                             gTCGCCCGTCCCGTTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCa               <  2:189173/139‑1 (MQ=255)                                                                                                     tcccgtTTGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAAGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGc       >  1:369692/1‑139 (MQ=255)                                                                                                          tttGCCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTc  >  1:547676/1‑139 (MQ=255)                                                                                                             |                                                                                                                                        ACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCATTCCACGGTTTAGGCACGTCGATAAGCACGTCGCCCGTCCCGTTTACCTGATAGTTACGGCCTTTCCAGGTGTAGTCACCCGCGCCGACAACCAGCGTATCGACGCTGATATATCCCTCCGCGAGCGTCGCAGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTC  >  NC_000913/1595319‑1595562

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