Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,858	C→T	*141* (TAG→TAA)	yqiB ←	DUF1249 protein YqiB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,858	0	C	T	100.0%	120.4 / NA	38	*141* (TAG→TAA)	yqiB	DUF1249 protein YqiB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/16);  total (22/16)

ACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176727‑3176990                                                                                                                                    |                                                                                                                                     aCCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATaaa                                                                                                                               >  1:381689/1‑139 (MQ=255)                  gTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              >  2:54467/1‑139 (MQ=255)                  gTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              <  1:54467/139‑1 (MQ=255)                        ttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  2:844980/1‑135 (MQ=255)                        ttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  1:844980/135‑1 (MQ=255)                                aCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                <  2:381689/139‑1 (MQ=255)                                         tgtCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              <  2:864007/116‑1 (MQ=255)                                         tgtCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              >  1:864007/1‑116 (MQ=255)                                                                           aCCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  1:1038202/1‑84 (MQ=255)                                                                           aCCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:1038202/84‑1 (MQ=255)                                                                                          ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    >  1:401028/1‑61 (MQ=255)                                                                                          ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    <  2:401028/61‑1 (MQ=255)                                                                                              aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACt                                  >  2:1070093/1‑139 (MQ=255)                                                                                                cAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                                 <  1:568240/122‑1 (MQ=255)                                                                                                cAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                                 >  2:568240/1‑122 (MQ=255)                                                                                                 aGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                       >  2:253138/1‑67 (MQ=255)                                                                                                 aGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                       <  1:253138/67‑1 (MQ=255)                                                                                                    cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttat                            <  1:850668/139‑1 (MQ=255)                                                                                                     tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAAt                                                                                                                        <  2:1215192/46‑1 (MQ=255)                                                                                                     tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAAt                                                                                                                        >  1:1215192/1‑46 (MQ=255)                                                                                                        ccAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:1029438/55‑1 (MQ=255)                                                                                                        ccAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  1:1029438/1‑55 (MQ=255)                                                                                                              ggTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCg                                                                                                                >  2:1084130/1‑45 (MQ=255)                                                                                                              ggTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCg                                                                                                                <  1:1084130/45‑1 (MQ=255)                                                                                                              ggTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         >  2:327006/1‑132 (MQ=255)                                                                                                              ggTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  1:327006/132‑1 (MQ=255)                                                                                                                 gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCTCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCaa                                    <  2:150581/118‑1 (MQ=255)                                                                                                                 gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCTCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCaa                                    >  1:150581/1‑118 (MQ=255)                                                                                                                 gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACtttta                              <  1:720903/124‑2 (MQ=255)                                                                                                                 gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACtttta                              >  2:720903/1‑123 (MQ=255)                                                                                                                 gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATa               >  2:1078794/1‑139 (MQ=255)                                                                                                                 gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATa               >  1:1259456/1‑139 (MQ=255)                                                                                                                  tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATAc              >  1:499921/1‑139 (MQ=255)                                                                                                                  tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATAc              >  1:404446/1‑139 (MQ=255)                                                                                                                            tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGa    >  2:359007/1‑139 (MQ=255)                                                                                                                              gcgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAg  >  2:942384/2‑139 (MQ=255)                                                                                                                              acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  1:185979/116‑1 (MQ=255)                                                                                                                              acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         >  2:185979/1‑116 (MQ=255)                                                                                                                                    |                                                                                                                                     ACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176727‑3176990

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