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

AACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTTGA  >  NC_000913/335836‑336103                                                                                                                                       |                                                                                                                                      aaCGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                   <  1:1416395/139‑1 (MQ=255)                                     ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                               >  2:237620/1‑139 (MQ=255)                                     ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                               >  1:1342643/1‑139 (MQ=255)                                     ccGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGACTTTGCCAGAGCGGTAAatat                                                                                               >  2:1355125/1‑139 (MQ=255)                                           gTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                   <  2:1301862/97‑1 (MQ=255)                                           gTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                   >  1:1301862/1‑97 (MQ=255)                                           gTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACt                                                                                         >  1:561382/1‑139 (MQ=255)                                                     cGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                               >  1:1142784/1‑123 (MQ=255)                                                     cGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAatat                                                                                               <  2:1142784/123‑1 (MQ=255)                                                       ggAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCca                                                                             >  2:94438/1‑139 (MQ=255)                                                        gAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCcag                                                                            >  2:197477/1‑139 (MQ=255)                                                        gAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCcag                                                                            >  1:954610/1‑139 (MQ=255)                                                              tGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTg                                                                      >  2:672379/1‑139 (MQ=255)                                                                      acacacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCTGGCTCGCCAGCAGTTGCTGACgct                                                              <  1:612158/139‑1 (MQ=255)                                                                      acacacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgct                                                              >  2:644192/1‑139 (MQ=255)                                                                      acacacCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACgct                                                              <  1:783994/139‑1 (MQ=255)                                                                          acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                   <  2:316934/66‑1 (MQ=255)                                                                          acCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAgg                                                                                                                                   >  1:316934/1‑66 (MQ=255)                                                                            cAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATcc                                                                                             <  1:1006235/102‑2 (MQ=255)                                                                            cAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAAGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCg                                                                                             >  2:1006235/1‑102 (MQ=255)                                                                             aGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAg                                                                                                                      <  2:1169990/76‑1 (MQ=255)                                                                             aGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAg                                                                                                                      >  1:1169990/1‑76 (MQ=255)                                                                              gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                      <  2:696929/139‑1 (MQ=255)                                                                              gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                      <  1:197477/139‑1 (MQ=255)                                                                              gAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                      <  2:1342643/139‑1 (MQ=255)                                                                                aaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGaa                                                    <  1:94438/139‑1 (MQ=255)                                                                                aaGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGaa                                                    <  1:644192/139‑1 (MQ=255)                                                                                      gtgtgtgtGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTAc                                              >  2:1264543/1‑139 (MQ=255)                                                                                        gtgtgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGa                                            >  2:60685/1‑139 (MQ=255)                                                                                         tgtgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGAc                                           >  2:977134/1‑139 (MQ=255)                                                                                         tgtgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGAc                                           >  2:18193/1‑139 (MQ=255)                                                                                           tgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                      >  2:1339167/1‑126 (MQ=255)                                                                                           tgtgGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTtgtg                                                      <  1:1339167/126‑1 (MQ=255)                                                                                             tgGGCGCTAACTGCGGAAGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGCGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATg                                       >  2:254189/1‑139 (MQ=255)                                                                                                gcgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATgcgc                                    >  1:1272482/1‑139 (MQ=255)                                                                                                 cgcTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGcc                                   >  2:1012266/1‑139 (MQ=255)                                                                                                   cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                        <  2:1288702/52‑1 (MQ=255)                                                                                                   cTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                        >  1:1288702/1‑52 (MQ=255)                                                                                                    tAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaaa                                <  1:18193/139‑1 (MQ=255)                                                                                                    tAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCaaa                                <  1:672379/139‑1 (MQ=255)                                                                                                     aaCTGCGGATGCGGATTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAAc                               >  2:584919/1‑139 (MQ=255)                                                                                                          cggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                        >  1:908482/1‑45 (MQ=255)                                                                                                          cggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAAc                                                                                                                        <  2:908482/45‑1 (MQ=255)                                                                                                           ggatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGt                         <  2:1272482/139‑1 (MQ=255)                                                                                                            gatgcggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTc                        >  2:950450/1‑139 (MQ=255)                                                                                                                cggatTTTCTGGCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGtt                    >  1:49879/1‑139 (MQ=255)                                                                                                                          ggCGCAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcact          <  1:1264543/139‑1 (MQ=255)                                                                                                                            cgcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATAGCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactca        <  2:358943/139‑1 (MQ=255)                                                                                                                             gcAGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGcactcactcac       >  1:1334017/1‑139 (MQ=255)                                                                                                                               aGCGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACtt     >  2:469351/1‑139 (MQ=255)                                                                                                                                 cGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGca                >  2:162985/1‑126 (MQ=255)                                                                                                                                 cGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGca                <  1:162985/126‑1 (MQ=255)                                                                                                                                 cGGGGGTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTtg   <  2:561382/139‑1 (MQ=255)                                                                                                                                  gggggTAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTtga  <  2:6930/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      AACGGGGCTTTGAAGAGATTGCGCAGTTACTGATTGCCGCAGGTGCATAAACCGGGAGGCTTGCTATCAACACACCAGAAAGACGGTGTGTGTGGGCGCTAACTGCGGATGCGGATTTTCTGGCGCAGCGGGGGCAAGGACAGGTTGAACAGGTCTTTGCCAGAGCGGTAAATATCGCACTCCCGGCTCGCCAGCAGTTGCTGACGCTGCTTTGTGAAGAGTACGACAATGCGCCAAACAGTTGTCGGTTGGCACTCACTCACTTTGA  >  NC_000913/335836‑336103

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