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

TCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACA  >  NC_000913/1081926‑1082187                                                                                                                                       |                                                                                                                                tcgCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGtttt                                                                                                                             <  2:1337670/139‑1 (MQ=255) tcgCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGtttt                                                                                                                             <  1:1562307/139‑1 (MQ=255)    ccgcAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTgg                                                                                                                          <  2:26644/139‑1 (MQ=255)                              cgcgCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCg                                                                                                >  2:1228482/1‑139 (MQ=255)                              cgcgCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCg                                                                                                <  2:726899/139‑1 (MQ=255)                                       cATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                       <  2:1416332/139‑1 (MQ=255)                                        aTGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  1:1279793/139‑1 (MQ=255)                                        aTGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  2:862612/139‑1 (MQ=255)                                        aTGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  2:68600/139‑1 (MQ=255)                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGGAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttt                                                                               >  2:533926/1‑139 (MQ=255)                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAg                                                                                                        >  1:705671/1‑114 (MQ=255)                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAg                                                                                                        <  2:705671/114‑1 (MQ=255)                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttt                                                                               >  1:222133/1‑139 (MQ=255)                                                 gATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                       >  1:225396/1‑129 (MQ=255)                                                 gATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                       <  2:225396/129‑1 (MQ=255)                                                 gATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      >  1:376638/1‑130 (MQ=255)                                                 gATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  2:376638/130‑1 (MQ=255)                                                        gAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttta                                                                              >  2:862200/1‑131 (MQ=255)                                                        gAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCGCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttta                                                                              <  1:862200/131‑1 (MQ=255)                                                           ccACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATc                                                                   >  2:489105/1‑139 (MQ=255)                                                              cTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCg                                                                >  1:538800/1‑139 (MQ=255)                                                                ttCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGgc                                                              >  1:368049/1‑139 (MQ=255)                                                                        aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  1:362233/107‑1 (MQ=255)                                                                        aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      >  2:362233/1‑107 (MQ=255)                                                                                        gTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacc                                      <  1:489105/139‑1 (MQ=255)                                                                                         tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTAt                                                                                                                      <  1:230801/58‑1 (MQ=255)                                                                                         tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTAt                                                                                                                      >  2:230801/1‑58 (MQ=255)                                                                                         tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCg                                                                                          <  1:411886/86‑1 (MQ=255)                                                                                         tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacca                                     >  1:1368410/1‑139 (MQ=255)                                                                                         tGCGGTGCTCTCAACTCACACGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCg                                                                                          >  2:411886/1‑86 (MQ=255)                                                                                           cGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  2:628733/88‑1 (MQ=255)                                                                                           cGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      >  1:628733/1‑88 (MQ=255)                                                                                           cGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGTCGCTGGTGGCATTTGCTCTGCCACCAcg                                   >  2:967641/1‑139 (MQ=255)                                                                                            ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGttttt                                                                                                                            >  2:186940/1‑49 (MQ=255)                                                                                            ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGttttt                                                                                                                            <  1:186940/49‑1 (MQ=255)                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGTTCTGCCACCAcgcg                                 >  2:31437/1‑139 (MQ=255)                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgcg                                 >  2:1529754/1‑139 (MQ=255)                                                                                                ctctCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttat                                                                             >  2:1300752/1‑92 (MQ=255)                                                                                                ctctCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttat                                                                             <  1:1300752/92‑1 (MQ=255)                                                                                                  ctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAggg                            >  2:975535/1‑139 (MQ=255)                                                                                                        tcactcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTCCTCTGCCACCACGCGCAGGGGCGACa                      <  1:1228482/139‑1 (MQ=255)                                                                                                          actcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGc                    >  1:1374999/1‑139 (MQ=255)                                                                                                          actcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGc                    >  1:1403107/1‑139 (MQ=255)                                                                                                               cTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtctc               >  2:1542088/1‑139 (MQ=255)                                                                                                                  tttGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCt            <  2:1374999/139‑1 (MQ=255)                                                                                                                      gCACGCTGATGGGAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc        <  1:533926/139‑1 (MQ=255)                                                                                                                      gCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc        >  2:461178/1‑139 (MQ=255)                                                                                                                       cACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgcg       <  2:368049/139‑1 (MQ=255)                                                                                                                          gCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTg                                                >  2:1494695/1‑95 (MQ=255)                                                                                                                          gCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTg                                                <  1:1494695/95‑1 (MQ=255)                                                                                                                           cTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAAc   <  2:538800/139‑1 (MQ=255)                                                                                                                            tGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  1:355173/55‑1 (MQ=255)                                                                                                                            tGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      >  2:355173/1‑55 (MQ=255)                                                                                                                            tGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACa  >  2:1523079/1‑139 (MQ=255)                                                                                                                              aTGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCt            >  2:131070/1‑127 (MQ=255)                                                                                                                              aTGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCt            <  1:131070/127‑1 (MQ=255)                                                                                                                                ggAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      >  2:1268943/1‑51 (MQ=255)                                                                                                                                ggAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                                      <  1:1268943/51‑1 (MQ=255)                                                                                                                                       |                                                                                                                                TCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACA  >  NC_000913/1081926‑1082187

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