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% 177.8 / NA 44K662I (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 (29/15);  total (29/15)

TGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGC  >  REL606/3762670‑3762812
                                                                       |                                                                       
tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                                                                      <  1:532697/75‑1 (MQ=255)
 gcgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat                                                                     >  1:842446/1‑75 (MQ=255)
   gcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                                                                   >  1:674244/1‑75 (MQ=255)
   gcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                                                                   >  2:693792/1‑75 (MQ=255)
     tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTc                                                                 <  2:798396/75‑1 (MQ=255)
     tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTc                                                                 >  2:796678/1‑75 (MQ=255)
      ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                                                                >  1:1039529/1‑75 (MQ=255)
        cAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                                                              >  1:95276/1‑75 (MQ=255)
          aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                                                            >  1:1237268/1‑75 (MQ=255)
          aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                                                            <  1:703598/75‑1 (MQ=255)
          aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                                                            >  2:1140697/1‑75 (MQ=255)
           aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                                                           <  2:823317/75‑1 (MQ=255)
              tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                                                        >  1:249428/1‑75 (MQ=255)
               gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg                                                       >  2:776779/1‑75 (MQ=255)
               gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg                                                       >  1:801544/1‑75 (MQ=255)
                 cggcggCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg                                                     >  1:1614611/1‑75 (MQ=255)
                 cggcggCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg                                                     >  2:1303344/1‑75 (MQ=255)
                    cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt                                                  <  2:549764/75‑1 (MQ=255)
                     ggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt                                                 <  2:1104218/75‑1 (MQ=255)
                       cAATTACCCCCCCGACTCCG‑ATATTCCAAGTTTGAATTCGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                              <  2:1543182/75‑1 (MQ=255)
                        aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                              <  1:568310/75‑1 (MQ=255)
                        aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                              >  1:746215/1‑75 (MQ=255)
                         aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt                                             >  1:1165419/1‑75 (MQ=255)
                         aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt                                             >  1:845459/1‑75 (MQ=255)
                         aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt                                             >  1:1207286/1‑75 (MQ=255)
                           tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                           >  1:1439296/1‑75 (MQ=255)
                            aaCACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                          >  1:1413169/1‑75 (MQ=255)
                            aaCACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt                                          >  1:313057/1‑75 (MQ=255)
                              caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCt                                        <  2:37799/75‑1 (MQ=255)
                               accacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg                                       <  2:735291/75‑1 (MQ=255)
                                  acGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAcc                                    <  2:1225804/75‑1 (MQ=255)
                                       ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg                               <  2:1374619/75‑1 (MQ=255)
                                        tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGt                              >  2:754803/1‑75 (MQ=255)
                                           aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAc                           <  1:305413/75‑1 (MQ=255)
                                            atatTCAAAGTTTGAATACGGAAGAGATCGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                          >  2:161238/1‑75 (MQ=255)
                                            atatTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc                          >  2:229541/1‑75 (MQ=255)
                                              atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgt                        <  2:335666/75‑1 (MQ=255)
                                                tCAAAGTTTGAATACGGAAGAGATAGATGGCCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtgg                      >  2:838658/1‑74 (MQ=255)
                                                   aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCa                   >  2:1478786/1‑75 (MQ=255)
                                                      tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt                >  1:835154/1‑75 (MQ=255)
                                                            tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGaa          >  2:1191468/1‑75 (MQ=255)
                                                              cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAata        >  2:1033001/1‑75 (MQ=255)
                                                              cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAata        <  1:211989/75‑1 (MQ=255)
                                                                gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc      >  2:67993/1‑75 (MQ=255)
                                                                   gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATg   >  2:958554/1‑75 (MQ=255)
                                                                   gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATg   <  2:258428/75‑1 (MQ=255)
                                                                    agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATgc  <  1:776779/75‑1 (MQ=255)
                                                                       |                                                                       
TGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGC  >  REL606/3762670‑3762812

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