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
*REL6063,762,7410AT100.0% 75.8 / NA 21K662I (AAA→ATA) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (9/12);  total (9/12)

GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762668‑3762813
                                                                         |                                                                        
ggTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAg                                                                        <  1:307738/76‑1 (MQ=255)
         gCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                                                               >  1:464298/1‑76 (MQ=255)
         gCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                                                               >  2:345167/1‑76 (MQ=255)
           aaaCCTGACGGCGGGTATCAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                                                             <  1:8957/76‑1 (MQ=255)
                   cggcggCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc                                                     >  1:350371/1‑76 (MQ=255)
                       ggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCttt                                                 <  2:350371/76‑1 (MQ=255)
                            ttAAAAGCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                            >  1:791085/1‑76 (MQ=255)
                              aaCACAACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTc                                          >  2:236169/1‑76 (MQ=255)
                                caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg                                        <  1:197682/76‑1 (MQ=255)
                                        cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg                                <  1:164862/76‑1 (MQ=255)
                                            gAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAc                            >  1:760840/1‑76 (MQ=255)
                                                atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTCGAGCGCCCTTCTTCGTCTGACCTCTCNTGACCgtg                        >  2:133950/1‑76 (MQ=255)
                                                  tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGtgtg                      <  2:799072/76‑1 (MQ=255)
                                                    aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGACTTTANGCGTCAGACCGCTCGTGACCGTGTGCa                    <  2:3251/76‑1 (MQ=255)
                                                     aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCAt                   <  2:278830/76‑1 (MQ=255)
                                                       gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt                 <  1:629501/76‑1 (MQ=255)
                                                       gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt                 >  1:593809/1‑76 (MQ=255)
                                                          tGNATACGTAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTNGc              >  2:691529/1‑76 (MQ=255)
                                                           gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg             <  2:464298/76‑1 (MQ=255)
                                                           gAATACGGAAGAGATAGATGGTCGCGTATACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg             <  1:83383/76‑1 (MQ=255)
                                                                cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat        <  1:677225/76‑1 (MQ=255)
                                                                      agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATgcg  >  2:396160/1‑76 (MQ=255)
                                                                         |                                                                        
GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762668‑3762813

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