Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,762,741	A→T	K662I (AAA→ATA)	spoT →	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,762,741	0	A	T	100.0%	236.2 / NA	60	K662I (AAA→ATA)	spoT	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (34/26);  total (34/26)

GGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762676‑3762813                                                                  |                                                                         ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTctc                                                                >  1:1516107/1‑74 (MQ=255) ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTctc                                                                >  1:773067/1‑74 (MQ=255) ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTctc                                                                >  1:1381161/1‑74 (MQ=255) ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                                                                >  2:150422/1‑76 (MQ=255) ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                                                                >  1:1442754/1‑76 (MQ=255)    aaaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                                                             >  1:837578/1‑76 (MQ=255)     aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                                                            >  1:1246836/1‑76 (MQ=255)     aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                                                            <  2:775778/76‑1 (MQ=255)      aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                                                           <  1:1213399/76‑1 (MQ=255)        cTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                                                         <  1:1084688/76‑1 (MQ=255)         tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg                                                        <  1:1662891/76‑1 (MQ=255)           aCGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg                                                      >  2:878337/1‑76 (MQ=255)            cggcggCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc                                                     >  1:639986/1‑76 (MQ=255)                  cAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                               <  1:44791/76‑1 (MQ=255)                   aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt                                              >  2:676772/1‑76 (MQ=255)                      tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                           >  1:556302/1‑76 (MQ=255)                      tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                           >  1:2298/1‑76 (MQ=255)                      tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                           >  2:545748/1‑76 (MQ=255)                         caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg                                        <  2:583167/76‑1 (MQ=255)                          accacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa                                       <  1:550466/76‑1 (MQ=255)                           ccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAc                                      >  2:1190924/1‑76 (MQ=255)                           ccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAc                                      >  1:1163783/1‑76 (MQ=255)                                aCTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTc                                 >  1:882029/1‑76 (MQ=255)                                  ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGt                               >  1:1621078/1‑76 (MQ=255)                                  ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACAGCTCGt                               >  1:1621095/1‑76 (MQ=255)                                   tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              >  1:22999/1‑76 (MQ=255)                                   tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              <  1:539319/76‑1 (MQ=255)                                   tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              >  1:165419/1‑76 (MQ=255)                                   tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              >  1:473608/1‑76 (MQ=255)                                   tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              >  1:1671123/1‑76 (MQ=255)                                   ncGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              >  1:630678/2‑76 (MQ=255)                                    cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                             >  1:1154034/1‑76 (MQ=255)                                    cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                             >  1:830688/1‑76 (MQ=255)                                    cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                             >  1:87530/1‑76 (MQ=255)                                      aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                           <  1:532892/76‑1 (MQ=255)                                      aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                           <  1:587780/76‑1 (MQ=255)                                         atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                        >  1:413927/1‑76 (MQ=255)                                            cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGc                     <  2:209342/76‑1 (MQ=255)                                               aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATc                  <  1:567510/76‑1 (MQ=255)                                                gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt                 <  2:906872/76‑1 (MQ=255)                                                gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt                 <  2:921351/76‑1 (MQ=255)                                                 tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                <  2:864389/76‑1 (MQ=255)                                                 tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                <  2:693358/76‑1 (MQ=255)                                                 tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                >  2:45215/1‑76 (MQ=255)                                                  ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg               >  2:1320947/1‑76 (MQ=255)                                                  ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg               <  2:451598/76‑1 (MQ=255)                                                  ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg               >  2:324067/1‑76 (MQ=255)                                                  ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg               <  2:1258516/76‑1 (MQ=255)                                                   tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGc              <  1:1143278/76‑1 (MQ=255)                                                    gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg             <  1:1445974/76‑1 (MQ=255)                                                     aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa            >  2:1511997/1‑76 (MQ=255)                                                     aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa            >  1:429641/1‑76 (MQ=255)                                                      aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGaa           <  2:1608002/76‑1 (MQ=255)                                                       tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGCGACCGTGTGCATCTGGCGAat          <  2:211189/76‑1 (MQ=255)                                                        aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGGCTGACCGCTCGTGACCGTGTGCATCTGGCGAata         >  1:45976/1‑76 (MQ=255)                                                         cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat        <  1:1280681/76‑1 (MQ=255)                                                           gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa      >  1:538558/1‑76 (MQ=255)                                                           gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa      >  1:1112154/1‑76 (MQ=255)                                                            aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     <  2:1553726/76‑1 (MQ=255)                                                            aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     <  2:1483783/76‑1 (MQ=255)                                                            aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     <  1:1113930/76‑1 (MQ=255)                                                               agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATgcg  <  1:1258462/76‑1 (MQ=255)                                                                  |                                                                         GGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762676‑3762813

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