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

TTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176725‑3176990                                                                                                                                      |                                                                                                                                     ttACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaata                                                                                                                                 >  1:1096592/1‑139 (MQ=255)    ccccTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAAc                                                                                                                              >  1:1099556/1‑139 (MQ=255)              tgCTTGGTGTTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAGAATCCTGACTCGGGCCTCACCAGCTAGAGGAAGGGTTAACAGGGTTTCCAAATGGTGGCCTTAGGTTTCACGACGTTAATAAACCTGAATCGcc                                                                                                                    <  2:779477/138‑1 (MQ=255)                cTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCAt                                                                                                                  <  2:1096592/139‑1 (MQ=255)                 ttCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCa                                                                                                                   <  2:543342/137‑1 (MQ=255)                 ttCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCa                                                                                                                   >  1:543342/1‑137 (MQ=255)                    gTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              >  1:60265/1‑139 (MQ=255)                      gCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:60265/139‑1 (MQ=255)                                           tgtCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCaa                                                                                       >  2:313596/1‑139 (MQ=255)                                              cGGTAATTTGTAAAATCCTGACTCTGTCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATc                                                                                                                 <  1:605081/110‑1 (MQ=255)                                              cGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATc                                                                                                                 >  2:605081/1‑110 (MQ=255)                                                   aTTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:206726/110‑1 (MQ=255)                                                   aTTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  1:206726/1‑110 (MQ=255)                                                                    tctGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAAtt                                                              <  1:313596/139‑1 (MQ=255)                                                                       ggCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATc                                                                                                                 >  2:929001/1‑85 (MQ=255)                                                                       ggCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATc                                                                                                                 <  1:929001/85‑1 (MQ=255)                                                                        gCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGAt                                                          >  1:641081/1‑139 (MQ=255)                                                                            cACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              <  1:284209/83‑1 (MQ=255)                                                                            cACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              >  2:284209/1‑83 (MQ=255)                                                                                       ggAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGtt                                           >  2:43781/1‑139 (MQ=255)                                                                                            ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGc                                      >  1:113887/1‑139 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    <  2:727621/57‑1 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    >  1:727621/1‑57 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCAt                                                                                                                  <  2:989608/59‑1 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCAt                                                                                                                  >  1:989608/1‑59 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACt                                  >  2:1107368/1‑139 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACt                                  >  2:893808/1‑139 (MQ=255)                                                                                                     gCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTtta                             <  1:893808/139‑1 (MQ=255)                                                                                                      cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttat                            <  1:1107368/139‑1 (MQ=255)                                                                                                       tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                >  1:412687/1‑70 (MQ=255)                                                                                                       tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                <  2:412687/70‑1 (MQ=255)                                                                                                       tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                <  2:59259/70‑1 (MQ=255)                                                                                                       tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                >  1:59259/1‑70 (MQ=255)                                                                                                       tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                           >  2:220295/1‑139 (MQ=255)                                                                                                         tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAATAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  2:962958/139‑1 (MQ=255)                                                                                                         tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGTAACTTTTTATTAg                         <  2:1099556/139‑1 (MQ=255)                                                                                                         tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  1:43781/139‑1 (MQ=255)                                                                                                                  tgtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         >  2:405670/1‑130 (MQ=255)                                                                                                                  tgtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  1:405670/130‑1 (MQ=255)                                                                                                                   gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                                 <  2:1062921/105‑1 (MQ=255)                                                                                                                   gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                                 >  1:1062921/1‑105 (MQ=255)                                                                                                                              tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                           <  2:322936/116‑1 (MQ=255)                                                                                                                              tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                           >  1:322936/1‑116 (MQ=255)                                                                                                                                acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAg  >  2:87288/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                     TTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176725‑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: ‑