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%	221.9 / NA	68	*233* (TAG→TAA)	thiQ	thiamine/thiamine pyrophosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/36);  total (32/36)

AAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGC  >  NC_000913/72098‑72362                                                                                                                                    |                                                                                                                                      aaaGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                <  1:1023586/139‑1 (MQ=255) aaaGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                <  2:1451686/139‑1 (MQ=255)   aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                   <  1:260299/134‑1 (MQ=255)   aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                   >  2:260299/1‑134 (MQ=255)      ggTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATc                                                                                                                           <  2:732210/139‑1 (MQ=255)           gggTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCaa                                                                                                                      <  2:559096/139‑1 (MQ=255)                 aagagaTAGTCTTGTGCGGGTTTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCaa                                                                                                                      <  1:1040589/131‑1 (MQ=255)                 aagagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCaa                                                                                                                      >  2:1040589/3‑133 (MQ=255)                       tAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGc                                                                                                          <  2:704578/139‑1 (MQ=255)                          tCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACt                                                                                                       <  1:726680/139‑1 (MQ=255)                               tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                              <  1:1240910/111‑1 (MQ=255)                               tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                              >  2:1240910/1‑111 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                >  1:1270058/1‑100 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                <  2:1270058/100‑1 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACt                                                                                                       >  2:1157168/1‑125 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACt                                                                                                       <  1:1157168/125‑1 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTca                                                                                         >  1:1637775/1‑139 (MQ=255)                                           gATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                >  1:123440/1‑97 (MQ=255)                                           gATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                <  2:123440/97‑1 (MQ=255)                                             tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCaacaac                                                                                    >  1:348594/1‑139 (MQ=255)                                               cGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTc                                                                                  <  1:1429464/139‑1 (MQ=255)                                                gACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                 <  2:348594/139‑1 (MQ=255)                                                gACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                 <  1:617765/139‑1 (MQ=255)                                                  cGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAtt                                                                               >  1:729062/1‑139 (MQ=255)                                                    cTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTgg                                                                             <  1:885499/139‑1 (MQ=255)                                                       gCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCa                                                                          >  1:983105/1‑139 (MQ=255)                                                           cgTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAcc                                                                      <  1:874043/139‑1 (MQ=255)                                                                           cAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGt                                                                            <  1:1048366/117‑1 (MQ=255)                                                                           cAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGt                                                                            >  2:1048366/1‑117 (MQ=255)                                                                           cAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATgcgc                                                      >  1:104616/1‑139 (MQ=255)                                                                            aaaaCGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTa                                                                                                >  1:253975/1‑96 (MQ=255)                                                                            aaaaCGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTa                                                                                                <  2:253975/96‑1 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTc                                                 >  2:427289/1‑139 (MQ=255)                                                                                 gCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                                <  2:1637775/139‑1 (MQ=255)                                                                                     tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCGCCAATAGTGCCg                                                                                                              >  2:925705/1‑73 (MQ=255)                                                                                     tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCg                                                                                                              <  1:925705/73‑1 (MQ=255)                                                                                     tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAcc                                                                      <  1:1218166/113‑1 (MQ=255)                                                                                     tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAcc                                                                      >  2:1218166/1‑113 (MQ=255)                                                                                     tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGa                                            >  1:1442463/1‑139 (MQ=255)                                                                                      aCCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGAt                                                          <  2:802133/124‑1 (MQ=255)                                                                                      aCCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGAt                                                          >  1:802133/1‑124 (MQ=255)                                                                                               tCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                              >  1:1517176/1‑47 (MQ=255)                                                                                               tCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                              <  2:1517176/47‑1 (MQ=255)                                                                                               tCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGa                                  >  1:1019041/1‑139 (MQ=255)                                                                                               tCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGa                                  >  1:344649/1‑139 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                <  1:1367100/43‑1 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                >  2:1367100/1‑43 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                              >  1:386363/1‑45 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                              <  2:386363/45‑1 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGac                                           >  2:1561450/1‑128 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGac                                           <  1:1561450/128‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                                >  2:1043496/1‑122 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgcg                               >  1:160620/1‑139 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATACCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                                <  1:1043496/122‑1 (MQ=255)                                                                                                        gCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                                <  2:1579811/116‑1 (MQ=255)                                                                                                        gCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                                >  1:1579811/1‑116 (MQ=255)                                                                                                         cGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCg                                                                                                     <  2:918662/62‑1 (MQ=255)                                                                                                         cGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCg                                                                                                     >  1:918662/1‑62 (MQ=255)                                                                                                            tCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATc                     >  2:1499935/1‑139 (MQ=255)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTg                                                                                                                 >  2:747291/1‑46 (MQ=255)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTg                                                                                                                 <  1:747291/46‑1 (MQ=255)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATg                                                         <  1:23679/102‑1 (MQ=255)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATg                                                         >  2:23679/1‑102 (MQ=255)                                                                                                                            aGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCa     <  1:427289/139‑1 (MQ=255)                                                                                                                            aGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCa     <  2:1442463/139‑1 (MQ=255)                                                                                                                            aGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCa     <  1:878503/139‑1 (MQ=255)                                                                                                                              tGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACg   <  2:160620/139‑1 (MQ=255)                                                                                                                               gCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGc  >  1:1014898/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                      AAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGC  >  NC_000913/72098‑72362

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