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.6 / NA 22K662I (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 (11/11);  total (11/11)

GTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCG  >  REL606/3762669‑3762813
                                                                        |                                                                        
gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                                                                       <  2:554489/76‑1 (MQ=255)
 tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat                                                                      <  1:747633/76‑1 (MQ=255)
          aaaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                                                             >  1:738003/1‑76 (MQ=255)
           aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                                                            <  2:471510/76‑1 (MQ=255)
                     cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt                                                  >  2:388339/1‑76 (MQ=255)
                        cAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                               >  1:387848/1‑76 (MQ=255)
                         aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt                                              <  1:156503/76‑1 (MQ=255)
                           ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                            >  2:505707/1‑76 (MQ=255)
                            tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                           >  2:542382/1‑76 (MQ=255)
                                 ccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAc                                      >  1:323809/1‑76 (MQ=255)
                                  cacGACTTCGAATNTTCAAAGTTTNAATACGGACGAGATNGATGGTCGCGTCTACNGCGCCTNTATTCGTCTGAcc                                     >  2:675347/1‑76 (MQ=255)
                                         tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                              <  1:60730/76‑1 (MQ=255)
                                              tagtCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgt                         <  1:367979/73‑1 (MQ=255)
                                                ttCAAAGTTTGAATACGGAAGAGATAGGTGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTTACCgtgt                       <  2:705281/76‑1 (MQ=255)
                                                ttCAAAGTTTGAATACGGAAGAGATAGATGGTCTCGTCTACAGCGCCATTATTCGTCTGACCGCTCGTGACCgtgt                       <  2:742060/76‑1 (MQ=255)
                                                 tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGtgtg                      >  1:522142/1‑76 (MQ=255)
                                                     aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATc                  >  2:419394/1‑76 (MQ=255)
                                                      gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt                 <  1:497782/76‑1 (MQ=255)
                                                          gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg             >  1:470427/1‑76 (MQ=255)
                                                             tACGGAAGAGATGGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat          <  2:19708/76‑1 (MQ=255)
                                                             tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat          >  2:540972/1‑76 (MQ=255)
                                                                 gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa      <  1:435254/76‑1 (MQ=255)
                                                                 gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa      <  1:390065/76‑1 (MQ=255)
                                                                     agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATgcg  >  1:786545/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 < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: