Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,750,344	T→C	Y87H (TAT→CAT)	sgbE →	L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,750,344	0	T	C	100.0%	206.2 / NA	57	Y87H (TAT→CAT)	sgbE	L‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/31);  total (26/31)

TGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCTATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGACCGCAGAGG  >  NC_000913/3750197‑3750491                                                                                                                                                    |                                                                                                                                                    tggatggTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATg                                                                                                                                                  >  1:339045/1‑151 (MQ=255)      ggTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGa                                                                                                                                              <  1:486043/150‑1 (MQ=255)        tAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGAGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAAt                                                                                                                                           >  2:899408/1‑151 (MQ=255)            cAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGc                                                                                                                                       >  1:336790/1‑151 (MQ=255)             aaaCCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCg                                                                                                                                      >  2:587953/1‑151 (MQ=255)             aaaCCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCg                                                                                                                                      <  1:302793/151‑1 (MQ=255)               aCCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCgg                                                                                                                                     >  2:648426/1‑150 (MQ=255)                   tCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTAtt                                                                                                                                 <  1:1054827/150‑1 (MQ=255)                   tCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTAtt                                                                                                                                 <  2:432448/150‑1 (MQ=255)                   tCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTAt                                                                                                                                  <  1:386014/149‑1 (MQ=255)                    ccGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATtg                                                                                                                                <  2:597004/150‑1 (MQ=255)                     cGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATtgtg                                                                                                                              <  1:934528/151‑1 (MQ=255)                     cGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATtgtg                                                                                                                              >  1:851064/1‑151 (MQ=255)                      ggCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGc                                                                                                                             <  1:860426/151‑1 (MQ=255)                        cGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGACAGCAGAAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGcct                                                                                                                           >  2:1046255/1‑149 (MQ=255)                            gAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATAccc                                                                                                                       >  1:865514/1‑151 (MQ=255)                              gTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCAc                                                                                                                     <  1:727761/151‑1 (MQ=255)                               tACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACt                                                                                                                    >  1:693580/1‑151 (MQ=255)                                acgacgTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTc                                                                                                                   <  2:970467/151‑1 (MQ=255)                                 cgacgTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTc                                                                                                                   >  2:408745/1‑150 (MQ=255)                                  gacgTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTcgc                                                                                                                 >  1:1067183/1‑151 (MQ=255)                                  gacgTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGGGCATACCCACTc                                                                                                                   >  2:380524/1‑149 (MQ=255)                                   acgTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTcgcg                                                                                                                <  1:1065470/151‑1 (MQ=255)                                     gTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCgcgc                                                                                                               <  2:854384/150‑1 (MQ=255)                                      tgatgaCCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGcgcc                                                                                                              >  1:901195/1‑150 (MQ=255)                                                cgacgaTATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATc                                                                                                   >  1:846348/1‑151 (MQ=255)                                                   cgaTATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTgg                                                                                                >  1:731554/1‑151 (MQ=255)                                                       atGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTca                                                                                             <  1:781367/150‑1 (MQ=255)                                                        tggtggtggTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTcac                                                                                            >  1:562354/1‑150 (MQ=255)                                                              tggtTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCgg                                                                                     >  2:698354/1‑151 (MQ=255)                                                                          ccAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTccc                                                                         <  1:389125/151‑1 (MQ=255)                                                                            aGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTcccc                                                                        >  2:469325/1‑150 (MQ=255)                                                                             gCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGc                                                                      <  2:392880/151‑1 (MQ=255)                                                                                  aaGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTggg                                                                 <  2:672043/151‑1 (MQ=255)                                                                                          ggAAGGCAGCAAAAAACCCTCTTCCGATCCACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCaccacc                                                         <  1:807343/151‑1 (MQ=255)                                                                                          ggAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCaccacc                                                         >  1:57469/1‑151 (MQ=255)                                                                                             aGGCAGCAAAAAACCCTCTTCCGATACCCCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCAccacccac                                                      <  1:388701/151‑1 (MQ=255)                                                                                              ggCAGCAAAAACCCCTCTTCCTATACACCAACGCATCTGGCGCTCTCCCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCCCGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCAccacccac                                                      <  2:959407/150‑1 (MQ=255)                                                                                              ggCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACg                                                     >  1:325604/1‑151 (MQ=255)                                                                                              ggCAGCAAAAAACCCTCTCCCGATACACCAACGCATCTGGCGCTCTACCTTCCCCATGCCGAATTTGGCGGTATTTTGCATCCCCACCCGCGCCACGCCACCATCTGGTCCCAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACg                                                     <  2:465970/151‑1 (MQ=255)                                                                                                    aaaaaaCCCTCTTCCGATACCCCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATCCCCCCTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGatt                                               <  2:596286/151‑1 (MQ=255)                                                                                                      aaaaCCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGattat                                             <  1:268433/151‑1 (MQ=255)                                                                                                        aaCCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTAttt                                           <  2:47022/151‑1 (MQ=255)                                                                                                            ctctTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTAc                                       >  1:384798/1‑151 (MQ=255)                                                                                                                 ccGATACACCAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGc                                  <  2:851016/151‑1 (MQ=255)                                                                                                                         ccAACGCATCTGGCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATcccc                           >  1:366921/1‑150 (MQ=255)                                                                                                                          cAACGCATCTGGCGCTCTCCCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACgg                                     <  2:216010/139‑1 (MQ=255)                                                                                                                                    ggCGCTCTACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAg               <  2:336765/151‑1 (MQ=255)                                                                                                                                          ctACCGTCGCCATGCCGAAATTGGCGGTTTTGTGCATACCCCCTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATg            <  2:214042/148‑1 (MQ=255)                                                                                                                                           tACCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGAc          >  1:448827/1‑149 (MQ=255)                                                                                                                                            aCCGTCGCCATGCCGAAATTGGCGGTATTGTGCATCCCCCCTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCCCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGAc          <  2:1022999/148‑1 (MQ=255)                                                                                                                                            aCCGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGACCg        <  1:24567/150‑1 (MQ=255)                                                                                                                                            aCCATCGCCAGGCCAAAGTTGGGGGTTTTTTCCATACCCACCCGCGCCACGCAGCCATCTGGCCGCAGGCCGGGTTGGATGGCCCCGCCTGGGGCCCCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGACCg        <  2:1011181/150‑1 (MQ=255)                                                                                                                                             ccGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCCCGCGACAGATGAc          <  2:921126/147‑1 (MQ=255)                                                                                                                                             ccGTCGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGACCGc       >  1:299926/1‑150 (MQ=255)                                                                                                                                                 cGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGCCAGATGACCGCAgagg  >  2:287417/1‑151 (MQ=255)                                                                                                                                                 cGCCATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGACCGCAgagg  >  1:948529/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACACCAACGCATCTGGCGCTCTACCGTCGCTATGCCGAAATTGGCGGTATTGTGCATACCCACTCGCGCCACGCCACCATCTGGTCACAGGCCGGGCTGGATCTCCCCGCCTGGGGCACCACCCACGCCGATTATTTTTACGGTGCCATCCCCTGCACGCGACAGATGACCGCAGAGG  >  NC_000913/3750197‑3750491

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