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

TTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGA  >  NC_000913/3176725‑3176988                                                                                                                                      |                                                                                                                                   ttACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaata                                                                                                                               >  1:75825/1‑139 (MQ=255) ttACCCCTAACAGGGCTTCCTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCACACACTGGCGTCACCAGCCAGAGGAAGGGTTAACAGACTTTCCAAATGAAGACCTAAGGATTCACGACGTtaata                                                                                                                               >  2:879497/1‑139 (MQ=255)                   cGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGc                                                                                                             >  2:748116/1‑139 (MQ=255)                      gCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  1:116584/139‑1 (MQ=255)                      gCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                          <  1:298173/139‑1 (MQ=255)                        ttttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAt                                                                                                        <  1:148470/139‑1 (MQ=255)                          ttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAGATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgt                                                                                                      >  1:561295/1‑139 (MQ=255)                           tgtgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaataa                                                                                                                              <  2:214428/114‑1 (MQ=255)                           tgtgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaataa                                                                                                                              >  1:214428/1‑114 (MQ=255)                           tgtgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                     <  2:474398/139‑1 (MQ=255)                           tgtgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                     <  1:702251/139‑1 (MQ=255)                                  aCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                              >  2:795809/1‑139 (MQ=255)                                  aCAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                              <  1:879137/139‑1 (MQ=255)                                     ggTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTa                                                                                           <  2:473539/139‑1 (MQ=255)                                                     ttGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCt                                                                           <  1:796164/139‑1 (MQ=255)                                                          aaaTCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGGTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTaaaaa                                                                      <  2:717184/139‑1 (MQ=255)                                                                      tGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTg                                                          <  2:955047/139‑1 (MQ=255)                                                                             aCCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGAt                                                                 >  2:667915/1‑125 (MQ=255)                                                                             aCCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGAt                                                                 <  1:667915/125‑1 (MQ=255)                                                                                      agGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGt                                          >  1:967717/1‑139 (MQ=255)                                                                                      agGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGt                                          <  1:685328/139‑1 (MQ=255)                                                                                           gggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGt                                                                                >  1:132508/1‑96 (MQ=255)                                                                                           gggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGt                                                                                <  2:132508/96‑1 (MQ=255)                                                                                           gggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATg                                     >  2:638076/1‑139 (MQ=255)                                                                                            ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaa                                                                                                                                 <  2:834905/46‑1 (MQ=255)                                                                                            ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaa                                                                                                                                 >  1:834905/1‑46 (MQ=255)                                                                                             gTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCa                                   >  1:16902/1‑139 (MQ=255)                                                                                              ttAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAGTTTGATGCTTTTCGTCGCGTTGATGCaa                                  >  1:632948/1‑139 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                  <  1:550277/57‑1 (MQ=255)                                                                                                aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                  >  2:550277/1‑57 (MQ=255)                                                                                                 aCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACtt                               <  2:894517/139‑1 (MQ=255)                                                                                                      cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttat                          >  1:70485/1‑139 (MQ=255)                                                                                                      cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttat                          <  2:16902/139‑1 (MQ=255)                                                                                                      cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttat                          >  2:290571/1‑139 (MQ=255)                                                                                                       tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                         <  2:70485/139‑1 (MQ=255)                                                                                                       tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                         <  1:290571/139‑1 (MQ=255)                                                                                                        ttCCAAATGGTGTCCTTAGGTTTCACGACCTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatta                        >  1:39772/1‑139 (MQ=255)                                                                                                         tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttattac                       <  2:659510/139‑2 (MQ=255)                                                                                                         tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       <  2:561295/139‑1 (MQ=255)                                                                                                           cAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGa                     <  1:816574/139‑1 (MQ=255)                                                                                                               tGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAATATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATGAGGATAAt                 >  1:628638/1‑139 (MQ=255)                                                                                                                  tgtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCAt              >  2:9975/1‑139 (MQ=255)                                                                                                                  tgtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCAt              >  1:588608/1‑139 (MQ=255)                                                                                                                   gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                  >  1:129486/1‑38 (MQ=255)                                                                                                                   gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                  <  2:129486/38‑1 (MQ=255)                                                                                                                   gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       <  2:509717/129‑1 (MQ=255)                                                                                                                   gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       >  1:509717/1‑129 (MQ=255)                                                                                                                   gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATa             <  2:632948/139‑1 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       <  2:63031/128‑1 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       >  2:160813/1‑128 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       >  1:63031/1‑128 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       <  2:23879/128‑1 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       >  1:23879/1‑128 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                       <  1:160813/128‑1 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATAc            >  2:593782/1‑139 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATAc            >  2:178443/1‑139 (MQ=255)                                                                                                                    tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATAc            >  1:743898/1‑139 (MQ=255)                                                                                                                              tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCg                                                                             >  1:550208/1‑64 (MQ=255)                                                                                                                              tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCg                                                                             <  2:550208/64‑1 (MQ=255)                                                                                                                              tCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGa  >  1:90917/1‑139 (MQ=255)                                                                                                                                acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                               <  1:501809/92‑1 (MQ=255)                                                                                                                                acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                               >  2:501809/1‑92 (MQ=255)                                                                                                                                      |                                                                                                                                   TTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGA  >  NC_000913/3176725‑3176988

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