Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	679,874	G→A	G123R (GGG→AGG)	ybeU →	DUF1266 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	679,874	0	G	A	100.0%	197.1 / NA	62	G123R (GGG→AGG)	ybeU	DUF1266 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/31);  total (31/31)

CTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAT  >  NC_000913/679741‑680002                                                                                                                                      |                                                                                                                                 ctctTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGCATCAAAGCCTGGGATTATGTTCGAATGAGGttt                                                                                                                             <  1:1252542/139‑1 (MQ=255) ctctTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGttt                                                                                                                             <  1:569510/139‑1 (MQ=255) ctctTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGttt                                                                                                                             >  1:87577/1‑139 (MQ=255)           gAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACggg                                                                                                                   <  1:173193/139‑1 (MQ=255)            aaGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGt                                                                                                                  >  2:1357349/1‑139 (MQ=255)                 tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCTCAAACCGCTATGTGCTGGGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                             >  1:1030836/1‑139 (MQ=255)                 tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                             <  2:672886/139‑1 (MQ=255)                 tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                             >  2:597232/1‑139 (MQ=255)                 tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                             >  1:1410437/1‑139 (MQ=255)                 tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                             <  2:471786/139‑1 (MQ=255)                 tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                             >  2:1407320/1‑139 (MQ=255)                         aTTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCg                                                                                                                <  2:578754/128‑1 (MQ=255)                         aTTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCg                                                                                                                >  1:578754/1‑128 (MQ=255)                         aTTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTcaac                                                                                                     <  1:597232/139‑1 (MQ=255)                               tGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTgg                                                                                               <  2:781435/139‑1 (MQ=255)                                   ggtggCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGa                                                                                           >  1:510304/1‑139 (MQ=255)                                       gCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGga                                                                                       >  1:1528661/1‑139 (MQ=255)                                        cTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGgaa                                                                                      <  2:855973/139‑1 (MQ=255)                                           ggAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGCAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGgaagaa                                                                                   <  1:747382/139‑1 (MQ=255)                                                          ttATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTg                                                                    >  1:997372/1‑139 (MQ=255)                                                             taCCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCaa                                                                 >  1:292536/1‑139 (MQ=255)                                                              aCCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAAt                                                                >  1:1361489/1‑139 (MQ=255)                                                                  aaTTCGTTGCTCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                                            <  2:1030836/139‑1 (MQ=255)                                                                  aaTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGTAAGAAGAGAGTTTATGGCTGCAATcgcg                                                            <  2:87577/139‑1 (MQ=255)                                                                  aaTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                                            >  1:1452431/1‑139 (MQ=255)                                                                   aTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcgc                                                           <  2:292536/139‑1 (MQ=255)                                                                    ttCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTgg                                                                                               >  1:913345/1‑102 (MQ=255)                                                                    ttCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTgg                                                                                               <  2:913345/102‑1 (MQ=255)                                                                       gTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGt                                                                                                                               <  1:613346/67‑1 (MQ=255)                                                                       gTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGt                                                                                                                               >  2:613346/1‑67 (MQ=255)                                                                        ttGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTa                                                      >  2:29457/1‑139 (MQ=255)                                                                         tGCCCATACCGCTATGTGCTGTAGGGATGGGGGAACCATTGCCTGGGACTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCCCAACAAGTGGCTGACGGAAGAAGAGCGTTTATGGATGCAATTGCGCGTCAAt                                                     >  1:1335211/1‑139 (MQ=255)                                                                          gCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAagaaga                                                                                 >  2:1015557/1‑110 (MQ=255)                                                                          gCCCAAACCGCGATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAagaaga                                                                                 <  1:1015557/110‑1 (MQ=255)                                                                                     tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGc                                                                   <  1:573533/113‑1 (MQ=255)                                                                                     tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGc                                                                   >  2:573533/1‑113 (MQ=255)                                                                                     tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAc                                         >  2:1279094/1‑139 (MQ=255)                                                                                      aTGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAcc                                        <  1:1357349/139‑1 (MQ=255)                                                                                        gtgCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCAt                                      <  2:1435592/139‑1 (MQ=255)                                                                                            tgtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgt                                                                                                              <  1:1137276/63‑1 (MQ=255)                                                                                            tgtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgt                                                                                                              >  2:1137276/1‑63 (MQ=255)                                                                                              tgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAcc                                        <  2:897185/131‑1 (MQ=255)                                                                                              tgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAcc                                        >  1:897185/1‑131 (MQ=255)                                                                                               ggggAAGGGGGAATCAAAGCCTGGGATTATG‑TCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAg                                                                                                  >  1:430743/1‑71 (MQ=37)                                                                                               ggggAAGGGGGAATCAAAGCCTGGGATTATG‑TCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAg                                                                                                  <  2:430743/71‑1 (MQ=37)                                                                                                 ggAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATcac                             >  1:448558/1‑139 (MQ=255)                                                                                                 ggAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATcac                             >  2:1453685/1‑139 (MQ=255)                                                                                                 ggAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATcac                             >  2:1179272/1‑139 (MQ=255)                                                                                                   aaGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAg                           <  1:351555/139‑1 (MQ=255)                                                                                                    aGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                            <  2:243113/121‑1 (MQ=255)                                                                                                    aGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                            >  1:243113/1‑121 (MQ=255)                                                                                                      ggggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                            >  1:397586/1‑119 (MQ=255)                                                                                                      ggggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                            <  2:397586/119‑1 (MQ=255)                                                                                                       gggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTgg                       <  1:1453685/139‑1 (MQ=255)                                                                                                        ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGa                                                                                           >  2:1533639/1‑70 (MQ=255)                                                                                                        ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGa                                                                                           <  1:1533639/70‑1 (MQ=255)                                                                                                          aaTCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATg                    <  1:200891/139‑1 (MQ=255)                                                                                                                aGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTAt              >  2:1051498/1‑139 (MQ=255)                                                                                                                      ggATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTcc        <  2:448558/139‑1 (MQ=255)                                                                                                                           aTGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTa   <  2:510304/139‑1 (MQ=255)                                                                                                                           aTGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTa   <  2:1361489/139‑1 (MQ=255)                                                                                                                            tGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAt  >  1:333764/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                 CTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAT  >  NC_000913/679741‑680002

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