Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,103,378	T→C	pseudogene (12/446 nt)	wbbL ←	interrupted rhamnosyltransferase WbbL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,103,378	0	T	C	100.0%	30.4 / NA	10	pseudogene (12/446 nt)	wbbL	interrupted rhamnosyltransferase WbbL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/8);  total (2/8)

ATTAAGATTTTCGAGTAATTTTTTGATGTAGTCTTCATGTCCGTGGGAAACGATTATTATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAACCCAACTAATACATTTTCATTAC  >  NC_000913/2103321‑2103452                                                          |                                                                           aTTAAGATTTTCGAGTAATTTTTTGATGTAGTCTTCATGTCCGTGGGAAACGATTATCATATATACCATTTCAATGttctt                                                     <  2:372718/81‑1 (MQ=255)                     ttttGATGTAGTCTTCATGTCCGTGGGAAACGATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGt                                 >  1:372736/1‑81 (MQ=255)                          aTGTAGTCTTCATGTCCGTGGGAAACGATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTtcatc                            <  1:372719/81‑1 (MQ=255)                          aTGTAGTCTTCATGTCCGTGGGAAACGATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTtcat                             <  2:372720/80‑1 (MQ=255)                              aGTCTTCATGTCCGTGGGAAACGATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAAc                        >  2:372745/1‑81 (MQ=255)                               gTCTTCATGTCCGTGGGAAACGATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAAcc                       <  1:372721/81‑1 (MQ=255)                                    cATGTCCGTGGGAAACGATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAACCCAAc                   <  1:372722/80‑1 (MQ=255)                                                  aCGATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAACCCAACTAATACATTTTCAtt    <  1:372723/81‑1 (MQ=255)                                                    gATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAACCCAACTAATACATTTTCATTAc  <  1:372724/81‑1 (MQ=255)                                                    gATTATCATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAACCCAACTAATACATTTTCATTAc  <  2:372725/81‑1 (MQ=255)                                                          |                                                                           ATTAAGATTTTCGAGTAATTTTTTGATGTAGTCTTCATGTCCGTGGGAAACGATTATTATATATACCATTTCAATGTTCTTCAGTAATAAAATTAACTAGTTCATCAAACCCAACTAATACATTTTCATTAC  >  NC_000913/2103321‑2103452

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑