Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,316	T→G	Q402H (CAA→CAC)	thiP ←	thiamine/thiamine pyrophosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,316	0	T	G	100.0%	187.8 / NA	53	Q402H (CAA→CAC)	thiP	thiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (23/30);  total (23/30)

TTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAT  >  NC_000913/73169‑73463                                                                                                                                                    |                                                                                                                                                    ttttAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAgtgt                                                                                                                                                  >  2:339909/1‑151 (MQ=255)  tttAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAgtgt                                                                                                                                                  <  1:346114/150‑1 (MQ=255)      aGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAg                                                                                                                                             <  2:438576/151‑1 (MQ=255)         gTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGcc                                                                                                                                          <  1:123760/151‑1 (MQ=255)               aGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAgt                                                                                                                                    <  2:494271/151‑1 (MQ=255)                    tCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGtgttg                                                                                                                                <  1:598568/150‑1 (MQ=255)                    tCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGtgtt                                                                                                                                 <  1:254434/149‑1 (MQ=255)                     cAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTgttgtt                                                                                                                              <  1:143833/151‑1 (MQ=255)                               cGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTaa                                                                                                                     >  1:408230/1‑150 (MQ=255)                                gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTaaa                                                                                                                    >  2:732147/1‑150 (MQ=255)                                gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTaaa                                                                                                                    >  1:151303/1‑150 (MQ=255)                                gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTaaa                                                                                                                    <  2:3812/150‑1 (MQ=255)                                 aCTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTaaaaa                                                                                                                  <  2:670555/151‑1 (MQ=255)                                 aCTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTaaaaa                                                                                                                  <  2:477808/151‑1 (MQ=255)                                    gACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAaagaa                                                                                                               <  2:148244/151‑1 (MQ=255)                                    gACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAaaga                                                                                                                <  1:25759/150‑1 (MQ=255)                                    gACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAACAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAaagaa                                                                                                               <  2:449393/151‑1 (MQ=255)                                         aaCATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCg                                                                                                           <  1:407808/150‑1 (MQ=255)                                                 gTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGc                                                                                                  <  2:98807/151‑1 (MQ=255)                                                             gATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCgg                                                                                       <  1:77488/150‑1 (MQ=255)                                                              atatCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCa                                                                                     <  1:173874/151‑1 (MQ=255)                                                              atatCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGACAGCAGAGTGTGGCAGGCCGTTAGTGTAGTTGAGCAGTAAAAAGAAGCCGGTAGCCGGCCCACTGCCCGGCa                                                                                     <  2:123797/151‑1 (MQ=255)                                                                atCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCAt                                                                                    <  2:480928/150‑1 (MQ=255)                                                                  cgcgCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGc                                                                                 >  2:521576/1‑151 (MQ=255)                                                                  cgcgCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGc                                                                                 <  1:23306/151‑1 (MQ=255)                                                                       aTCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAgg                                                                             >  1:168596/1‑150 (MQ=255)                                                                             ttttCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAAc                                                                       >  1:787517/1‑150 (MQ=255)                                                                                  cAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCg                                                                 >  2:861823/1‑151 (MQ=255)                                                                                  cAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCg                                                                 <  1:756662/151‑1 (MQ=255)                                                                                     cACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTc                                                              >  2:684447/1‑151 (MQ=255)                                                                                       cTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCAt                                                            >  2:563688/1‑151 (MQ=255)                                                                                       cTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCAt                                                            >  1:445184/1‑151 (MQ=255)                                                                                             gcgcATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAg                                                      <  1:339917/151‑1 (MQ=255)                                                                                              cgcATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGc                                                     <  1:483805/151‑1 (MQ=255)                                                                                                cATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCAc                                                   <  2:784332/151‑1 (MQ=255)                                                                                                 aTAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCAcc                                                  >  2:72586/1‑151 (MQ=255)                                                                                                 aTAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCAcc                                                  >  1:724107/1‑151 (MQ=255)                                                                                                   aaGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTg                                                <  2:13763/151‑1 (MQ=255)                                                                                                          cGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGcc                                         >  1:486947/1‑151 (MQ=255)                                                                                                            ccATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAg                                       <  1:575617/151‑1 (MQ=255)                                                                                                             cATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGc                                      >  2:336994/1‑151 (MQ=255)                                                                                                              aTTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGc                                      >  1:672270/1‑150 (MQ=255)                                                                                                                     cATTGGTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGAAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCt                               <  1:95435/150‑1 (MQ=255)                                                                                                                          gTGAAAATCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCcgcg                         >  2:455676/1‑151 (MQ=255)                                                                                                                                 tCACAATGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCGGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGc                    >  2:420011/1‑149 (MQ=255)                                                                                                                                     aaTGCCGTCAGCAGAGTGTGGCGGCCCGATAGTGTTGTGGAGCAAAAAAAAGAAGCCGGTAGCCACCACAACCCCCGGCAGGGCGAGAATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGtt                <  2:318537/149‑1 (MQ=255)                                                                                                                                     aaTGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcg              >  2:744738/1‑151 (MQ=255)                                                                                                                                     aaTGCCGCCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcg              >  2:255510/1‑151 (MQ=255)                                                                                                                                      aTGCCGTCAGCAGAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcgc             >  2:814967/1‑151 (MQ=255)                                                                                                                                         ccGTCAGCGGAGGGTGGCAGGCCAATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGa           <  2:447691/150‑1 (MQ=255)                                                                                                                                           gTCAGCAGAGTGTGGCAGGCCGATAGTTTTTTTGAGCAGTAAAAAGAGGCCGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt         <  1:684464/150‑1 (MQ=255)                                                                                                                                             cagcagAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCt      >  2:817940/1‑151 (MQ=255)                                                                                                                                              agcagGGTGGGGCAGGCCGATAGTGTTGTTGAGCAGAAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt         <  1:431079/147‑1 (MQ=255)                                                                                                                                              agcagAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTc     >  2:252909/1‑151 (MQ=255)                                                                                                                                              agcagAGTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTc     >  1:891707/1‑151 (MQ=255)                                                                                                                                                 agagTGGGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTGGCCAGCACAACCCCCGGCATGGCGGGGATCAACATGCCGCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAt  <  1:452499/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAT  >  NC_000913/73169‑73463

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