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	98.2%	179.0 / ‑5.6	57	*233* (TAG→TAA)	thiQ	thiamine/thiamine pyrophosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/28);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCT  >  NC_000913/72098‑72363                                                                                                                                    |                                                                                                                                        aaaGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                  <  2:223426/139‑1 (MQ=255)        ttGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                  >  1:221928/1‑132 (MQ=255)        ttGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                  <  2:221928/132‑1 (MQ=255)        ttGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑ccc                                                                                                                          >  1:666320/1‑139 (MQ=255)           gggTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCaa                                                                                                                       <  1:470041/139‑1 (MQ=255)               ttaagagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTaa                                                                                                                                      >  1:256321/5‑121 (MQ=255)               ttAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTaa                                                                                                                                      <  2:256321/121‑1 (MQ=255)                 aagagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGc                                                                                                                 >  1:766161/3‑139 (MQ=255)                     gaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGaa                                                                                                             <  2:666320/139‑1 (MQ=255)                        aGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCTGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCa                                                                                                          >  2:1307080/1‑139 (MQ=255)                        aGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCa                                                                                                          >  2:920750/1‑139 (MQ=255)                               tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGCGCATACTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCt                                                                                                   <  1:1130840/139‑1 (MQ=255)                                 cGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAGCCCg                                                                                                                                  >  1:1064087/1‑107 (MQ=255)                                 cGGGTTGCATGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                  <  2:1064087/107‑1 (MQ=255)                                  gggTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTAc                                                                                                <  2:766161/139‑1 (MQ=255)                                    gTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                              >  2:1354084/1‑139 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAAt                                                                                                                      >  1:1406969/1‑111 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAAt                                                                                                                      <  2:1406969/111‑1 (MQ=255)                                                       gCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCa                                                                           >  1:555945/1‑139 (MQ=255)                                                         cgcgTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                         <  2:555945/139‑1 (MQ=255)                                                                 aTCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑c                                                                                                                            >  2:1679427/1‑80 (MQ=255)                                                                 aTCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑c                                                                                                                            <  1:1679427/80‑1 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCa                                                                                                          <  2:641731/83‑1 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCa                                                                                                          >  1:641731/1‑83 (MQ=255)                                                                                 gCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTACT‑CCCCAATAGTGCCGAAGCACTCGCCTTAACGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGTGCCCGTCg                                                 <  1:920750/139‑1 (MQ=255)                                                                                 gCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                                 <  2:840452/139‑1 (MQ=255)                                                                                   aCTACCCGTAGGTCGGATAAGGCGTTCACGCTGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                               <  1:1307080/139‑1 (MQ=255)                                                                                      aCCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCt                                                                                            >  1:558288/1‑91 (MQ=255)                                                                                      aCCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCt                                                                                            <  2:558288/91‑1 (MQ=255)                                                                                        ccGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                  <  2:1076565/52‑1 (MQ=255)                                                                                        ccGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                                  >  1:1076565/1‑52 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTg                                                                                                                  >  1:495970/1‑58 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTg                                                                                                                  <  2:495970/58‑1 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgc                                 >  1:830050/1‑139 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGcc                                                                                                                >  2:124212/1‑59 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGcc                                                                                                                <  1:124212/59‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCa                                                                                                          <  1:685898/65‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCa                                                                                                          >  2:685898/1‑65 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTcaac                                                                                        >  2:1047204/1‑83 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTcaac                                                                                        <  1:1047204/83‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                  <  1:101276/89‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCa                                                                                  >  2:101276/1‑89 (MQ=255)                                                                                                        gCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGc                          >  1:820703/1‑139 (MQ=255)                                                                                                            tCACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATc                      <  2:830050/139‑1 (MQ=255)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCCAATAGTGCCGAAGCACTCGCCt                                                                                                   >  1:1314829/1‑62 (MQ=38)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCCAATAGTGCCGAAGCACTCGCCt                                                                                                   <  2:1314829/62‑1 (MQ=38)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAt                                                                                 >  1:1070074/1‑79 (MQ=255)                                                                                                             cACGCCGCATCCGACAGTGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAt                                                                                 <  2:1070074/79‑1 (MQ=255)                                                                                                                              tGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATc                      >  1:1609936/1‑121 (MQ=255)                                                                                                                              tGCATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATc                      <  2:1609936/121‑1 (MQ=255)                                                                                                                                cATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                         <  2:448187/68‑1 (MQ=255)                                                                                                                                cATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                         >  1:448187/1‑68 (MQ=255)                                                                                                                                cATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTGCCGAGCGCGTGGCGATCCGcgccgc              >  2:1079834/1‑127 (MQ=255)                                                                                                                                cATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCcgc                   <  2:246470/122‑1 (MQ=255)                                                                                                                                cATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCcgc                   >  1:246470/1‑122 (MQ=255)                                                                                                                                cATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgccgc              <  1:1079834/127‑1 (MQ=255)                                                                                                                                cATATTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCt  >  1:1264461/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                        AAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAAT‑CCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCT  >  NC_000913/72098‑72363

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