Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	72,229	C→T	*233* (TAG→TAA)	thiQ ←	thiamine/thiamine pyrophosphate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	72,229	0	C	T	100.0%	213.8 / NA	66	*233* (TAG→TAA)	thiQ	thiamine/thiamine pyrophosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/32);  total (34/32)

CGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACG  >  NC_000913/72095‑72361                                                                                                                                       |                                                                                                                                     cGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                  >  2:95878/1‑139 (MQ=255) cGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                  >  2:700656/1‑139 (MQ=255)      aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                  <  1:1087818/134‑1 (MQ=255)      aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                  >  2:1087818/1‑134 (MQ=255)      aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  1:178003/139‑1 (MQ=255)         ggTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATc                                                                                                                          <  1:976070/139‑1 (MQ=255)              gggTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCaa                                                                                                                     <  2:93524/139‑1 (MQ=255)                  ttAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  1:854137/125‑1 (MQ=255)                  ttAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  2:854137/1‑125 (MQ=255)                      gagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCg                                                                                                             <  1:604450/139‑1 (MQ=255)                         aTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAg                                                                                                          >  2:262119/1‑139 (MQ=255)                                 gtgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGcc                                                                                                  >  1:535704/1‑139 (MQ=255)                                  tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACTCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCt                                                                                                 <  2:764712/139‑1 (MQ=255)                                  tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACt                                                                                                      <  2:586952/134‑1 (MQ=255)                                  tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCt                                                                                                 >  1:1269614/1‑139 (MQ=255)                                  tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATACCCAATAGTGCCGAAGCACt                                                                                                      >  1:586952/1‑134 (MQ=255)                                       gTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  2:869424/104‑1 (MQ=255)                                       gTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  1:869424/1‑104 (MQ=255)                                       gTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            >  2:1058755/1‑139 (MQ=255)                                             tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTcaac                                                                                      <  2:1269614/139‑1 (MQ=255)                                                    aCGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAt                                                                               <  1:262119/139‑1 (MQ=255)                                                    aCGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAt                                                                               <  2:914393/139‑1 (MQ=255)                                                      gCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCa                                                                                                                      <  2:552981/98‑1 (MQ=255)                                                      gCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCa                                                                                                                      >  1:552981/1‑98 (MQ=255)                                                          gCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCa                                                                         >  2:1026035/1‑139 (MQ=255)                                                             gcgTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAc                                                                      <  1:1058755/139‑1 (MQ=255)                                                                         gCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                >  2:379857/1‑117 (MQ=255)                                                                         gCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                <  1:379857/117‑1 (MQ=255)                                                                          ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                  >  2:1081813/1‑66 (MQ=255)                                                                          ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                  <  1:1081813/66‑1 (MQ=255)                                                                          ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  2:963352/69‑1 (MQ=255)                                                                          ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  1:963352/1‑69 (MQ=255)                                                                           cTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATg                                                        >  1:608009/1‑139 (MQ=255)                                                                            tACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            >  1:646908/1‑102 (MQ=255)                                                                            tACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            <  2:646908/102‑1 (MQ=255)                                                                             aCAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             >  1:1207673/1‑68 (MQ=255)                                                                             aCAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  2:1207673/68‑1 (MQ=255)                                                                             aCAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGt                                                                           >  1:281179/1‑118 (MQ=255)                                                                             aCAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGt                                                                           <  2:281179/118‑1 (MQ=255)                                                                                      aCTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             <  1:700656/139‑1 (MQ=255)                                                                                        tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGa                                           >  1:701200/1‑139 (MQ=255)                                                                                                    ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  1:433675/43‑1 (MQ=255)                                                                                                    ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  2:433675/1‑43 (MQ=255)                                                                                                    ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCa                                                                                                        <  2:470166/66‑1 (MQ=255)                                                                                                    ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCa                                                                                                        >  1:470166/1‑66 (MQ=255)                                                                                                    ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgc                               >  1:1011908/1‑139 (MQ=255)                                                                                                    ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgc                               >  2:357045/1‑139 (MQ=255)                                                                                                     gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                <  1:1050288/89‑1 (MQ=255)                                                                                                     gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                >  2:1050288/1‑89 (MQ=255)                                                                                                     gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                       <  1:38893/98‑1 (MQ=255)                                                                                                     gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                       >  2:38893/1‑98 (MQ=255)                                                                                                     gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgcg                              >  1:518451/1‑139 (MQ=255)                                                                                                      aTAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgcgc                             >  1:126867/1‑139 (MQ=255)                                                                                                                cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCTATAGTGCCGAAGCACTCGCCTTACCGCTCAAGAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             <  1:794836/113‑1 (MQ=255)                                                                                                                cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             >  2:794836/1‑113 (MQ=255)                                                                                                                cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATcc                   >  1:318538/1‑139 (MQ=255)                                                                                                                     cgcATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTAc                                                                                              <  2:25916/59‑1 (MQ=255)                                                                                                                     cgcATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTAc                                                                                              >  1:25916/1‑59 (MQ=255)                                                                                                                       cATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgccgc            <  2:535704/139‑1 (MQ=255)                                                                                                                        aTCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAcc                                                                     <  2:532261/81‑1 (MQ=255)                                                                                                                        aTCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAcc                                                                     >  1:532261/1‑81 (MQ=255)                                                                                                                               aGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCa    <  2:126867/139‑1 (MQ=255)                                                                                                                                 tGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACg  <  2:318538/139‑1 (MQ=255)                                                                                                                                 tGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACg  >  1:230652/1‑139 (MQ=255)                                                                                                                                   cATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCACACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgccgc            <  2:92976/127‑1 (MQ=255)                                                                                                                                   cATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCACACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgccgc            >  1:92976/1‑127 (MQ=255)                                                                                                                                       |                                                                                                                                     CGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACG  >  NC_000913/72095‑72361

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