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%	182.7 / NA	50	pseudogene (601/726 nt)	efeU	ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/25);  total (25/25)

CCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGT  >  NC_000913/1081950‑1082183                                                                                                               |                                                                                                                            ccATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGt                                                                                                 <  1:197416/139‑1 (MQ=255)                   aaGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTg                                                                               <  1:611609/139‑1 (MQ=255)                    aGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTgg                                                                              <  2:145579/139‑1 (MQ=255)                     gCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                             <  1:468162/139‑1 (MQ=255)                     gCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                             >  1:28963/1‑139 (MQ=255)                       cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttt                                                                           >  1:1060733/1‑139 (MQ=255)                          aTTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttatt                                                                        >  2:458051/1‑139 (MQ=255)                              tgGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGttttt                                                                                                                        <  1:811592/87‑1 (MQ=255)                              tgGAACCACTTTCAGGAAAACGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGttttt                                                                                                                        >  2:811592/1‑87 (MQ=255)                                   ccACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttat                                                                         <  1:1023648/129‑1 (MQ=255)                                   ccACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttat                                                                         >  2:1023648/1‑129 (MQ=255)                                      cTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCg                                                            >  1:345689/1‑139 (MQ=255)                                      cTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCg                                                            <  1:744103/139‑1 (MQ=255)                                        ttCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGgc                                                          <  2:1060733/139‑1 (MQ=255)                                         tCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGgcg                                                         <  1:27932/139‑1 (MQ=255)                                           aGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCt                                                       >  2:1189884/1‑139 (MQ=255)                                                aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTgg                                                                              >  2:703374/1‑111 (MQ=255)                                                aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTCTTTGGCTATCAGGAATCGCCGAGCGTCAGCGAAGTCGCCGTCTgg                                                                              <  1:703374/111‑1 (MQ=255)                                                      ttCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTg                                            <  2:28963/139‑1 (MQ=255)                                                       tCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCa                                                                                               <  1:963574/87‑1 (MQ=255)                                                       tCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCa                                                                                               >  2:963574/1‑87 (MQ=255)                                                       tCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGc                                           >  1:486624/1‑139 (MQ=255)                                                        cGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGct                                          <  1:946618/139‑1 (MQ=255)                                                               aGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAt                                               <  1:942312/127‑1 (MQ=255)                                                               aGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAt                                               >  2:942312/1‑127 (MQ=255)                                                                gTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAtt                                              <  2:214030/127‑1 (MQ=255)                                                                gTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAtt                                              >  1:214030/1‑127 (MQ=255)                                                                gTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacc                                  >  2:340095/1‑139 (MQ=255)                                                                 tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacca                                 >  1:1194519/1‑139 (MQ=255)                                                                    ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgc                              >  2:1183606/1‑139 (MQ=255)                                                                     gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                  >  2:954325/1‑86 (MQ=255)                                                                     gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                  <  1:954325/86‑1 (MQ=255)                                                                         tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTgg                                                                              <  2:190813/86‑1 (MQ=255)                                                                         tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTgg                                                                              >  1:190813/1‑86 (MQ=255)                                                                         tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGct                                          >  1:933981/1‑122 (MQ=255)                                                                         tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGct                                          <  2:933981/122‑1 (MQ=255)                                                                           tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                   >  1:77721/1‑79 (MQ=255)                                                                           tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                   <  2:77721/79‑1 (MQ=255)                                                                                  actcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACgcgc                            <  2:295101/127‑1 (MQ=255)                                                                                  actcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACgcgc                            >  1:295101/1‑127 (MQ=255)                                                                                   ctcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCg               >  2:219173/1‑139 (MQ=255)                                                                                   ctcGCTGTTTGGCACGCTGATGGAAGGGGGTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCg               >  1:482202/1‑139 (MQ=255)                                                                                             ggCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcg     <  1:1189884/139‑1 (MQ=255)                                                                                              gCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc    >  2:567648/1‑139 (MQ=255)                                                                                              gCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc    <  1:365971/139‑1 (MQ=255)                                                                                                aCGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTg                                                                               <  2:4573/62‑1 (MQ=255)                                                                                                aCGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTg                                                                               >  1:4573/1‑62 (MQ=255)                                                                                                aCGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGt  <  1:219173/139‑1 (MQ=255)                                                                                                      aTGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGgcgc                                                        >  2:666958/1‑79 (MQ=255)                                                                                                      aTGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGgcgc                                                        <  1:666958/79‑1 (MQ=255)                                                                                                               |                                                                                                                            CCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGT  >  NC_000913/1081950‑1082183

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