Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,387	C→T	Q791* (CAG→TAG)	aceE →	pyruvate dehydrogenase E1 component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	125,387	0	C	T	100.0%	124.1 / NA	49	Q791* (CAG→TAG)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/24);  total (25/24)

CCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG  >  NC_000913/125251‑125524                                                                                                                                         |                                                                                                                                          cgtGGGGAAAGATGGCGTCGTGGGAACTGAGGTGGGTAGCGTGACCTCCTGCACCGAGCTGGCGCGGGAGGGGCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                     <  2:9303/141‑1 (MQ=255)   tGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGGACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                                                                   >  1:698540/1‑143 (MQ=255)   tGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTGTATCGCTTAGGtgatg                                                                                                                                   >  1:103029/1‑143 (MQ=255)   tGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                                                                   >  1:316702/1‑143 (MQ=255)        aaaGGTCACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                                                              <  2:340708/143‑1 (MQ=255)         aaGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGTATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                                                             <  1:437523/143‑1 (MQ=255)                  gCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCa                                                                                                                    <  2:224801/143‑1 (MQ=255)                       ggTTCTGACGTTTATAGCGTGCCCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGGGACCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGc                                                                                                               <  2:517622/143‑1 (MQ=255)                        gTTCTGACGTTTATCGCGTGCCCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAACCTCCGCCAGTCCCGTATCTCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCa                                                                                                              <  1:793603/143‑1 (MQ=255)                                gTTCATAGCGTGACCTCCTTCACCGCGCTGGCGCGTGAGGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGa                                                                                                      <  1:120959/143‑1 (MQ=255)                                    ataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGGCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTAGATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCAGCTACCGCCtat                                                                                                  >  2:112851/1‑143 (MQ=255)                                               tcccTCAACGAGCTCGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGtt                                                                                       <  2:633885/143‑1 (MQ=255)                                                cctTTACCGGGCTTGCCCGGGATGGTCAGGATTGTGAACTCTGGAACATGCTGCGCCCGCTGGAAACTCCGAGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                      <  2:428623/143‑1 (MQ=255)                                                ccgtCACCGAGCCGGCGCGGGATGGGCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                      <  1:184047/140‑1 (MQ=255)                                                     aCCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATCTGAACCTGTTCGCTGa                                                                                 >  1:367668/1‑143 (MQ=255)                                                       cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGCAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                               >  1:12884/1‑143 (MQ=255)                                                         aGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAg                                                                             <  2:551115/143‑1 (MQ=255)                                                          gCTGGCGCGTGATGGTCAGGTTTGGGCACGCTGGATCATGGTGCGCCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAgg                                                                            <  2:675194/143‑1 (MQ=255)                                                          gCTGGCGCGTGATGGTCAGGATTGTGAAAGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCAGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCCTATCTGAAAATGTTCGCTAAGAAgg                                                                            >  2:1004762/1‑143 (MQ=255)                                                           cTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                                                                           <  1:417205/143‑1 (MQ=255)                                                             ggCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCCCCCGCCGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTcc                                                                         >  1:407549/1‑143 (MQ=255)                                                              gcgcGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTACGTATATCGCGTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCCCTGAGCAGGTCCg                                                                        >  2:622877/1‑143 (MQ=255)                                                                gcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGCTCCGTATATCGCTTCGGTCATGAACGACGCTCCGGCAGTGGCCTCTACCGACTATATCAAACTATTCTCGAAGCATGTCCGTa                                                                      >  2:257737/1‑143 (MQ=255)                                                                   tGATGGTCAGGATTGTGACCGCTGGAACATGCTGCACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACtt                                                                   <  2:235273/143‑1 (MQ=255)                                                                   tGATGGTCAGGATTGTGAACGCTGGAATATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACtt                                                                   >  1:211701/1‑143 (MQ=255)                                                                    gATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTAGATGAAACTGTTCGCTGAGCAGGTCCGTAcctc                                                                  >  2:112640/1‑140 (MQ=255)                                                                             gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGCAAGCCGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGc                                                         >  1:904002/1‑143 (MQ=255)                                                                                  tgACCGCTGGCACATGCTCCACCCGCTGGCACCTCCGCGCGTTCCGTATATCGCTTAGGTCATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacg                                                    <  1:269047/143‑1 (MQ=255)                                                                                     aCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATAAACGACTCTCCGGCAGTGGCATCTACCGACTATATGAAGCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACt                                                 >  1:152902/1‑143 (MQ=255)                                                                                                tgctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCCGGTCCGTACTTACGTACCGGCTGACCACTACCGCGtccgg                                      >  2:832821/1‑139 (MQ=255)                                                                                                 gcggcACGCGCTGTCGACTCCGCGCGTTCCGTATATCGCCTGGGTGAGGGACGACGCTCCGGCAGTGGCATCTACCGACTATATTAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                     <  2:12950/140‑1 (MQ=255)                                                                                                  ctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTggg                                    >  1:573966/1‑143 (MQ=255)                                                                                                    gcACCCGCTGGAAACTCCGCGCGTACCGTATATCGCTTAGGTGATGAACGACGCGCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTa                                  <  1:415191/143‑1 (MQ=255)                                                                                                       cccGCTGGAGCCTCCGCGCCTGCCGTCTATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGGCGACTACCGCGTACTGGGTACTg                               <  2:765241/143‑1 (MQ=255)                                                                                                        ccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGCGATGAACGACGCCCCGGCAGTGGCATCTACCGACTATATGAAACTATTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGCCGACTACCGCGTACGGGGTCCTGa                              >  2:433510/1‑143 (MQ=255)                                                                                                             ggAAACTCCGCGCGGTCCGAATATCGCGTAGGTGATGAACGAGGCTCCGGGAGTGGGATCTGGCGAATATATCGGACTGCTAGCTGAGCAGGACCGTACTTACGTACCGGCTGACGACTACCGCGTCCTGTGTACTGATGTCt                         >  2:492574/1‑143 (MQ=255)                                                                                                              gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCtt                        >  2:485326/1‑143 (MQ=255)                                                                                                                aaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCCCGTGCTGGGTACTCATGGCGTCg                      >  2:283210/1‑143 (MQ=255)                                                                                                                   tcggcgcgcTCTGTACCTCGCTTCGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCGTGGTCGTTCCGCCCCACTGCGCAGTTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTc                   <  1:248441/140‑1 (MQ=255)                                                                                                                   tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGGGCAGGTCCGTACTTCCGTACCGGCTGGCGACTACCGCGTACTGGGTACTGATGGCTTCGGTc                   >  1:73822/1‑143 (MQ=255)                                                                                                                      gcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGtt                <  1:1031341/143‑1 (MQ=255)                                                                                                                      gcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGtt                >  1:372076/1‑143 (MQ=255)                                                                                                                       cgcgTGCCGTAGAGCGCTTAGGTGATGGACGACGCCCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTc               <  2:228385/143‑1 (MQ=255)                                                                                                                        gggtgccGGATATCGCTTAGGTGTTGGAGGGCGCGCCGGCAGGGGCATGTGCCGACGGTATGAAGCTGTTCGCTGAGCGGGTCCGTACTTACGTACCGGCTGGCGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTcc              <  2:344150/138‑1 (MQ=255)                                                                                                                           ttCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGAc           >  2:922317/1‑143 (MQ=255)                                                                                                                                tattTCGCGTAGGTGTTGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCg      <  2:407549/143‑1 (MQ=255)                                                                                                                                tataTCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGCCTTCGGTCGCTCCGACAGCCg      >  2:304297/1‑143 (MQ=255)                                                                                                                                   accGCTTAGGCGCTGAACGACGCTCCGGGAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTga   <  1:543679/141‑1 (MQ=255)                                                                                                                                    tCGCTTCGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTGCTGATGGCTCCGGTCGTTCCGACAGCCGTgag  >  1:981998/1‑143 (MQ=255)                                                                                                                                         |                                                                                                                                          CCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG  >  NC_000913/125251‑125524

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