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

TTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762709‑3762773
                                |                                
ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                               >  1:7749665/1‑36 (MQ=255)
ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                               >  1:5194425/1‑36 (MQ=255)
 tCGAATATTCAAAGTTTGAATACGGAAGAGATAGAt                              <  1:3093674/36‑1 (MQ=255)
  cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             <  1:653637/36‑1 (MQ=255)
  cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             >  1:3192630/1‑36 (MQ=255)
   gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            >  1:2406460/1‑36 (MQ=255)
   gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            <  1:797071/36‑1 (MQ=255)
   gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            >  1:2088841/1‑36 (MQ=255)
   gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            <  1:2593642/36‑1 (MQ=255)
    aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           <  1:6870824/36‑1 (MQ=255)
    aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           <  1:4741776/36‑1 (MQ=255)
    aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           >  1:7018007/1‑36 (MQ=255)
     atgttCAAAGTTTGAATACGGAAGAGATAGATGGTc                          >  1:7491995/4‑36 (MQ=37)
     atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                          >  1:167534/1‑36 (MQ=255)
     atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                          <  1:1230295/36‑1 (MQ=255)
     atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                          <  1:6090156/36‑1 (MQ=255)
     atatTCAAAGTTAGAATACGGAAGAGATAGATGGTc                          >  1:4504190/1‑36 (MQ=37)
      tatTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                         >  1:3868745/1‑36 (MQ=255)
      tatTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                         >  1:2430085/1‑36 (MQ=255)
      tatTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                         <  1:2208383/36‑1 (MQ=255)
       atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        <  1:90952/36‑1 (MQ=255)
          cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                     >  1:719078/1‑36 (MQ=255)
          cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                     >  1:5544521/1‑36 (MQ=255)
          cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                     >  1:727821/1‑36 (MQ=255)
           aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    >  1:4902943/1‑36 (MQ=255)
           aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    <  1:1358905/36‑1 (MQ=255)
            aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                   <  1:6039469/36‑1 (MQ=255)
            aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                   <  1:1723192/36‑1 (MQ=255)
             aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  <  1:1076209/36‑1 (MQ=255)
              gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                 <  1:7542726/36‑1 (MQ=255)
              gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                 <  1:4174161/36‑1 (MQ=255)
              gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                 <  1:3534316/36‑1 (MQ=255)
               tttGAATACGGAAGAGATAGATGGTCGCGTCTACAg                <  1:132355/36‑1 (MQ=255)
                ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc               >  1:270663/1‑36 (MQ=255)
                 tGAATACGGAAGAGATAGATGGTCGCGTCTATAgcg              >  1:6426706/1‑36 (MQ=37)
                 tGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg              <  1:3228951/36‑1 (MQ=255)
                 tGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg              <  1:4799631/36‑1 (MQ=255)
                 tGAATACGGAAGAGATAGATGGTCGAGTCTACAgcg              >  1:3406210/1‑36 (MQ=37)
                  gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             >  1:2964470/1‑36 (MQ=255)
                  gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             <  1:7818611/36‑1 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           >  1:1904734/1‑36 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           >  1:1996709/1‑36 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           <  1:7504046/36‑1 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           >  1:6945219/1‑36 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           >  1:3431303/1‑36 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           <  1:3105144/36‑1 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           <  1:2787840/36‑1 (MQ=255)
                    aTACGGAAGAGATAGATGGTCGCGTATACAGCGCCt           >  1:1286552/1‑36 (MQ=37)
                      aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCttt         >  1:2078402/1‑36 (MQ=255)
                       cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa        <  1:3875317/36‑1 (MQ=255)
                        ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       >  1:1461489/1‑36 (MQ=255)
                        ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       >  1:5429543/1‑36 (MQ=255)
                         gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      >  1:2487910/1‑36 (MQ=255)
                         gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      >  1:765416/1‑36 (MQ=255)
                         gAAGAGATAGATGGTCGCGTCTACAGCGCATTTAtt      >  1:2811571/1‑36 (MQ=38)
                          aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc     >  1:1051837/1‑36 (MQ=255)
                          aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc     <  1:2001243/36‑1 (MQ=255)
                          aaGAGATAGATGGCCGGGTCTACAGCGCCTTTAttt     >  1:2155648/1‑35 (MQ=25)
                            gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGt   >  1:5662519/1‑36 (MQ=255)
                            gagaTAGATGGTCGCGTCTACAGCGCCTTTATTCGt   <  1:5314735/36‑1 (MQ=255)
                             agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  <  1:6721195/36‑1 (MQ=255)
                             agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  <  1:2327906/36‑1 (MQ=255)
                             agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  <  1:5341156/36‑1 (MQ=255)
                                |                                
TTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762709‑3762773

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