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

CGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762711‑3762773
                              |                                
cGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             >  1:5843955/1‑36 (MQ=255)
agAATATTCAAAGTTTGAATACGGAAGAGATAGATg                             <  1:1294295/35‑1 (MQ=255)
 gAATATTCAAAGTTTGAATACGGAAGAGATAGATgg                            >  1:546493/1‑36 (MQ=255)
  aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           <  1:1189641/36‑1 (MQ=255)
  aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           <  1:72267/36‑1 (MQ=255)
  aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           >  1:2363772/1‑36 (MQ=255)
  aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGt                           <  1:4949405/36‑1 (MQ=255)
    tatTCAAAGTTTGAATACGGAAGAGATAGATGGTcg                         <  1:4563410/36‑1 (MQ=255)
     atTCAAAGTTTGAATACGGAAGAGATAGATGGTcgc                        <  1:3613471/36‑1 (MQ=255)
      ttCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                       <  1:1123827/36‑1 (MQ=255)
        cAAAGTTTGAAGACGGAAGAGATAGATGGTCGCGTc                     <  1:2145213/36‑1 (MQ=37)
         aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCt                    <  1:5665493/36‑1 (MQ=255)
          aaGTTTGAATATGGAAGAGATAGATGGTCGCGTCTa                   <  1:3766991/36‑1 (MQ=37)
          aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                   >  1:5170764/1‑36 (MQ=255)
          aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                   >  1:4658389/1‑36 (MQ=255)
          aaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                   >  1:3750786/1‑36 (MQ=255)
           cgTTTGAATGCGGAAGAGATAGATGGTCGCGTCTAc                  <  1:1253112/35‑1 (MQ=37)
           aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                  >  1:689558/1‑36 (MQ=255)
            gTTTGAATACGGAAGAGAGAGATGGTCGCGTCTaaa                 >  1:3164979/1‑34 (MQ=38)
              ttGAATACGGAAGAGATAGATGGTCGCGTCTACAgc               >  1:93079/1‑36 (MQ=255)
               tGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg              >  1:3774134/1‑36 (MQ=255)
               tGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg              >  1:3893038/1‑36 (MQ=255)
                gAATACGGAGGAGATAGAGGTTCGCGTCTACAgcgc             <  1:1250646/36‑1 (MQ=19)
                 aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            <  1:3292965/36‑1 (MQ=255)
                 aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGcc            >  1:2331952/1‑36 (MQ=255)
                 aaTACGGAAGAGATAGATGGACGAGTCTACAgcgct            >  1:291088/1‑35 (MQ=21)
                  aTACGGACGCGATCGATGGTCGCGTCTACAGCGCCt           <  1:1242213/36‑1 (MQ=19)
                  aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt           >  1:5895058/1‑36 (MQ=255)
                   tACGGAAGAGATAGATGGTCGCTTCTACAGCGCCtt          <  1:5714016/36‑1 (MQ=37)
                   tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt          >  1:1722163/1‑36 (MQ=255)
                     cgtaaGAGATAGATGGTCGCGTCTACAGCGCCTTTa        <  1:704955/33‑1 (MQ=38)
                     cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa        <  1:2142539/36‑1 (MQ=255)
                     cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa        <  1:3892936/36‑1 (MQ=255)
                      ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt       >  1:5175601/1‑36 (MQ=255)
                       gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt      <  1:2235807/36‑1 (MQ=255)
                        aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc     <  1:77631/36‑1 (MQ=255)
                         agagATAGATGGTCGCGTCTACAGCGCCTTTATTCg    <  1:4270322/36‑1 (MQ=255)
                         agagATAGATGGTCGCGTCTACAGCGCCTTTATGCg    >  1:2053894/1‑36 (MQ=38)
                         agagATAGATGGTCGCGCCTACAGCGCCTTTATTCg    >  1:1957898/1‑36 (MQ=37)
                           agaTAGTTGGCCGCGTCTACAGCGCCTTTATTCGTc  <  1:5212875/36‑1 (MQ=25)
                           agaTAGATGGTCGCGTCTACAGCGCCTTTATTCGTc  <  1:833653/36‑1 (MQ=255)
                              |                                
CGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTC  >  REL606/3762711‑3762773

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG