Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,465,733	A→G	Y150H (TAC→CAC)	treB ←	trehalose‑specific PTS enzyme IIBC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,465,733	0	A	G	100.0%	128.7 / NA	39	Y150H (TAC→CAC)	treB	trehalose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/11);  total (28/11)

CGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCT  >  NC_000913/4465598‑4465867                                                                                                                                        |                                                                                                                                        cGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCAtt                                                                                                                                     >  1:1131482/1‑140 (MQ=255) cGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCAtt                                                                                                                                     >  1:509205/1‑140 (MQ=255)   ccAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTg                                                                                                                                   >  2:10554/1‑140 (MQ=255)      aGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTgcgc                                                                                                                                <  1:10554/140‑1 (MQ=255)       gCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGcc                                                                                                                               <  2:1465710/140‑1 (MQ=255)       gCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGcc                                                                                                                               <  1:156517/140‑1 (MQ=255)            aTACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGt                                                                                                                          <  2:996269/140‑1 (MQ=255)            aTACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGt                                                                                                                          >  2:329315/1‑140 (MQ=255)             tACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTc                                                                                                                         >  2:1305707/1‑140 (MQ=255)                 aaGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGAc                                                                                                                     <  2:88905/140‑1 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAg                                                                                                                             <  2:361570/126‑1 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAg                                                                                                                             >  1:361570/1‑126 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGc                                                                                                               >  1:53728/1‑140 (MQ=255)                               cgccCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGc                                                                                                       <  2:932974/140‑1 (MQ=255)                                gccCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCa                                                                                                      <  1:329315/140‑1 (MQ=255)                                          ttCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGAt                                                                                           <  2:509205/140‑1 (MQ=255)                                          ttCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGAt                                                                                           <  2:1131482/140‑1 (MQ=255)                                             aCCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATcac                                                                                        >  1:1158547/1‑140 (MQ=255)                                              ccGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATcaca                                                                                       <  2:473585/140‑1 (MQ=255)                                               cGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACAt                                                                                      >  1:851175/1‑140 (MQ=255)                                                                  accgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCaa                                                                   >  1:219987/1‑140 (MQ=255)                                                                   ccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCaaa                                                                  >  1:1242212/1‑140 (MQ=255)                                                                    gggCAGGTAGAAGAAGATCGTTTCACCGATCAACCACAAAAAATCGTAGATCGTTTGCAGGGAGGGGTGCTTTTGCGCCAGCGTCTGACCTTTGCTCATGGGAAAAT‑CGCCGATCACATTGCAAAACCCGAGAATAAAAc                                                                 <  2:219987/139‑1 (MQ=255)                                                                             gaagaagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTaa                                                        >  1:293458/1‑140 (MQ=255)                                                                                  agaTCGCTTCCCCAATCACCCAAAAAAAATCGTAGATCTTTTGCAGGAAAGGGTGCTTTGGCGCCAGCGTCTGCCCTTTGCTCAGGGGAAAAT‑CGCCGATCACATTGAAAAACCCGAGAAAAAAACCGCCGCTAATCAAc                                                   <  2:888804/140‑1 (MQ=255)                                                                                    aTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATcctgtc                                                 >  1:1399844/1‑135 (MQ=255)                                                                                    aTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTACAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAAcgc                                                 >  1:1365686/1‑140 (MQ=255)                                                                                      cGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCgg                                               >  1:1204152/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCATCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                             >  1:257146/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                             >  2:1076071/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                             >  2:313192/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCAAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                             <  2:53728/140‑1 (MQ=255)                                                                                         ttCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgca                                            >  1:922796/1‑140 (MQ=255)                                                                                           cACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgcagc                                          >  1:1245492/1‑140 (MQ=255)                                                                                              cGACCACCCACAGAAAATCGTAGTTCGTTTGGGGGGAAGGGTGCTTTTGCGCCGGCTTCTGACCGTTGCTCATGGGAAAAT‑CGCCGATCACATTGCAAAACCCGAGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAAc                                       <  2:1245492/140‑1 (MQ=255)                                                                                                aTCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACgg                                     >  1:35793/1‑140 (MQ=255)                                                                                                          aaaaaTCGTAGATCGTTTGAAGGGAGGGGTGCATTTGCCCCAGCGTCTGACCGTTGCTCAGGGAAAAAC‑CGCCGATCACATTGCAAAACCCGAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATc                           <  2:293458/139‑1 (MQ=255)                                                                                                              atggTAGATCGTTTGGAGGGAGGGGTGTTTTTGCGCCGGCGTCTGCCTGTTGCTCAGGGAAAAAC‑CGCCGATCACATTGCAAAACCCGAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcg                       <  2:992241/137‑1 (MQ=255)                                                                                                                cGTAGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcgcg                     >  1:149939/1‑140 (MQ=255)                                                                                                                 agaGATGGTTTGCAGGGAAGGGTGCTTTCGCGCCGGCGTCTGCCCGTTGCTCAGGGGAAAAC‑CGCCGATCAAATTGAAAAAACCGAGAAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGa                    <  2:70156/138‑1 (MQ=255)                                                                                                                   aGATCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGaaa                  >  1:393275/1‑140 (MQ=255)                                                                                                                    gATCGTTTGCAGGGAGGGGTGCTTTTGCGCCAGCGTCTGCCCGTTGCTCAGGGGAAAAC‑CGCCGATCAAATGGCAAAACCCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAAt                 <  2:875726/140‑1 (MQ=255)                                                                                                                     aTCGTTTGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATg                >  1:1166351/1‑140 (MQ=255)                                                                                                                         tttGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag            >  2:947241/1‑140 (MQ=255)                                                                                                                         tttGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag            >  1:992241/1‑140 (MQ=255)                                                                                                                         ttggCAGGGAAGGGTGCTTTTGCGCCAGCGTCTGACCGTTGCTCAGGGGAAAAC‑CGCCGATCACATTGCAAAAACCGAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag            <  2:1166351/137‑1 (MQ=255)                                                                                                                          ttGCAGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagaga           >  2:551890/1‑140 (MQ=255)                                                                                                                              ggggAGGGGTGCTTTGGCGCCGGCGTCTGCCCGTTGCTCAGGGGCAAAC‑CGCCGATCACATTGCAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCaa       <  2:539984/139‑1 (MQ=255)                                                                                                                              aGGGAAGGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATcct       >  1:70156/1‑138 (MQ=255)                                                                                                                                   ggggTGCTTTGGCCCCGGCGTCGGCCCGTTCCTCA‑GGGAAAATCCCCCGACCAAATGGAAAAAACGGAGAACAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCt  <  2:331782/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        CGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAAT‑CGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCT  >  NC_000913/4465598‑4465867

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: ‑