Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,270,523	C→A	L96M (CTG→ATG)	yjbQ →	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,270,523	0	C	A	100.0%	184.4 / NA	53	L96M (CTG→ATG)	yjbQ	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/26);  total (27/26)

CCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCAT  >  NC_000913/4270377‑4270665                                                                                                                                                   |                                                                                                                                               ccTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAatgatg                                                                                                                                              <  2:109336/149‑1 (MQ=255)    ctctGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATggg                                                                                                                                            >  7:100558/1‑148 (MQ=255)       tGACACTTAATGAGAACTGCGATCCCACCGTACGCCTCGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATTAAATCCTCAATGATGGGAACa                                                                                                                                        >  3:67085/1‑149 (MQ=255)            cTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCt                                                                                                                                   <  8:100558/149‑1 (MQ=255)             ttAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCtt                                                                                                                                  >  8:61584/1‑149 (MQ=255)                aTGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAatgatg                                                                                                                                              <  3:88050/134‑1 (MQ=255)                aTGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAatgatg                                                                                                                                              >  4:88050/1‑134 (MQ=255)                      aCTGCGATCCCACCGTACGCCTCGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATTAAATCCTCAATGATGGGAACATCGCTTGTATTGCCg                                                                                                                         <  4:67085/149‑1 (MQ=255)                             tCCCACCGAACGCCACGACATGGAGCGTTTTTTCCTCAGCACCGTTACCGACAACGGAAATTATGAGCATGACTATGAGGCAGCAGACGATATGCCATCTCATATCAAAACCTCAATGATGGGAACATCGCTAGTATTGCCGGTGCAt                                                                                                                   >  7:133716/1‑148 (MQ=255)                                cACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAg                                                                                                               >  3:98373/1‑149 (MQ=255)                                 aCCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATggg                                                                                                                                            <  6:92220/119‑1 (MQ=255)                                 aCCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATggg                                                                                                                                            >  5:92220/1‑119 (MQ=255)                                  ccGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAggg                                                                                                             <  5:96737/149‑1 (MQ=255)                                        gCCACTAAATGGAGCGTTTTTTCCTCCGCTCCGTTCCCGACAACGGAAATTATGTTCATGACTTTGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGATCATCGCTTGTATTGCCGGTGCATAAAGGGCGTAtt                                                                                                       <  8:133716/149‑1 (MQ=255)                                         ccACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAg                                                                                                               >  2:146669/1‑140 (MQ=255)                                         ccACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAg                                                                                                               <  1:146669/140‑1 (MQ=255)                                          caccacATGGCGCGTATTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGCGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCACAATGATGGGAACATCGCTTGTATTGCCTGTGCATAAAGGGCGTAtt                                                                                                       >  2:151779/1‑147 (MQ=255)                                          cACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTg                                                                                                                      <  7:54660/132‑1 (MQ=255)                                          cACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTg                                                                                                                      >  8:54660/1‑132 (MQ=255)                                                  ggAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCa                                                                                             >  8:32586/1‑149 (MQ=255)                                                   gAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCAc                                                                                            <  7:61584/149‑1 (MQ=255)                                                        ttttttCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGc                                                                                       <  7:32586/149‑1 (MQ=255)                                                            ttCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAg                                                                                    <  1:9810/148‑1 (MQ=255)                                                                tccGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCtt                                                                                                                                  <  3:63154/98‑1 (MQ=255)                                                                tccGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCtt                                                                                                                                  >  4:63154/1‑98 (MQ=255)                                                                  cGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACa                                                                                                                                        <  4:59732/90‑1 (MQ=255)                                                                  cGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACa                                                                                                                                        >  3:59732/1‑90 (MQ=255)                                                                       cGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTgg                                                                        >  4:73696/1‑149 (MQ=255)                                                                        gTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGGGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTggg                                                                       <  1:151779/149‑1 (MQ=255)                                                                         ttCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTgggg                                                                      >  1:492/1‑149 (MQ=255)                                                                         ttCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTgggg                                                                      <  4:98373/149‑1 (MQ=255)                                                                            ccGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGaa                                                                   >  2:134104/1‑149 (MQ=255)                                                                                    ggAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCAt                                                           >  5:78815/1‑149 (MQ=255)                                                                                     gAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCAt                                                           >  4:129603/1‑148 (MQ=255)                                                                                     gAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCAt                                                           <  3:129603/148‑1 (MQ=255)                                                                                              aGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGa                                                 <  1:134104/149‑1 (MQ=255)                                                                                               gCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATggg                                                                                                                                            <  1:58522/57‑1 (MQ=255)                                                                                               gCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATggg                                                                                                                                            >  2:58522/1‑57 (MQ=255)                                                                                               gCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTgg                                                                            >  5:29000/1‑121 (MQ=255)                                                                                               gCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTgg                                                                            <  6:29000/121‑1 (MQ=255)                                                                                                       aTGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCg                                         >  2:140514/1‑148 (MQ=255)                                                                                                         gAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGca                                       <  1:140514/148‑1 (MQ=255)                                                                                                             gAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCatcatc                                  <  2:492/149‑1 (MQ=255)                                                                                                               gCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATcgc                                <  3:73696/149‑1 (MQ=255)                                                                                                                       tatGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCAc                                                       >  1:99998/1‑118 (MQ=255)                                                                                                                       tatGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCAc                                                       <  2:99998/118‑1 (MQ=255)                                                                                                                           ccTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGacac                           >  4:50264/1‑142 (MQ=255)                                                                                                                           ccTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGacac                           <  3:50264/142‑1 (MQ=255)                                                                                                                                   tatCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTggggg                                                                     <  6:90458/92‑1 (MQ=255)                                                                                                                                   tatCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTggggg                                                                     >  5:90458/1‑92 (MQ=255)                                                                                                                                       aaaTCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAAt        >  3:1511/1‑149 (MQ=255)                                                                                                                                            cTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATcg                                 >  6:86717/1‑119 (MQ=255)                                                                                                                                            cTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATcg                                 <  5:86717/119‑1 (MQ=255)                                                                                                                                            cTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACa                       >  7:107764/1‑129 (MQ=255)                                                                                                                                            cTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACa                       <  8:107764/129‑1 (MQ=255)                                                                                                                                             tCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCAt  >  2:14406/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               CCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCAT  >  NC_000913/4270377‑4270665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑