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%	117.3 / NA	37	*141* (TAG→TAA)	yqiB	DUF1249 protein YqiB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/18);  total (19/18)

TTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGA  >  NC_000913/3176724‑3176988                                                                                                                                       |                                                                                                                                   tttACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaat                                                                                                                                >  1:922957/1‑139 (MQ=255)  ttACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaata                                                                                                                               >  1:1089544/1‑139 (MQ=255)  ttACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaata                                                                                                                               >  1:861539/1‑139 (MQ=255)        cTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCgg                                                                                                                         <  1:25172/139‑1 (MQ=255)           aCAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAAt                                                                                                                      >  2:864700/1‑139 (MQ=255)               ggCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                  >  2:904993/1‑139 (MQ=255)               ggCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                  >  2:795071/1‑139 (MQ=255)               ggCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                  <  2:1089544/139‑1 (MQ=255)                       gCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  2:922957/139‑1 (MQ=255)                       gCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  1:795071/139‑1 (MQ=255)                                   aCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                              <  1:471951/139‑1 (MQ=255)                                   aCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                              <  1:904993/139‑1 (MQ=255)                                       gTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTAt                                                                                          <  1:864700/139‑1 (MQ=255)                                             gtCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAAc                                                                                    <  1:452300/139‑1 (MQ=255)                                               cGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCa                                                                                  >  1:1069588/1‑139 (MQ=255)                                                                cTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATaaa                                                                                                                             <  1:545793/79‑1 (MQ=255)                                                                cTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATaaa                                                                                                                             >  2:545793/1‑79 (MQ=255)                                                                   aCTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTaa                                                              >  2:1084901/1‑139 (MQ=255)                                                                         gCCTCACCAGCCAGAGGAAGTGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGAt                                                        >  1:234613/1‑139 (MQ=255)                                                                             cACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  2:366211/85‑1 (MQ=255)                                                                             cACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          >  1:366211/1‑85 (MQ=255)                                                                                          aaGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGa                                       >  1:490509/1‑139 (MQ=255)                                                                                                 aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  1:169604/65‑1 (MQ=255)                                                                                                 aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          >  2:169604/1‑65 (MQ=255)                                                                                                      gCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTtta                           <  2:234613/139‑1 (MQ=255)                                                                                                        tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtg                                                                                                       <  2:172844/61‑1 (MQ=255)                                                                                                        tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtg                                                                                                       >  1:172844/1‑61 (MQ=255)                                                                                                          tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATACACCTGAATCGCCATCTCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTGTTATTAg                       <  2:469343/139‑1 (MQ=255)                                                                                                            cAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGa                     <  2:490509/139‑1 (MQ=255)                                                                                                                     tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  1:897082/45‑1 (MQ=255)                                                                                                                     tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  2:951928/45‑1 (MQ=255)                                                                                                                     tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          >  2:897082/1‑45 (MQ=255)                                                                                                                     tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          >  1:951928/1‑45 (MQ=255)                                                                                                                     tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATAc            >  1:557611/1‑139 (MQ=255)                                                                                                                               tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATg                                     <  2:682342/104‑1 (MQ=255)                                                                                                                               tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATg                                     >  1:682342/1‑104 (MQ=255)                                                                                                                               tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGa  >  2:766733/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                   TTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGA  >  NC_000913/3176724‑3176988

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: ‑