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%	175.4 / NA	45	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 (25/20);  total (25/20)

CGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG  >  REL606/3762672‑3762811                                                                      |                                                                       cgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                                                                  <  1:830122/76‑1 (MQ=255) cgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                                                                  >  1:989837/1‑76 (MQ=255)  gcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGt                                                                 <  2:103279/76‑1 (MQ=255)    tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                               <  2:990018/76‑1 (MQ=255)    tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                               >  1:961531/1‑76 (MQ=255)    tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                               >  2:1157353/1‑76 (MQ=255)     ggCAAACNTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                                                              >  2:1187529/1‑76 (MQ=255)         aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                                                          <  2:694642/76‑1 (MQ=255)         aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                                                          <  2:461738/76‑1 (MQ=255)          aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTAGCGTCTa                                                         <  1:245686/76‑1 (MQ=255)           ccTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                                                        <  1:637270/76‑1 (MQ=255)            cTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                                                       <  1:407929/76‑1 (MQ=255)             tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg                                                      <  1:250821/76‑1 (MQ=255)                      cAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                             >  1:756737/1‑76 (MQ=255)                         ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                          >  2:1396161/1‑76 (MQ=255)                             caccacGACTTCGANTATTCAAAGCTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGAGCCTTTATTCGTCTg                                      >  1:7785/1‑76 (MQ=255)                             caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg                                      <  1:970060/76‑1 (MQ=255)                              accacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa                                     <  2:990954/76‑1 (MQ=255)                               ccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAc                                    <  2:1308455/76‑1 (MQ=255)                                cacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAcc                                   >  2:1351831/1‑76 (MQ=255)                                    aCTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTc                               <  2:110895/76‑1 (MQ=255)                                    aCTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTc                               >  2:1119360/1‑76 (MQ=255)                                       tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                            >  1:458025/1‑76 (MQ=255)                                       tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                            <  2:843832/76‑1 (MQ=255)                                        cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTTGTCTGACCGCTCGTGa                           >  1:249968/1‑76 (MQ=255)                                          aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                         >  1:206853/1‑76 (MQ=255)                                            taatCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgt                       >  1:1361204/4‑76 (MQ=255)                                             atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                      >  2:1326617/1‑76 (MQ=255)                                             atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                      >  2:604032/1‑76 (MQ=255)                                             atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                      >  2:672806/1‑76 (MQ=255)                                               tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGtgtg                    <  1:604034/76‑1 (MQ=255)                                               tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGngtg                    <  1:1326619/76‑1 (MQ=255)                                                 aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCa                  >  2:703941/1‑76 (MQ=255)                                                    gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAACGGCTTTATTCNTCTGACCGCTCGTGACCGTGTGCATCt               <  1:1298741/76‑1 (MQ=255)                                                     tttGAATACGGAAGCGATAGATGGTCGCGGCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg              <  2:25080/76‑1 (MQ=255)                                                      ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg             >  2:687929/1‑76 (MQ=255)                                                       tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACAGCTCGTGACCGTGTGCATCTGGc            >  2:437490/1‑76 (MQ=255)                                                         aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa          >  2:144349/1‑76 (MQ=255)                                                         aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa          <  1:1236601/76‑1 (MQ=255)                                                           tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat        >  2:887885/1‑76 (MQ=255)                                                           tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat        <  1:1009147/76‑1 (MQ=255)                                                            aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAata       >  2:264225/1‑76 (MQ=255)                                                             cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat      >  2:744905/1‑76 (MQ=255)                                                                aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt   >  2:786229/1‑76 (MQ=255)                                                                 agagATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATg  >  2:721277/1‑76 (MQ=255)                                                                      |                                                                       CGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG  >  REL606/3762672‑3762811

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