Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	335,970	C→T	Q16* (CAA→TAA)	yahE →	DUF2877 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	335,970	0	C	T	100.0%	189.7 / NA	59	Q16* (CAA→TAA)	yahE	DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/27);  total (32/27)

GGAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCAC  >  NC_000913/335834‑336098                                                                                                                                         |                                                                                                                                 ggAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTaa                                                                                                                                >  2:742016/1‑139 (MQ=255)   aaCGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                              <  1:742016/139‑1 (MQ=255)                agagATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAg                                                                                                                 >  1:290667/1‑139 (MQ=255)                       gcgcAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAACGACAGGTTGAAc                                                                                                                   >  2:594494/1‑130 (MQ=255)                       gcgcAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCAGGGGTAAGGACAGGTTGAAc                                                                                                                   <  1:594494/130‑1 (MQ=255)                              tACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGc                                                                                                   >  1:974010/1‑139 (MQ=255)                               aCTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGt                                                                                                               >  1:174104/1‑126 (MQ=255)                               aCTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGt                                                                                                               <  2:174104/126‑1 (MQ=255)                                       ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                          >  2:823601/1‑139 (MQ=255)                                        cGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                       >  1:420528/1‑125 (MQ=255)                                        cGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                       <  2:420528/125‑1 (MQ=255)                                          cAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGc                                                                                       >  1:280203/1‑139 (MQ=255)                                             gTGCATAATCCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACt                                                                                    >  2:657089/1‑139 (MQ=255)                                                     aCCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTc                                                                            >  1:609954/1‑139 (MQ=255)                                                      ccGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCg                                                                           <  2:609954/139‑1 (MQ=255)                                                         ggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGGGGATTTTCTGGCGCAGCGGGGGTACGGACAGGTTGAACCGGTTGTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCca                                                                        >  1:29167/1‑139 (MQ=255)                                                              cTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGt                                                                   >  1:132539/1‑139 (MQ=255)                                                               ttGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGtt                                                                  >  1:940558/1‑139 (MQ=255)                                                               ttGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGtt                                                                  >  1:442993/1‑139 (MQ=255)                                                                 gCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGc                                                                >  1:964149/1‑139 (MQ=255)                                                                           cacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgctgct                                                      <  2:974010/139‑1 (MQ=255)                                                                           cacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgctgct                                                      <  2:280203/139‑1 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                       <  2:839660/89‑1 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                       >  1:839660/1‑89 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgct                                                         <  1:976604/135‑1 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgct                                                         >  2:976604/1‑135 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                     >  1:248689/1‑139 (MQ=255)                                                                             ccAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCttt                                                    >  1:6668/1‑139 (MQ=255)                                                                                  aaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGaa                                               <  1:823601/139‑1 (MQ=255)                                                                                  aaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGaa                                               <  2:328710/139‑1 (MQ=255)                                                                                      cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                     <  2:373199/81‑1 (MQ=255)                                                                                      cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCaga                                                                                                     >  1:373199/1‑81 (MQ=255)                                                                                      cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGt                                           >  2:429005/1‑139 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                       <  1:466830/77‑1 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                       >  2:466830/1‑77 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                           >  2:343537/1‑121 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                           <  1:343537/121‑1 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTAc                                         >  2:813960/1‑139 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                        <  1:370322/81‑1 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                        >  2:370322/1‑81 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGcc                              >  2:699092/1‑139 (MQ=255)                                                                                                    gcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCa                             <  1:483160/139‑1 (MQ=255)                                                                                                     cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAgg                                                                                                                         >  2:196242/1‑46 (MQ=255)                                                                                                     cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAgg                                                                                                                         <  1:196242/46‑1 (MQ=255)                                                                                                     cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaa                            >  2:455001/1‑139 (MQ=255)                                                                                                      tAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaaa                           <  2:29167/139‑1 (MQ=255)                                                                                                            cggatgcggatTTTCTGCCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTTCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTg                     >  1:505940/1‑139 (MQ=255)                                                                                                              gatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGt                                                                                                <  2:445183/62‑1 (MQ=255)                                                                                                              gatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGt                                                                                                >  1:445183/1‑62 (MQ=255)                                                                                                              gatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTc                   <  1:429005/139‑1 (MQ=255)                                                                                                                  cggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                   >  2:589780/1‑39 (MQ=255)                                                                                                                  cggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                   <  1:589780/39‑1 (MQ=255)                                                                                                                          cTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactca       <  1:813960/139‑1 (MQ=255)                                                                                                                            ggCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcact     <  2:964149/139‑1 (MQ=255)                                                                                                                            ggCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcact     <  2:556232/139‑1 (MQ=255)                                                                                                                            ggCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcact     <  2:6668/139‑1 (MQ=255)                                                                                                                              cgcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactca   <  1:699092/139‑1 (MQ=255)                                                                                                                              cgcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactca   <  2:248689/139‑1 (MQ=255)                                                                                                                               gcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactcac  >  1:971741/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                 GGAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCAC  >  NC_000913/335834‑336098

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