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

AAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACG  >  NC_000913/72098‑72361                                                                                                                                    |                                                                                                                                     aaaGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  2:620452/139‑1 (MQ=255) aaaGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  2:604442/139‑1 (MQ=255)   aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  1:811412/139‑1 (MQ=255)            ggTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAAt                                                                                                                    <  1:705527/139‑1 (MQ=255)                   gagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCg                                                                                                             <  1:1059988/139‑1 (MQ=255)                      aTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAg                                                                                                          >  2:384467/1‑139 (MQ=255)                       tAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGc                                                                                                         >  1:558768/1‑139 (MQ=255)                               tgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCt                                                                                                 >  1:966736/1‑139 (MQ=255)                                    gTTGCCTGATGTGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            >  1:1222347/1‑139 (MQ=255)                                    gTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            <  2:966736/139‑1 (MQ=255)                                    gTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            >  1:137793/1‑139 (MQ=255)                                       gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCa                                                                                                        <  1:1306635/124‑1 (MQ=255)                                       gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCa                                                                                                        >  2:1306635/1‑124 (MQ=255)                                        ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTca                                                                                        >  2:101739/1‑139 (MQ=255)                                         cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTcaa                                                                                       <  1:1232386/139‑1 (MQ=255)                                                  cGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  1:64777/90‑1 (MQ=255)                                                  cGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  2:64777/1‑90 (MQ=255)                                                         cgcgTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                       <  2:717913/139‑1 (MQ=255)                                                         cgcgTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                       <  1:384467/139‑1 (MQ=255)                                                           cgTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAcc                                                                     <  1:1229550/139‑1 (MQ=255)                                                             tCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  1:487739/81‑1 (MQ=255)                                                             tCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             >  2:487739/1‑81 (MQ=255)                                                                       ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTc                                                                                 >  2:527702/1‑115 (MQ=255)                                                                       ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTc                                                                                 <  1:527702/115‑1 (MQ=255)                                                                          aCAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGt                                                                           >  1:614617/1‑118 (MQ=255)                                                                          aCAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGt                                                                           <  2:614617/118‑1 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTATCCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            <  1:905135/95‑1 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  1:951801/60‑1 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  2:951801/1‑60 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             >  1:1158420/1‑62 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  2:1158420/62‑1 (MQ=255)                                                                                cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                            >  2:905135/1‑95 (MQ=255)                                                                                   aCTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             <  2:441524/139‑1 (MQ=255)                                                                                   aCTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             <  2:1290766/139‑1 (MQ=255)                                                                                         cGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCg                                                                                                             <  1:401883/69‑1 (MQ=255)                                                                                         cGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCg                                                                                                             >  2:401883/1‑69 (MQ=255)                                                                                           tAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             >  2:1183534/1‑51 (MQ=255)                                                                                           tAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  1:1183534/51‑1 (MQ=255)                                                                                               tCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCATGCCAGGCGATGCGCCCGTCGGCGACTACTACCGa                                 >  2:108243/1‑139 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  2:1288329/1‑43 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  1:1288329/43‑1 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  2:1198375/45‑1 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             >  1:1198375/1‑45 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGc                                                                                                         >  2:1315163/1‑65 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGc                                                                                                         <  1:1315163/65‑1 (MQ=255)                                                                                                 ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgc                               <  1:196719/139‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             >  2:433328/1‑124 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             >  1:776186/1‑124 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             <  2:776186/124‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                             <  1:433328/124‑1 (MQ=255)                                                                                                  gATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgcg                              >  1:268279/1‑139 (MQ=255)                                                                                                    tAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCATGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGcgcg                            <  1:108243/139‑1 (MQ=255)                                                                                                        gCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  1:1025044/1‑36 (MQ=21)                                                                                                        gCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  2:1025044/36‑1 (MQ=21)                                                                                                           ttCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTAc                                                                                              >  2:632402/1‑66 (MQ=255)                                                                                                           ttCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTAc                                                                                              <  1:632402/66‑1 (MQ=255)                                                                                                                  cgcATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgcc              >  1:613786/1‑139 (MQ=255)                                                                                                                     aTCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCa           <  2:268279/139‑1 (MQ=255)                                                                                                                              tGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACg  <  2:613786/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     AAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACG  >  NC_000913/72098‑72361

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