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%	151.4 / NA	40	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 (20/20);  total (20/20)

GTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCA  >  REL606/3762669‑3762809                                                                         |                                                                     gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                                                                   <  1:788733/76‑1 (MQ=255)  tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat                                                                  <  1:988158/76‑1 (MQ=255)  tGCGCTGGCAAACCTGACGGCGGCAATAAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat                                                                  >  1:1277018/1‑76 (MQ=255)    cgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                                                                <  1:554490/76‑1 (MQ=255)       tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGNGATAGATGGTcg                                                             <  1:20488/76‑1 (MQ=255)       tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                             <  1:851105/76‑1 (MQ=255)               cTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                                                     <  1:645070/76‑1 (MQ=255)                     gcggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt                                               <  1:426797/76‑1 (MQ=255)                      cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt                                              <  2:437825/76‑1 (MQ=255)                       ggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCttt                                             <  2:479446/76‑1 (MQ=255)                       ggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCttt                                             <  2:1204906/76‑1 (MQ=255)                        gcAATTAACACCACGACTTCGAATATTCAAAGTGTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa                                            <  2:318336/76‑1 (MQ=255)                         cAATTAACACCTCGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                           <  2:454122/76‑1 (MQ=255)                            ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                        >  1:615059/1‑76 (MQ=255)                             tAACNCCACGACTTCGAATATTCTAAGTCTGCATACGGANGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                       <  1:158099/76‑1 (MQ=255)                             tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                       >  1:188328/1‑76 (MQ=255)                                caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg                                    >  1:526750/1‑76 (MQ=255)                                caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg                                    <  2:83746/76‑1 (MQ=255)                                         ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGt                           >  2:760490/1‑76 (MQ=255)                                          tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATNGTCGCGACTACAGCGCCTTTATTCGTCTGACCGCNCGTg                          >  1:159776/1‑76 (MQ=255)                                          tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                          >  2:968479/1‑76 (MQ=255)                                           cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                         >  1:1077412/1‑76 (MQ=255)                                           cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                         <  2:139064/76‑1 (MQ=255)                                           cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                         >  1:958944/1‑76 (MQ=255)                                           cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                         >  1:1077394/1‑76 (MQ=255)                                             aataatCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                       >  1:377755/1‑76 (MQ=255)                                             aaTATTCAAAGTTTGAATACGGAAGAGNTAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                       >  1:1295002/1‑76 (MQ=255)                                                atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                    <  1:758851/76‑1 (MQ=255)                                                 ttCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtgt                   <  2:994095/76‑1 (MQ=255)                                                  ttaaaGTTGGAATACGGAAGAGATANATGGTCGCGTCTACNGCGCCGTCATNCNTCTGACCGCGCGTGACCGtgtg                  >  2:22784/3‑76 (MQ=255)                                                   cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGCCCGTGTGc                 >  1:121643/1‑76 (MQ=255)                                                   cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGc                 >  2:303374/1‑76 (MQ=255)                                                     aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCAt               >  2:1199683/1‑76 (MQ=255)                                                      aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACACCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATc              >  2:85258/1‑76 (MQ=255)                                                         ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg           >  1:1154376/1‑76 (MQ=255)                                                            aaTGCGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa        <  1:306523/76‑1 (MQ=255)                                                               aCGGAAGAGCTAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAata     <  1:40192/76‑1 (MQ=255)                                                               aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAata     >  1:851665/1‑76 (MQ=255)                                                                cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGATCGCTCGTGACCGTGTGCATATGGCGAatat    <  1:1200718/76‑1 (MQ=255)                                                                cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat    <  2:411995/76‑1 (MQ=255)                                                                 ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc   >  2:1246547/1‑76 (MQ=255)                                                                 ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc   <  1:250185/76‑1 (MQ=255)                                                                 ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc   >  2:774548/1‑76 (MQ=255)                                                                  gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa  >  1:587508/1‑76 (MQ=255)                                                                         |                                                                     GTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCA  >  REL606/3762669‑3762809

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