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 change freq score reads annotation genes product
*REL6063,762,7410A→T100.0% 119.7 401985spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  new base (24/15):  ref base (0/0):  total (24/16)

CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762708‑3762773
                                 |                                
cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAg                                >  1:1541789/1‑36 (MQ=255)
cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAg                                >  2:233318/1‑36 (MQ=255)
  tCGAATATTCAAAGTTTGAATACGGAAGAGATAGAt                              <  1:349403/36‑1 (MQ=255)
   cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             <  1:51621/36‑1 (MQ=255)
   cGAATAATCAAAGTTTGAATACGGAAGAGATAGATg                             <  1:1429063/36‑1 (MQ=37)
   cGAAGAGTCAGAGTTGGAATACGGAAGAGATAGATg                             <  1:936246/36‑1 (MQ=12)
    gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            >  2:413795/1‑36 (MQ=255)
    gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            >  2:1813553/1‑36 (MQ=255)
      atatTCAAAGTTTGAATACGGAAGAGATAGATGGTc                          >  1:2959497/1‑36 (MQ=255)
      atatTCAAAGTTTGAATACCGAAGAGATACATggcc                          >  1:242603/1‑34 (MQ=21)
         ttCAAAGTTTGAATACGGAACAGGTAGATGGTcggt                       >  2:215954/1‑34 (MQ=21)
           cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                     <  2:528769/36‑1 (MQ=255)
           cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTc                     >  2:2608913/1‑36 (MQ=255)
            aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    <  2:1383509/36‑1 (MQ=255)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  <  1:2607106/36‑1 (MQ=255)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  <  2:2988054/36‑1 (MQ=255)
              aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  <  1:1407564/36‑1 (MQ=255)
              aGTTTGAATACGGAAGAGATAGATGGTCGCCTCTAc                  >  2:3043560/1‑36 (MQ=37)
               gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                 >  2:1465587/1‑36 (MQ=255)
                tttGAATACGGAAGAGATAGATGGTCGCGTCTACAg                >  2:2592203/1‑36 (MQ=255)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTATAgc               >  2:544213/1‑36 (MQ=37)
                 ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc               <  1:138356/36‑1 (MQ=255)
                  tGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg              >  2:357014/1‑36 (MQ=255)
                  tGAATACGGAAGAGAGAGATGGTCGCGTCTACAgcg              <  1:59172/36‑1 (MQ=255)
                   gAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc             <  2:3010471/36‑1 (MQ=255)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            >  2:2438558/1‑36 (MQ=255)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            >  1:2199246/1‑36 (MQ=255)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            >  2:3074487/1‑36 (MQ=255)
                    aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            >  1:1610083/1‑36 (MQ=255)
                     aTACGGAAGAGATAGATGGTCGCGTCTACAGTGCCt           >  2:353043/1‑36 (MQ=37)
                     aTACGAAAGAGATAGATGGTCGCGCCTACAGCGCCt           <  2:2906956/36‑1 (MQ=25)
                      tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt          >  1:3666513/1‑36 (MQ=255)
                       aCGGAAGAGATAGATGGTCGCGTCTACCACGGCttt         >  2:3125476/1‑36 (MQ=19)
                         ggAAGAGGTAGATGGTCGCGTCTACAGCGCCTTTAt       <  1:2404201/36‑1 (MQ=38)
                         ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       >  2:560346/1‑36 (MQ=255)
                          gAAGAGATAGATGGGCGCGTCTACAGCGCCTTTAtt      <  2:3720893/36‑1 (MQ=38)
                            agagATAGATGGTCGCGTCTACAGCGCCTTTATTCg    >  2:72052/1‑36 (MQ=255)
                            agagATAGATGGTCGCGTCTACAGCGCCTTTATTCg    >  1:1637628/1‑36 (MQ=255)
                              agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  <  2:2075044/36‑1 (MQ=255)
                              agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  >  2:1230231/1‑36 (MQ=255)
                                 |                                
CTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762708‑3762773

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG