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% 242.8 701985spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  new base (38/28):  ref base (2/0):  total (42/28)

CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762708‑3762773
                                 |                                
cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAg                                <  1:4463951/36‑1 (MQ=255)
 ttCGAATATTCAAAGTTTGAATAAGGAAGAGAGAGa                               >  1:3743765/1‑36 (MQ=37)
  tCGAATATTCAAAGTTTGAATACGGAAGAGATAGAt                              >  1:8283386/1‑36 (MQ=255)
   cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             <  1:55622/36‑1 (MQ=255)
   cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             >  1:1562286/1‑36 (MQ=255)
   cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             >  1:3938429/1‑36 (MQ=255)
    gAATATTCAAAGTTTGAATCCGGAAGAGATAGATgg                            <  1:147533/36‑1 (MQ=37)
    gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            >  1:8256213/1‑36 (MQ=255)
     aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           >  1:7907901/1‑36 (MQ=255)
      atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                          >  1:6148085/1‑36 (MQ=255)
      atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                          <  1:2979533/36‑1 (MQ=255)
      atatTCAAAGTTTGAATACGGAAGAGAAAGATGGTc                          >  1:1881222/1‑36 (MQ=255)
       tatTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                         <  1:6432268/36‑1 (MQ=255)
       tatTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                         <  1:5760586/36‑1 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        >  1:7159349/1‑36 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        <  1:6407711/36‑1 (MQ=255)
        atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        >  1:5526213/1‑36 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       >  1:5231384/1‑36 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       <  1:3176090/36‑1 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       >  1:2400248/1‑36 (MQ=255)
         ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       >  1:5662027/1‑36 (MQ=255)
          tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                      <  1:4419852/36‑1 (MQ=255)
           cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                     >  1:4381111/1‑36 (MQ=255)
            aaaGTTTGAATACGGAAGAGATAGATGGTCTCGTCt                    <  1:1490926/36‑1 (MQ=37)
            aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    <  1:2463070/36‑1 (MQ=255)
            aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    <  1:3717901/36‑1 (MQ=255)
            aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    >  1:2352385/1‑36 (MQ=255)
             aaGTTTGAGTACGGAAGAGATAGATGGTCGCGTCTa                   >  1:3797829/1‑36 (MQ=37)
             aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                   >  1:496580/1‑36 (MQ=255)
             aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                   >  1:7410649/1‑36 (MQ=255)
              aGTTTGAATACGGAANAGATAGATGGTCGCGTCTAc                  <  1:2294629/36‑1 (MQ=39)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  >  1:759563/1‑36 (MQ=255)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  >  1:412637/1‑36 (MQ=255)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  >  1:4754135/1‑36 (MQ=255)
               gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                 >  1:521267/1‑36 (MQ=255)
               gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                 <  1:563193/36‑1 (MQ=255)
                tttGAATACGGAAGAGATAGATGGTCGCGTCTACAg                >  1:5486757/1‑36 (MQ=255)
                tttGAATACGGAAGAGATAGATGGTCGCGTCTACAg                <  1:1812627/36‑1 (MQ=255)
                tttGAATACGGAAGAGATAGATGGTCGCGTCTACAg                >  1:1203863/1‑36 (MQ=255)
                tttGAATACGGAAGAGATAGATGGTCGCGTCTACAg                <  1:6888934/36‑1 (MQ=255)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc               <  1:7587773/36‑1 (MQ=255)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc               >  1:3005589/1‑36 (MQ=255)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc               <  1:1629/36‑1 (MQ=255)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             >  1:5208237/1‑36 (MQ=255)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             >  1:4349616/1‑36 (MQ=255)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             <  1:3815351/36‑1 (MQ=255)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             >  1:685773/1‑36 (MQ=255)
                   gAATACGGAAGAGAAAGATGGTCGCGTCTACAgcgc             >  1:5732980/1‑36 (MQ=255)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAgcgaa            >  1:7484819/1‑34 (MQ=37)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            >  1:4149568/1‑36 (MQ=255)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            <  1:771187/36‑1 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           <  1:3572956/36‑1 (MQ=255)
                       aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCttt         <  1:4348012/36‑1 (MQ=255)
                        cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa        >  1:7720711/1‑36 (MQ=255)
                        cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa        <  1:8364845/36‑1 (MQ=255)
                        cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa        >  1:5757250/1‑36 (MQ=255)
                        cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa        >  1:4178623/1‑36 (MQ=255)
                         tgAAGGGATAGATGGTCGCGTCTACAGCGCCTTTAt       <  1:7827778/35‑1 (MQ=37)
                         ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       >  1:2855051/1‑36 (MQ=255)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      <  1:6819986/36‑1 (MQ=255)
                          gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      <  1:4757823/36‑1 (MQ=255)
                           aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc     >  1:3037882/1‑36 (MQ=255)
                           aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc     <  1:1084975/36‑1 (MQ=255)
                           aaGAGATAGATGGTCGCGTCGACAGCGCCTTTATTc     >  1:2532628/1‑36 (MQ=37)
                            agagATAGATGGTCGCGTCTACAGCGCCTTTATTCg    >  1:2619198/1‑36 (MQ=255)
                             gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGt   >  1:4231699/1‑36 (MQ=255)
                             gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGt   <  1:3286009/36‑1 (MQ=255)
                              agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  >  1:7442890/1‑36 (MQ=255)
                              agaTAGATGGTCGCGTCTACAGCGACTTTATTCGTc  >  1:7734120/1‑36 (MQ=37)
                              agaGAGATGGGCGCGTCTACAGCGCCTTTGTTAGTc  >  1:6900286/1‑36 (MQ=19)
                                 |                                
CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762708‑3762773

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG