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%	177.3 / NA	54	Q16* (CAA→TAA)	yahE	DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/26);  total (28/26)

GGAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACT  >  NC_000913/335834‑336095                                                                                                                                         |                                                                                                                              ggAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTaa                                                                                                                             >  2:12788/1‑139 (MQ=255)   aaCGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCTGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAATACGGTGTGTGTGGGCGCTAACTGCGGATGTGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                           <  2:322905/139‑1 (MQ=255)              gAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                <  2:532892/139‑1 (MQ=255)                        cgcAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                           >  2:291662/1‑118 (MQ=255)                        cgcAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                           <  1:291662/118‑1 (MQ=255)                          cAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCa                                                                                                    <  2:472243/139‑1 (MQ=255)                                     tGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAat                                                                                         >  1:895294/1‑139 (MQ=255)                                       ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                       >  2:310851/1‑139 (MQ=255)                                       ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                       >  2:167721/1‑139 (MQ=255)                                             gTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACt                                                                                 >  1:552441/1‑139 (MQ=255)                                               gCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTcc                                                                               >  2:26772/1‑139 (MQ=255)                                                cATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTccc                                                                              <  1:167721/139‑1 (MQ=255)                                                        gggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGcc                                                                      >  1:534183/1‑139 (MQ=255)                                                               ttGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGt                                                                                                            >  1:359172/1‑94 (MQ=255)                                                               ttGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGt                                                                                                            <  2:359172/94‑1 (MQ=255)                                                                 gCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAAGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGc                                                             <  1:247904/139‑1 (MQ=255)                                                                      cAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACg                                                        <  1:310851/139‑1 (MQ=255)                                                                         cacacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCACACTCCCGGCTCGCCAGCAGTTGCTGACgctg                                                     >  1:791346/1‑139 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAgg                                                                                                                      <  2:761639/71‑1 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAgg                                                                                                                      >  1:761639/1‑71 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                  >  2:525603/1‑139 (MQ=255)                                                                            acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCtt                                                  >  2:585803/1‑139 (MQ=255)                                                                               aGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgt                                               >  1:666048/1‑139 (MQ=255)                                                                                gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                              <  1:585803/139‑1 (MQ=255)                                                                                gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                              <  2:344065/139‑1 (MQ=255)                                                                                  aaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGTGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGaa                                            <  2:285641/139‑1 (MQ=255)                                                                                   aaagacggtgtgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                              >  1:332927/8‑136 (MQ=255)                                                                                   aaagacggtgtgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                              <  2:332927/129‑1 (MQ=255)                                                                                      cGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGt                                        >  2:514413/1‑139 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGc                                                                                                      <  2:389427/75‑1 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGc                                                                                                      >  1:389427/1‑75 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCagag                                                                                                 <  2:239036/80‑1 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCagag                                                                                                 >  1:239036/1‑80 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGt                                                                                             >  2:972398/1‑84 (MQ=255)                                                                                        gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGt                                                                                             <  1:972398/84‑1 (MQ=255)                                                                                           tgtgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGAc                                   >  2:482655/1‑139 (MQ=255)                                                                                             tgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACaa                                 >  2:852007/1‑139 (MQ=255)                                                                                               tgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGc                                                                           >  1:1022053/1‑95 (MQ=255)                                                                                               tgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGc                                                                           <  2:1022053/95‑1 (MQ=255)                                                                                                   cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGcc                           >  2:969990/1‑139 (MQ=255)                                                                                                     cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaa                         >  2:572360/1‑139 (MQ=255)                                                                                                     cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCa                          <  2:867037/138‑1 (MQ=255)                                                                                                     cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCa                          >  1:867037/1‑138 (MQ=255)                                                                                                      tAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaaa                        <  2:666048/139‑1 (MQ=255)                                                                                                       aaCTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAAc                       >  1:278472/1‑139 (MQ=255)                                                                                                            cggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                              <  1:834424/111‑1 (MQ=255)                                                                                                            cggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                              >  2:834424/1‑111 (MQ=255)                                                                                                            cggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTg                  >  1:1040335/1‑139 (MQ=255)                                                                                                              gatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTc                <  2:534183/139‑1 (MQ=255)                                                                                                                      ttttCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGca        <  1:572360/139‑1 (MQ=255)                                                                                                                      ttttCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGca        <  1:525603/139‑1 (MQ=255)                                                                                                                          cTGGCGCAGCGTGGGTAAGGACAGGTTGAACAGGTCTTTGCTAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTTCTGACGCTGCTTTTTGCAGACTACGACAATTCGCCAAACAGTTGTCGGTTGGcactaa    <  2:681284/139‑3 (MQ=255)                                                                                                                          cTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactca    <  2:552441/139‑1 (MQ=255)                                                                                                                            ggCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcact  <  1:852007/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                              GGAACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACT  >  NC_000913/335834‑336095

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