Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,762,741 A→T 100% K662I (AAA→ATA)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,762,7410A→T100.0% 166.9 481985spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  new base (25/20):  ref base (2/0):  total (27/21)

CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGA  >  REL606/3762708‑3762776
                                 |                                   
cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAg                                   <  1:7868171/36‑1 (MQ=255)
  tCGAATATTCAAAGTTTGAATACGGAAGAGATAGAt                                 >  1:5299677/1‑36 (MQ=255)
   cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                                <  1:6185151/36‑1 (MQ=255)
    gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                               <  1:3609886/36‑1 (MQ=255)
     aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                              <  1:580964/36‑1 (MQ=255)
      atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                             <  1:2564652/36‑1 (MQ=255)
       tatTCAAAGTTTGAATACGGAAGAGAAAGATGGTcg                            >  1:4910740/1‑36 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                           >  1:5940582/1‑36 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                           >  1:6111153/1‑36 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                           >  1:5623634/1‑36 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                           >  1:7023167/1‑36 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                           >  1:2765901/1‑36 (MQ=255)
         ttCACAGTTTGAATACGGAAGAGATAGATGGTcgcg                          <  1:7399528/36‑1 (MQ=37)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                          >  1:2742566/1‑36 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                          <  1:4096676/36‑1 (MQ=255)
          tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                         <  1:14445/36‑1 (MQ=255)
          tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                         >  1:263992/1‑36 (MQ=255)
          tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                         >  1:6574115/1‑36 (MQ=255)
          tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                         >  1:5182962/1‑36 (MQ=255)
            aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                       >  1:288840/1‑36 (MQ=255)
            aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                       >  1:1505689/1‑36 (MQ=255)
               gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                    >  1:7580935/1‑36 (MQ=255)
               gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                    >  1:1534302/1‑36 (MQ=255)
                tttGAATACGGAAGAGATAGATGGTCGCGTCTCCAg                   >  1:3769613/1‑36 (MQ=37)
                tttGAATACGGAAGAGATAGATGGTCGCGTCTACAg                   <  1:1092976/36‑1 (MQ=255)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc                  <  1:4972136/36‑1 (MQ=255)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc                  >  1:6173403/1‑36 (MQ=255)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc                  <  1:5229916/36‑1 (MQ=255)
                 ttGAATACGGAAGAGAAAGATGGTCGCGTCTACAgc                  >  1:2652738/1‑36 (MQ=255)
                  tGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg                 <  1:5105658/36‑1 (MQ=255)
                  tGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg                 <  1:2972517/36‑1 (MQ=255)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc                >  1:1539441/1‑36 (MQ=255)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc                >  1:5309814/1‑36 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt              >  1:6121313/1‑36 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt              <  1:7346706/36‑1 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt              <  1:8633/36‑1 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt              >  1:4545517/1‑36 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt              >  1:1312962/1‑36 (MQ=255)
                      tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt             <  1:5709096/36‑1 (MQ=255)
                      tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt             <  1:3632260/36‑1 (MQ=255)
                      tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt             <  1:7651127/36‑1 (MQ=255)
                        cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa           >  1:2991196/1‑36 (MQ=255)
                         ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt          >  1:6638452/1‑36 (MQ=255)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt         >  1:4231829/1‑36 (MQ=255)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt         >  1:1308536/1‑36 (MQ=255)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt         <  1:1281938/36‑1 (MQ=255)
                             gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGt      <  1:3420855/36‑1 (MQ=255)
                              agaGAGATGGTCGCGTCTACAGCGCCTTTATTCGTc     <  1:36726/36‑1 (MQ=255)
                                 aaGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa  <  1:2648575/36‑1 (MQ=255)
                                 |                                   
CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGA  >  REL606/3762708‑3762776

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG