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%	147.2 / NA	38	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 (18/20);  total (18/20)

GTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762669‑3762813                                                                         |                                                                         gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTNGAATATTCAAAGTTTGAATGCGGAAGAGATAGa                                                                       <  2:1292043/76‑1 (MQ=255) gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                                                                       <  1:1147623/76‑1 (MQ=255)  tGCGCTGGCGAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAATagatagat                                                                      <  1:49489/76‑1 (MQ=255)   gcggTGGCGAACCTGGCGGCGGCAATTAACACCACGACGTCGAATATTCAAAGTTTGAATACGGATGAGATAGATg                                                                     <  1:1420191/76‑1 (MQ=255)   gcgcTGGCAAACCTGACGGCGGCAATTAACCCCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                                                                     <  1:721287/76‑1 (MQ=255)      cTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTc                                                                  >  1:1121619/1‑76 (MQ=255)       tGGCCATCCTGACGGCGGCAATTAACACCACGACTTCGGATATTCAAATTTTGAATACGGAAGAGATAGATGGTcg                                                                 <  1:1508721/76‑1 (MQ=255)       tGGCAAACCTGAAGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                                 >  2:1524987/1‑76 (MQ=255)           aaaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                                                             >  1:554960/1‑76 (MQ=255)            aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                                                            <  2:430284/76‑1 (MQ=255)             aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                                                           <  2:928329/76‑1 (MQ=255)             aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                                                           <  2:424340/76‑1 (MQ=255)                 gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgc                                                       <  2:793473/76‑1 (MQ=255)                 gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgc                                                       <  2:1197538/76‑1 (MQ=255)                               acacCACGACTTCGAATATTCAAAGTTTGTATACGGTAGAGATAGATGGTCGCGTCTNCAGCGCCTTTATTCGTCt                                         <  1:1468108/76‑1 (MQ=255)                               acacCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCt                                         >  2:625466/1‑76 (MQ=255)                                 accacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa                                       <  2:801978/76‑1 (MQ=255)                                 accacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa                                       <  1:325683/76‑1 (MQ=255)                                  ccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAc                                      <  1:671871/76‑1 (MQ=255)                                       aCTTCGAATATTCAAAGGTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTc                                 <  1:177644/76‑1 (MQ=255)                                          tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              >  2:694167/1‑76 (MQ=255)                                            gAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAc                            <  1:892118/76‑1 (MQ=255)                                             aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                           >  2:922362/1‑76 (MQ=255)                                             aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                           >  2:1342028/1‑76 (MQ=255)                                              atatTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCg                          >  1:951260/1‑76 (MQ=255)                                                      aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATc                  <  1:810631/76‑1 (MQ=255)                                                        tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                <  2:380976/76‑1 (MQ=255)                                                        tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                <  2:1028753/76‑1 (MQ=255)                                                         ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg               <  2:159350/76‑1 (MQ=255)                                                          tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGc              >  1:903215/1‑76 (MQ=255)                                                              tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat          >  2:184518/1‑76 (MQ=255)                                                              tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat          >  2:665562/1‑76 (MQ=255)                                                              tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat          >  2:924830/1‑76 (MQ=255)                                                               aCGGAAGATATAGATGGTCGCGCCTACAGCGCNTTAATTCGNCTGACCGCTCCTGACCGTGNGCATCTGGCGAata         >  2:1271210/1‑76 (MQ=255)                                                               aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAata         >  2:829802/1‑76 (MQ=255)                                                                  gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa      >  2:697481/1‑76 (MQ=255)                                                                   aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     >  1:892992/1‑76 (MQ=255)                                                                   aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     >  1:892982/1‑76 (MQ=255)                                                                    agagATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATg    >  1:995558/1‑76 (MQ=255)                                                                      agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATgcg  >  2:536411/1‑76 (MQ=255)                                                                         |                                                                         GTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762669‑3762813

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