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

ACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176734‑3176990                                                                                                                             |                                                                                                                                     aCAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAAt                                                                                                                        >  1:1053391/1‑139 (MQ=255)     ggCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCAAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGTCGTTAATAAACCGGAATCGcc                                                                                                                    >  1:855603/1‑139 (MQ=255)      gCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCa                                                                                                                   >  2:93940/1‑139 (MQ=255)       cTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCAt                                                                                                                  <  2:700720/139‑1 (MQ=255)       cTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAAGCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCAt                                                                                                                  <  1:573596/139‑1 (MQ=255)           gTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              <  1:151383/139‑1 (MQ=255)            tGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTc                                                                                                             >  1:395060/1‑139 (MQ=255)             gCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  1:93940/139‑1 (MQ=255)             gCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:841338/139‑1 (MQ=255)               ttttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAt                                                                                                          <  1:1300989/139‑1 (MQ=255)                tttGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtg                                                                                                         >  1:1414373/1‑139 (MQ=255)                    tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTGACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  2:555180/1‑132 (MQ=255)                    tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTGACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  1:555180/132‑1 (MQ=255)                    tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCt                                                                                                     <  2:855603/139‑1 (MQ=255)                     gCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTa                                                                                                    <  1:552839/139‑1 (MQ=255)                             gTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    <  2:160022/115‑1 (MQ=255)                             gTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    >  1:160022/1‑115 (MQ=255)                             gTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTAt                                                                                            >  1:1112904/1‑139 (MQ=255)                             gTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGATTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTAt                                                                                            >  2:974131/1‑139 (MQ=255)                             gTGAGTGTCGGTAATTTGGAAAATCCTGGCGCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTGTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCCTGTGCAAAACAGTAt                                                                                            <  1:20291/139‑1 (MQ=255)                              tGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  1:415054/122‑1 (MQ=255)                              tGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  2:415054/1‑122 (MQ=255)                                 gtgtCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCa                                                                                        >  1:552395/1‑139 (MQ=255)                                       gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGc                                                                                                                     >  2:655325/1‑104 (MQ=255)                                       gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGc                                                                                                                     <  1:655325/104‑1 (MQ=255)                                       gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAt                                                                                                          <  1:894495/115‑1 (MQ=255)                                       gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAAc                                                                                                 <  2:830624/124‑1 (MQ=255)                                       gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAAc                                                                                                 >  1:830624/1‑124 (MQ=255)                                       gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGt                                                                                  <  2:395060/139‑1 (MQ=255)                                       gTAATTTGTAAAATACTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTAAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAt                                                                                                          >  2:894495/1‑115 (MQ=255)                                             tGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTa                                                                            >  2:1262703/1‑139 (MQ=255)                                               tAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTaaa                                                                          >  1:1472928/1‑139 (MQ=255)                                                 aaaTCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTaaaaa                                                                        <  2:265284/139‑1 (MQ=255)                                                                    aCCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCttt                                                     >  2:1364486/1‑139 (MQ=255)                                                                             agGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGt                                            <  2:928349/139‑1 (MQ=255)                                                                             agGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGt                                            <  2:1472928/139‑1 (MQ=255)                                                                                 aGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGAt                                        <  2:352546/139‑1 (MQ=255)                                                                                   ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaat                                                                                                                                  >  2:1110859/1‑47 (MQ=255)                                                                                   ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaat                                                                                                                                  <  1:1110859/47‑1 (MQ=255)                                                                                      tAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAAc                                   <  1:248999/139‑1 (MQ=255)                                                                                         cAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                       >  1:122838/1‑68 (MQ=255)                                                                                         cAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                       <  2:122838/68‑1 (MQ=255)                                                                                              tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATaaa                                                                                                                               <  1:120744/39‑1 (MQ=255)                                                                                              tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATaaa                                                                                                                               >  2:120744/1‑39 (MQ=255)                                                                                              tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                                 <  2:1255497/117‑1 (MQ=255)                                                                                              tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgt                                                 >  1:1255497/1‑117 (MQ=255)                                                                                              tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                           <  2:1414373/139‑1 (MQ=255)                                                                                                tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  2:1205929/139‑1 (MQ=255)                                                                                                  cAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGa                       <  2:1112904/139‑1 (MQ=255)                                                                                                       ggTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATc                  >  2:1149424/1‑139 (MQ=255)                                                                                                        gtgtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAAtt                                                              <  2:483212/94‑1 (MQ=255)                                                                                                        gtgtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAAtt                                                              >  1:483212/1‑94 (MQ=255)                                                                                                          gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATa               >  1:1014832/1‑139 (MQ=255)                                                                                                           tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    <  2:220762/37‑1 (MQ=255)                                                                                                           tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGcc                                                                                                                    >  1:220762/1‑37 (MQ=255)                                                                                                           tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAAtt                                                              <  2:1345198/91‑1 (MQ=255)                                                                                                           tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAAtt                                                              >  1:1345198/1‑91 (MQ=255)                                                                                                            ccTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAAc                                   <  2:364995/117‑1 (MQ=255)                                                                                                            ccTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAAc                                   >  1:364995/1‑117 (MQ=255)                                                                                                                  gTTTCACGACGTTAATAAACCGGAATCACCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGCTGCAACTTTTTATTAGGATAATCATACCGCGCtt       >  2:1277502/1‑139 (MQ=255)                                                                                                                    ttCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         >  1:132359/1‑119 (MQ=255)                                                                                                                    ttCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  2:132359/119‑1 (MQ=255)                                                                                                                       acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                <  1:1301657/45‑1 (MQ=255)                                                                                                                       acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                >  2:1301657/1‑45 (MQ=255)                                                                                                                       acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAg  >  2:1090646/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                     ACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176734‑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: ‑