Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,454,187	A→C	pseudogene (815/2646 nt)	yfcU ←	putative outer membrane usher protein YfcU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,454,187	0	A	C	100.0%	87.3 / NA	25	pseudogene (815/2646 nt)	yfcU	putative outer membrane usher protein YfcU
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (13/12);  total (13/12)

TCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCC  >  NC_000913/2454068‑2454308                                                                                                                        |                                                                                                                          tCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAat                                                                                                              >  2:42461/1‑133 (MQ=255)   tGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAatat                                                                                                            <  2:386047/133‑1 (MQ=255)              cACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAgg                                                                                                 >  1:536193/1‑133 (MQ=255)               aCTTTTGCTGTGGTGTGTGCCACGCCCGCAATGTCTGGCGCATAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGt                                                                                                <  2:273559/133‑1 (MQ=255)                         tggtgTGTGCCACGCCGGAAATGTTTGGTGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGcc                                                                                      <  2:189303/133‑1 (MQ=255)                                   cACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGt                                                                            >  2:408598/1‑133 (MQ=255)                                     cGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGttt                                                                          >  2:476761/1‑133 (MQ=255)                                     cGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGttt                                                                          >  1:400863/1‑133 (MQ=255)                                     cGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGttt                                                                          >  2:378241/1‑133 (MQ=255)                                               gTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAAt                                                                <  2:536193/133‑1 (MQ=255)                                                  tGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGAc                                                             >  2:103757/1‑133 (MQ=255)                                                              gcgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAAcgcccgc                                                 <  2:458527/133‑1 (MQ=255)                                                               cgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTGGCCCAGCGACAGGTTGCCTTTTAATAACGGGAAcgcccgcc                                                >  2:164332/1‑133 (MQ=255)                                                               cgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAAcgcccgcc                                                <  1:19417/133‑1 (MQ=255)                                                               cgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAAcgcccgcc                                                <  1:514616/133‑1 (MQ=255)                                                               cgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAAcgcccgcc                                                <  1:378241/133‑1 (MQ=255)                                                                      tGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCata                                         >  1:323733/1‑133 (MQ=255)                                                                           ggCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGc                                    <  1:476761/133‑1 (MQ=255)                                                                             cAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGa                                  <  2:537453/133‑1 (MQ=255)                                                                                  tgtGATCGTTAGTGCGGACACTGCCAGCTACATAATTACAGACGTCGAAAATATCGGAATTGAGGTAATCTTCGCCTAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactcca                             <  2:289917/131‑1 (MQ=255)                                                                                          tCAGTGCTGACACTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGTCGGTAACGCCCGCCCGGCATAATAGCGACTTCACTCCCAtt                     >  2:293755/1‑133 (MQ=255)                                                                                                   acacTGCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTtgtg            >  2:425026/1‑133 (MQ=255)                                                                                                        gCCACCAACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCa       >  2:51377/1‑133 (MQ=255)                                                                                                         ccaccaACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCAt      <  2:323733/133‑1 (MQ=255)                                                                                                         ccaccaACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCAt      <  1:42461/133‑1 (MQ=255)                                                                                                             caACATAATTACAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATcgcc  >  2:602070/1‑133 (MQ=255)                                                                                                                        |                                                                                                                          TCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCC  >  NC_000913/2454068‑2454308

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