Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,082,060	A→G	pseudogene (601/726 nt)	efeU →	ferrous iron permease (pseudogene)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,082,060	0	A	G	100.0%	97.4 / NA	28	pseudogene (601/726 nt)	efeU	ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/11);  total (17/11)

TTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACAAA  >  NC_000913/1081961‑1082189                                                                                                    |                                                                                                                                  tttCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCC‑TCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGc                                                                                           <  1:127249/139‑1 (MQ=255)          gCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGGTTGGCACGGTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                                   <  1:20217/139‑1 (MQ=255)            cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCg                                                                                            <  1:239892/128‑1 (MQ=255)            cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCg                                                                                            >  2:239892/1‑128 (MQ=255)            cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTc                                                                                       <  2:46864/133‑1 (MQ=255)            cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTc                                                                                       >  1:46864/1‑133 (MQ=255)            cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttt                                                                                 >  1:154135/1‑139 (MQ=255)               aTTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttatt                                                                              >  2:187/1‑139 (MQ=255)                           cTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCg                                                                  >  2:82661/1‑139 (MQ=255)                                                      tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTgg                                                                                                                            <  2:735357/54‑1 (MQ=255)                                                      tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTgg                                                                                                                            >  1:735357/1‑54 (MQ=255)                                                      tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacca                                       >  1:10252/1‑139 (MQ=255)                                                         ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgc                                    >  2:106615/1‑139 (MQ=255)                                                         ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgc                                    >  2:543616/1‑139 (MQ=255)                                                          gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTa                                                                                                                         >  1:331787/1‑53 (MQ=255)                                                          gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTa                                                                                                                         <  2:331787/53‑1 (MQ=255)                                                             ctctCAACTCACTCGCTGTTTGGCACGCTGATGGCAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAg                                >  2:65971/1‑139 (MQ=255)                                                                     tcactcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACa                        <  2:400141/139‑1 (MQ=255)                                                                          cGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacc                                        >  2:362073/1‑118 (MQ=255)                                                                          cGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacc                                        <  1:362073/118‑1 (MQ=255)                                                                          cGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc                   >  1:631222/1‑139 (MQ=255)                                                                            cTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtctc                 >  2:761686/1‑139 (MQ=255)                                                                                 tGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTcc            >  1:566201/1‑139 (MQ=255)                                                                                  ggCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcg           <  2:631222/139‑1 (MQ=255)                                                                                   gCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc          <  1:684841/139‑1 (MQ=255)                                                                                         tGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGt        >  2:644563/1‑135 (MQ=255)                                                                                         tGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGt        <  1:644563/135‑1 (MQ=255)                                                                                           aTGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACaaa  >  2:198475/1‑139 (MQ=255)                                                                                                    |                                                                                                                                  TTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACAAA  >  NC_000913/1081961‑1082189

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: ‑