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% 148.8 411985spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  new base (22/19):  ref base (0/0):  total (22/19)

CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762708‑3762773
                                 |                                
cTTCGAATATTCAAAGTTTGAATACGGAAGAGaaat                                >  1:2115776/1‑35 (MQ=255)
cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAg                                <  1:1028561/36‑1 (MQ=255)
 ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                               <  1:1363168/36‑1 (MQ=255)
  tCGAATATTCAAAGTTTGAATACGGAAGAGATAGAt                              >  1:3853032/1‑36 (MQ=255)
  tCGAATATTCAAAGTTTGAATACGGAAGAGATAGAt                              <  1:5971621/36‑1 (MQ=255)
   cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             >  1:2391502/1‑36 (MQ=255)
     aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           >  1:633861/1‑36 (MQ=255)
      atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                          <  1:4802687/36‑1 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        <  1:5432118/36‑1 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        >  1:5965647/1‑36 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        <  1:1575410/36‑1 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       <  1:2469417/36‑1 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       >  1:308562/1‑36 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       >  1:6016460/1‑36 (MQ=255)
          tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                      <  1:4515178/36‑1 (MQ=255)
            aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    >  1:4938590/1‑36 (MQ=255)
            aaaGTTTGAATACGGAAGAGATACATGGTCGCGTCt                    >  1:2919936/1‑36 (MQ=38)
             aaGTTTGAATACGGAAGGGATAGATGGTCGCGTCTa                   <  1:3307129/36‑1 (MQ=38)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  >  1:1410660/1‑36 (MQ=255)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  >  1:3334921/1‑36 (MQ=255)
               gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                 >  1:4965243/1‑36 (MQ=255)
                  tGAATACGGAAGACATAGATGGTCGCGTCTACAgcg              >  1:4304847/1‑36 (MQ=38)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             <  1:7379247/36‑1 (MQ=255)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            >  1:6135381/1‑36 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           <  1:6452917/36‑1 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           >  1:1696353/1‑36 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           <  1:5822832/36‑1 (MQ=255)
                      tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt          >  1:1577485/1‑36 (MQ=255)
                       aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCttt         <  1:6300317/36‑1 (MQ=255)
                       aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCttt         <  1:1128006/36‑1 (MQ=255)
                         ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       <  1:5279185/36‑1 (MQ=255)
                         ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       <  1:223867/36‑1 (MQ=255)
                         ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       >  1:1340605/1‑36 (MQ=255)
                         ggAAGAGATAGATGGGCGCGTCTACAGCGCCTTTAt       >  1:5428085/1‑36 (MQ=38)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      >  1:7306555/1‑36 (MQ=255)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      <  1:900661/36‑1 (MQ=255)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      <  1:4453191/36‑1 (MQ=255)
                           aaGAGATAGATGGTCGCGTCTACATCGCCTTTATTc     >  1:2053363/1‑36 (MQ=37)
                            agagATAGATGGTCGCGTCTACAGCGCCTTTATTCg    >  1:6697046/1‑36 (MQ=255)
                             gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGt   <  1:4855960/36‑1 (MQ=255)
                             gagaTAGATGGTCGCGTCTACAGCGCCTGTATTCGt   >  1:34819/1‑36 (MQ=37)
                              agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  >  1:1545289/1‑36 (MQ=255)
                                 |                                
CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762708‑3762773

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG