Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,246,924	T→C	G151G (GGA→GGG)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,246,924	0	T	C	100.0%	220.1 / NA	61	G151G (GGA→GGG)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/27);  total (34/27)

GGAATATGACCGTAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATAT  >  NC_000913/1246793‑1247059                                                                                                                                    |                                                                                                                                        ggAATATGACCGTAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCgg                                                                                                                                 <  2:1142008/140‑1 (MQ=255)   aaTATGACCGTAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCa                                                                                                                               <  1:1062663/140‑1 (MQ=255)   aaTATGACCGTAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCa                                                                                                                               <  2:88012/140‑1 (MQ=255)   aaTATGACCGTAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCa                                                                                                                               <  1:907070/140‑1 (MQ=255)             tAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAgg                                                                                                                     >  2:1178236/1‑140 (MQ=255)             tAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAgg                                                                                                                     >  2:303081/1‑140 (MQ=255)               aGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGtt                                                                                                                   >  1:737758/1‑140 (MQ=255)                 tgtCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTcc                                                                                                                 >  2:1205843/1‑140 (MQ=255)                   tCTATTTCATGAGTAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCgg                                                                                                               >  1:691984/1‑140 (MQ=255)                    cTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGc                                                                                                              <  2:167555/140‑1 (MQ=255)                             gAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAAc                                                                                                     >  2:576226/1‑140 (MQ=255)                                        gCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGaa                                                                                               <  2:943042/135‑1 (MQ=255)                                        gCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGaa                                                                                               >  1:943042/1‑135 (MQ=255)                                                     cgACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGgt                                                                                 <  2:1230509/136‑1 (MQ=255)                                                     cgACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGgt                                                                                 >  1:1230509/1‑136 (MQ=255)                                                           tATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTa                                                                       >  2:667568/1‑140 (MQ=255)                                                             tCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGa                                                                     <  1:721909/140‑1 (MQ=255)                                                              cccAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGaa                                                                    <  2:1164044/140‑1 (MQ=255)                                                              cccAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGaa                                                                    <  2:608286/140‑1 (MQ=255)                                                                  gTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGaa                                                                    >  2:594027/1‑136 (MQ=255)                                                                  gTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGaa                                                                    <  1:594027/136‑1 (MQ=255)                                                                        gCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGaa                                                                    >  1:674012/1‑130 (MQ=255)                                                                        gCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGaa                                                                    <  2:674012/130‑1 (MQ=255)                                                                                   gTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGt                                               <  1:980273/140‑1 (MQ=255)                                                                                          cATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCAGTCaa                                        >  1:840721/1‑140 (MQ=255)                                                                                              gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATg                                    >  1:1104475/1‑140 (MQ=255)                                                                                              gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATg                                    >  2:347608/1‑140 (MQ=255)                                                                                              gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATg                                    >  1:1200723/1‑140 (MQ=255)                                                                                              gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATg                                    >  1:13260/1‑140 (MQ=255)                                                                                              gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATg                                    >  1:237507/1‑140 (MQ=255)                                                                                                  aGTAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTcg                                <  2:691984/140‑1 (MQ=255)                                                                                                    tAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTcgcg                              <  1:303081/140‑1 (MQ=255)                                                                                                    tAACTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAACCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCTAAAGTCCGTCAATATGTTcgcg                              <  1:1178236/140‑1 (MQ=255)                                                                                                      aCTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCa                            >  1:541468/1‑140 (MQ=255)                                                                                                      aCTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCa                            <  2:1200723/140‑1 (MQ=255)                                                                                                      aCTGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCAGTCAATATGTTCGCGCa                            <  2:840721/140‑1 (MQ=255)                                                                                                        tGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGt                          >  2:714723/1‑140 (MQ=255)                                                                                                        tGTCCCAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGt                          <  2:541468/140‑1 (MQ=255)                                                                                                           cccAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGAc                       <  1:714723/140‑1 (MQ=255)                                                                                                           cccAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGAc                       >  1:1120981/1‑140 (MQ=255)                                                                                                           cccAGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGAc                       >  1:1122680/1‑140 (MQ=255)                                                                                                              aGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTcgcg                              <  1:1205646/130‑1 (MQ=255)                                                                                                              aGTAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTcgcg                              >  2:1205646/1‑130 (MQ=255)                                                                                                                tAATATACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGccc                  >  1:506082/1‑140 (MQ=255)                                                                                                                  atatACCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTTGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCct                >  1:349576/1‑140 (MQ=255)                                                                                                                   tataCCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCctc               <  2:1104475/140‑1 (MQ=255)                                                                                                                    ataCCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGt                          <  2:440026/128‑1 (MQ=255)                                                                                                                    ataCCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGt                          >  1:440026/1‑128 (MQ=255)                                                                                                                    ataCCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCctct              >  1:806885/1‑140 (MQ=255)                                                                                                                    ataCCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCctct              >  2:947345/1‑140 (MQ=255)                                                                                                                      aCCTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTgg            <  2:13260/140‑1 (MQ=255)                                                                                                                        cTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCg          >  1:386539/1‑140 (MQ=255)                                                                                                                        cTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCg          >  2:172986/1‑140 (MQ=255)                                                                                                                        cTCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCg          >  1:493259/1‑140 (MQ=255)                                                                                                                         tCGCGAAAGTGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCgg         >  2:403588/1‑140 (MQ=255)                                                                                                                         tCGCGAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCgg         <  1:576226/140‑1 (MQ=255)                                                                                                                          cgcgAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGa        >  1:158451/1‑140 (MQ=255)                                                                                                                          cgcgAAAGCGCCCGCCCGGCACGACATAAGGATCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGa        >  2:1116145/1‑140 (MQ=255)                                                                                                                           gcgAAAGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGaa       >  2:885546/1‑140 (MQ=255)                                                                                                                             gAAAGTGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACa     <  1:403588/140‑1 (MQ=255)                                                                                                                                aGCGCCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACatat  <  2:1122680/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        GGAATATGACCGTAAGTGTCTATTTCATGAGCAAAATTGGCCACCATATCCGCGACTTTATCCCAGTGACCGCTTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATAT  >  NC_000913/1246793‑1247059

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