Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,081,596	T→C	pseudogene (241/830 nt)	efeU →	inactive ferrous iron permease EfeU, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,081,596	0	T	C	100.0%	59.7 / NA	21	pseudogene (241/830 nt)	efeU	inactive ferrous iron permease EfeU, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/15);  total (6/15)

GTGATGTGGATTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGTGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAATCATCATGGCTGGGCGCTGGTGATGATGG  >  NC_000913/1081463‑1081721                                                                                                                                       |                                                                                                                              gtgATGTGGATTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATc                                                                                                                           <  1:3802/138‑1 (MQ=255) gtgATGTGGATTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATc                                                                                                                           >  2:3802/1‑138 (MQ=255)     tgtgGATTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTAcc                                                                                                                     >  1:38291/1‑140 (MQ=255)         gATTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGGATCGTGGTGGGGAT‑GGCCGGGAGGGTCCTTACCTggg                                                                                                                 >  1:186230/1‑139 (MQ=255)          aTTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGGT‑CGCCGTGGCGGTCCTTACCTgggt                                                                                                                >  2:108337/1‑138 (MQ=255)          aTTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGAt                                                                                                                >  2:115260/1‑140 (MQ=255)          aTTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCCTGGCGATCCTTACCTGGAt                                                                                                                >  1:129326/1‑140 (MQ=255)                   ttgcttgcCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGGTCCGCCG‑GGCGGTCCTTTCCTGGATGGgtttttt                                                                                                       >  1:474052/1‑137 (MQ=255)                    tgcttgcCGCTGCGTTGTGCCTGGGCTTGGGTATCCTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTgg                                                                                                      <  2:382430/140‑1 (MQ=255)                             cTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCaaa                                                                                             <  2:394862/140‑1 (MQ=255)                                             cTTGGGTATCTTCAGTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCa                                                                             <  1:614000/140‑1 (MQ=255)                                               tGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCaaa                                                                           >  1:385565/1‑140 (MQ=255)                                                gggTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAg                                                                          <  2:644509/140‑1 (MQ=255)                                                 ggTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGt                                                                         <  2:113746/140‑1 (MQ=255)                                                 ggTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGt                                                                         >  1:515675/1‑140 (MQ=255)                                                     tCTTCATTAACGAAAGCACCGGGGATTTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCaa                                                                     <  1:248524/140‑1 (MQ=255)                                                      cTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGCGGGGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAATGCAAc                                                                    >  2:511076/1‑140 (MQ=255)                                                           ttAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGCGGCGATCCTTACCTTGGTAGGTTTTTGGGTGCGAAAAGTGTCGCACAACTTCAAAGAGGAACTggga                                                               >  2:304301/1‑138 (MQ=255)                                                                      accGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGa                                                    <  2:656831/140‑1 (MQ=255)                                                                                ttCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAAAGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTg                                          <  2:385565/140‑1 (MQ=255)                                                                                     cAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCg                                     <  2:117650/140‑1 (MQ=255)                                                                                         aaGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGa                                 <  2:578569/140‑1 (MQ=255)                                                                                                 gAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAAtcatca                         <  1:115260/140‑1 (MQ=255)                                                                                                  aaCTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGGTCCTTACCTGGATGGTTTTCTGGATGGGCAAAGTGTCGCGCAACGTCAAAATGCAACTGGAACAGGCAGTCGATAGCGCATTTCAGCGGGGAAATcatcat                        >  2:134847/1‑140 (MQ=255)                                                                                                    cTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAATCATCATgg                      <  2:69925/140‑1 (MQ=255)                                                                                                      gTTTGAAGGTATCGTGGCGGTGAT‑CGCCGGGGCGGTCCTTACCTGGGTGGGTTTCTGGGTGCGGAAAGGGTCGCGCAACGGCAAAGTGCAAATGGGACAGGGAGGCGGTAGCGCCTTGCAGCGGGGGAATTATCATGGcg                    >  2:462818/1‑139 (MQ=255)                                                                                                          gAAGGTATCGTGGCGGTCAT‑CGCCGTGGCGATCCTCACCCGGATGGTTGTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAATCATCATGGCTGGgc                <  2:186230/140‑1 (MQ=255)                                                                                                                 tCGTGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGGTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAATCATCATGGCTGGGCGCTgggg         >  1:596007/1‑138 (MQ=255)                                                                                                                    tGGCGGTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAATCATCATGGCTGGGCGCTGGtgatg      <  2:515675/140‑1 (MQ=255)                                                                                                                        ggTGAT‑CGCCGTGGCGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGGGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAATCATCATGGCTGGGCGCTGGGGATGATgg  >  2:350603/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                              GTGATGTGGATTGGCGTGTTGCTTGCCGCTGCGTTGTGCCTGGGCTTGGGTATCTTCATTAACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGAT‑CGCCGTGGTGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAGCGTGGAAATCATCATGGCTGGGCGCTGGTGATGATGG  >  NC_000913/1081463‑1081721

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