Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,790,550	G→A	P266S (CCG→TCG)	tdh ←	threonine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,790,550	0	G	A	100.0%	197.3 / NA	56	P266S (CCG→TCG)	tdh	threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/24);  total (32/24)

GATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTCCGCCATCACGTC  >  NC_000913/3790413‑3790682                                                                                                                                          |                                                                                                                                     gATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGacgacg                                                                                                                                    <  2:308631/140‑1 (MQ=255)     ggCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGacgacg                                                                                                                                    >  1:319061/1‑136 (MQ=255)     ggCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGAcggaatacc                                                                                                                                >  1:973443/1‑134 (MQ=255)     ggCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAAt                                                                                                                                >  1:30556/1‑140 (MQ=255)     ggCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAAt                                                                                                                                >  2:691821/1‑140 (MQ=255)     ggCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAAt                                                                                                                                >  1:447467/1‑140 (MQ=255)      gCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATa                                                                                                                               <  2:447998/140‑1 (MQ=255)        gAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATAcc                                                                                                                             <  2:96494/140‑1 (MQ=255)              tCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCAt                                                                                                                       <  1:930328/140‑1 (MQ=255)                 aGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGc                                                                                                                    >  2:773506/1‑140 (MQ=255)                   gCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCaa                                                                                                                  <  1:691821/140‑1 (MQ=255)                        aCTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACgg                                                                                                             >  2:346353/1‑140 (MQ=255)                                    ccgccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACgg                                                                                                             <  2:321827/128‑1 (MQ=255)                                    ccgccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACgg                                                                                                             >  1:321827/1‑128 (MQ=255)                                     cgccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCa                                                                                                <  2:968914/140‑1 (MQ=255)                                      gccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGt                                                                                                      >  1:581139/1‑133 (MQ=255)                                      gccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGt                                                                                                      <  2:581139/133‑1 (MQ=255)                                       ccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATg                                                                                              <  2:344471/140‑1 (MQ=255)                                           cTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGgtgt                                                                                          >  2:698903/1‑140 (MQ=255)                                                 ccAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCa                                                                                    <  1:698903/140‑1 (MQ=255)                                                         cAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGa                                                                            >  2:788773/1‑140 (MQ=255)                                                              aTCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGc                                                                       <  1:788773/140‑1 (MQ=255)                                                               tctcACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCt                                                                      >  2:202476/1‑140 (MQ=255)                                                                      aCCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGgcggcg                                                               >  2:230144/1‑140 (MQ=255)                                                                         gTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCAc                                                            <  2:30556/140‑1 (MQ=255)                                                                          tAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCAcc                                                           >  2:354678/1‑140 (MQ=255)                                                                          tAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCAcc                                                           >  2:147358/1‑140 (MQ=255)                                                                             aTACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGa                                                        <  1:364863/140‑1 (MQ=255)                                                                               aCCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACa                                                      >  1:111088/1‑140 (MQ=255)                                                                                  tttAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGa                                                        <  1:481392/135‑1 (MQ=255)                                                                                  tttAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGa                                                        >  2:481392/1‑135 (MQ=255)                                                                                  tttAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACAttt                                                   >  2:335258/1‑140 (MQ=255)                                                                                    tAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTcc                                                 >  1:913003/1‑140 (MQ=255)                                                                                         aaCAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGAcc                                            >  1:772700/1‑140 (MQ=255)                                                                                         aaCAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGAcc                                            <  2:447467/140‑1 (MQ=255)                                                                                           cAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGa                                          >  2:471089/1‑140 (MQ=255)                                                                                              gCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACAt                                       >  1:87938/1‑140 (MQ=255)                                                                                                   tAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCaaaa                                  >  1:236989/1‑140 (MQ=255)                                                                                                   tAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCaaaa                                  >  2:786029/1‑140 (MQ=255)                                                                                                       gATCACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCtt                              <  1:230144/140‑1 (MQ=255)                                                                                                          cACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCgg                           <  2:972773/140‑1 (MQ=255)                                                                                                          cACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCgg                           <  2:913003/140‑1 (MQ=255)                                                                                                          cACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCgg                           <  2:87938/140‑1 (MQ=255)                                                                                                          cACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCgg                           >  1:686238/1‑140 (MQ=255)                                                                                                          cACTTTGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCgg                           <  1:621915/140‑1 (MQ=255)                                                                                                             tttGGTCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCa                        >  1:523203/1‑140 (MQ=255)                                                                                                                  tCCAGTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCt                   <  1:513402/140‑1 (MQ=255)                                                                                                                      gTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACt               >  1:588803/1‑140 (MQ=255)                                                                                                                      gTCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACt               <  2:772700/140‑1 (MQ=255)                                                                                                                       tCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGTCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTc                         <  2:388128/129‑1 (MQ=255)                                                                                                                       tCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGTCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTc                         >  1:388128/1‑129 (MQ=255)                                                                                                                       tCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTc                         <  2:744398/129‑1 (MQ=255)                                                                                                                       tCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTc                         >  1:744398/1‑129 (MQ=255)                                                                                                                       tCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAAc                <  1:244173/138‑1 (MQ=255)                                                                                                                       tCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAAc                >  2:244173/1‑138 (MQ=255)                                                                                                                       tCGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTc              >  1:712940/1‑140 (MQ=255)                                                                                                                        cGATAGACATATCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTcc             <  1:354678/140‑1 (MQ=255)                                                                                                                                 tatCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTCCGCCATCACg    >  1:50975/1‑140 (MQ=255)                                                                                                                                 tatCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTCCGCCATCACg    >  2:454305/1‑140 (MQ=255)                                                                                                                                   tCAGACGACGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTCCGCCATCACGTc  >  1:57431/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                     GATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTCCGCCATCACGTC  >  NC_000913/3790413‑3790682

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