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%	262.9 / NA	72	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 (39/33);  total (39/33)

GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762668‑3762813                                                                          |                                                                         ggTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAg                                                                        <  1:1261513/76‑1 (MQ=255)  gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                                                                       <  1:1658976/76‑1 (MQ=255)  gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                                                                       <  2:329946/76‑1 (MQ=255)        tGGCAAACNTGACGGCGGCAATTAACACCACGACTTCGAATATTCAACGTTTGAATACGGAAGAGATAGATGGTcg                                                                 >  1:1760505/1‑76 (MQ=255)        tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                                 >  1:688563/1‑76 (MQ=255)        tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                                 <  2:951406/76‑1 (MQ=255)        tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                                 >  1:982194/1‑76 (MQ=255)        tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                                 <  2:1731300/76‑1 (MQ=255)        tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                                 >  2:1601523/1‑76 (MQ=255)         ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                                                                >  2:1416642/1‑76 (MQ=255)             aaCCTGACGGCGGCAATTAACATCANGTCTTCGAATGTTCAAAGTTANAATTCGGAAGAGATAGATGGTCGCGTCt                                                            <  1:265700/76‑1 (MQ=255)              aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                                                           <  1:1225822/76‑1 (MQ=255)               ccTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                                                          <  2:318149/76‑1 (MQ=255)               ccTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                                                          <  1:375150/76‑1 (MQ=255)               ccTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                                                          <  2:747176/76‑1 (MQ=255)                 tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg                                                        >  2:216207/1‑76 (MQ=255)                 tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg                                                        >  1:781579/1‑76 (MQ=255)                  gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgc                                                       >  1:931088/1‑76 (MQ=255)                    cggcggCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc                                                     <  2:242738/76‑1 (MQ=255)                        ggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACGNCGCCttt                                                 <  2:1469794/76‑1 (MQ=255)                         gtaaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa                                                <  2:632603/74‑1 (MQ=255)                          cAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                               >  1:1017110/1‑76 (MQ=255)                          cAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                               >  1:1156228/1‑76 (MQ=255)                           aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt                                              >  2:864040/1‑76 (MQ=255)                           aaTTAACACCACCACCTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCNACAGCGCCTTTAtt                                              >  1:1570384/1‑76 (MQ=255)                            aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc                                             >  1:206251/1‑76 (MQ=255)                            aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc                                             <  2:950897/76‑1 (MQ=255)                             ttAACACCAGGAATTTGAACATTCAAAGTGTGAGTACGCAAGATATACTGGGTCGCGTCTACAGCGTCTTTGTTCg                                            >  1:1771958/1‑76 (MQ=255)                             ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                            >  1:334907/1‑76 (MQ=255)                             ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                            >  1:1231189/1‑76 (MQ=255)                             ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                            >  1:847509/1‑76 (MQ=255)                             ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                            >  1:1271231/1‑76 (MQ=255)                               aaCACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTc                                          >  1:1345974/1‑76 (MQ=255)                               aaCACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTc                                          >  2:1401193/1‑76 (MQ=255)                                acacCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCt                                         <  2:665363/76‑1 (MQ=255)                                acacCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCt                                         <  1:718398/76‑1 (MQ=255)                                 caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACCGCGCCTTTATTCGTCTg                                        >  1:40139/1‑76 (MQ=255)                                 caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg                                        <  1:839547/76‑1 (MQ=255)                                  accacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa                                       <  1:1281453/76‑1 (MQ=255)                                    cacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAcc                                     >  2:743386/1‑76 (MQ=255)                                     acGACTTCGAATATTTAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCg                                    <  1:1506177/76‑1 (MQ=255)                                         cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg                                <  2:679878/76‑1 (MQ=255)                                         cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg                                >  1:657732/1‑76 (MQ=255)                                          ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCgc                               <  1:552090/76‑2 (MQ=255)                                          ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGt                               >  1:539558/1‑76 (MQ=255)                                           tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              >  2:1205392/1‑76 (MQ=255)                                            cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                             >  1:1370456/1‑76 (MQ=255)                                             gAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAc                            <  2:600338/76‑1 (MQ=255)                                               atatTCAAAGTTTGAATACGGAAGCGATNAATGGTCGCGACTACAGCGACNTTATTCGTCTGACCGCTCGTGACCg                          >  1:266737/1‑76 (MQ=255)                                                   tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGtgtg                      >  1:1394323/1‑76 (MQ=255)                                                   tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGtgtg                      <  1:929444/76‑1 (MQ=255)                                                     aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCa                    >  2:649093/1‑76 (MQ=255)                                                       aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATACGTCTGACCGCTCGTGACCGTGTGCATc                  >  2:1369172/1‑76 (MQ=255)                                                        gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt                 <  2:950333/76‑1 (MQ=255)                                                        gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGCGTGCATCt                 >  1:249322/1‑76 (MQ=255)                                                         tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                <  2:982193/76‑1 (MQ=255)                                                         tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                <  2:688562/76‑1 (MQ=255)                                                         tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                >  1:1183462/1‑76 (MQ=255)                                                         tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg                <  2:1511536/76‑1 (MQ=255)                                                          ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg               >  1:104766/1‑76 (MQ=255)                                                          ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg               <  1:1305373/76‑1 (MQ=255)                                                           tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGc              >  2:1408788/1‑76 (MQ=255)                                                             aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa            >  2:728754/1‑76 (MQ=255)                                                             aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa            >  2:704229/1‑76 (MQ=255)                                                               tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat          >  1:978067/1‑76 (MQ=255)                                                               tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAACGTGTGCATCTGGCGAat          <  1:1721570/76‑1 (MQ=255)                                                                 cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat        <  2:1304735/76‑1 (MQ=255)                                                                  ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc       >  2:109688/1‑76 (MQ=255)                                                                  ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc       >  1:474525/1‑76 (MQ=255)                                                                  ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc       >  2:897962/1‑76 (MQ=255)                                                                   gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa      >  2:346115/1‑76 (MQ=255)                                                                    aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     <  2:215927/76‑1 (MQ=255)                                                                    aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     <  2:928453/76‑1 (MQ=255)                                                                    aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     >  2:1367170/1‑76 (MQ=255)                                                                    aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt     <  2:344827/76‑1 (MQ=255)                                                                       agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATgcg  >  1:332106/1‑76 (MQ=255)                                                                       agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAANATCATgcg  <  1:1761708/76‑1 (MQ=255)                                                                          |                                                                         GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762668‑3762813

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