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 ref new freq score (cons/poly) reads annotation genes product
*REL6063,762,7410AT98.2% 207.7 / NA 55K662I (AAA→ATA) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (31/23);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG  >  REL606/3762668‑3762811
                                                                         |                                                                      
ggTGCGCGGGCAAACCTGACGGAGGAAATTAACACCGCGACTTTGAATATTCAAAGTTTGAATACGGAAGAGATAg                                                                      <  1:1295264/76‑1 (MQ=255)
 gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa                                                                     <  1:648597/76‑1 (MQ=255)
  tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat                                                                    <  1:398833/76‑1 (MQ=255)
  ggcgcTTGCAAACCTGACGGCGGCAATTAACCCCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat                                                                    <  1:398832/75‑1 (MQ=255)
     gcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGt                                                                 <  2:414558/76‑1 (MQ=255)
      cTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTc                                                                >  1:768990/1‑76 (MQ=255)
         gCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg                                                             <  2:1499464/76‑1 (MQ=255)
          cAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                                                            >  1:1241850/1‑76 (MQ=255)
          cAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt                                                            >  1:1241827/1‑76 (MQ=255)
            aaCCTGACGGCGGCTATTAACACCGCGACTTNGAATATTCGAAGTTNGAATANGGAAGAGATAGATGGTCGCGTCt                                                          <  2:203502/76‑1 (MQ=255)
             aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                                                         <  1:1338242/76‑1 (MQ=255)
             aCCTGACGGCGACAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                                                         <  1:568918/76‑1 (MQ=255)
              ccTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc                                                        <  1:557041/76‑1 (MQ=255)
               cTGACGGCTGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                                                       <  1:1370494/76‑1 (MQ=255)
                tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg                                                      <  1:355990/76‑1 (MQ=255)
                 gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgc                                                     >  1:458763/1‑76 (MQ=255)
                  aCGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg                                                    <  2:744898/76‑1 (MQ=255)
                      cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt                                                <  2:868450/76‑1 (MQ=255)
                      cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt                                                <  1:1124265/76‑1 (MQ=255)
                        gcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa                                              >  1:519671/1‑76 (MQ=255)
                         cAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt                                             >  2:1164935/1‑76 (MQ=255)
                          aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt                                            <  1:699829/76‑1 (MQ=255)
                           aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc                                           >  2:960921/1‑76 (MQ=255)
                            ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                          <  2:439135/76‑1 (MQ=255)
                            ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg                                          >  1:699042/1‑76 (MQ=255)
                                  ccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAc                                    >  2:118046/1‑76 (MQ=255)
                                       aCTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTc                               >  1:373485/1‑76 (MQ=255)
                                        cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg                              >  2:1488639/1‑76 (MQ=255)
                                        cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg                              >  2:1070762/1‑76 (MQ=255)
                                         ttCGAATATTCAAAGTTTGNATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTCTTCGTCTGACCGCTCGt                             >  1:1432706/1‑76 (MQ=255)
                                          tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg                            <  2:97713/76‑1 (MQ=255)
                                           cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGGGa                           >  1:1008135/1‑76 (MQ=255)
                                           cGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGATCGTGa                           >  1:1390686/1‑76 (MQ=255)
                                               tatTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgt                       >  1:544469/1‑76 (MQ=255)
                                                nttCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGANCGCTCGTGACCgtg                      >  2:219758/2‑76 (MQ=255)
                                                atTCAAAGTTTGAATACGGAAGAGATAGTTGGTCCCATCTACAGCGCNTTAATTCGTCTGACCGCTCGTGACCgtg                      >  1:142200/1‑76 (MQ=255)
                                                atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                      >  1:574445/1‑76 (MQ=255)
                                                 ttCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtgt                     >  1:524912/1‑76 (MQ=255)
                                                  tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGtgtg                    <  1:559312/76‑1 (MQ=255)
                                                    aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCa                  >  1:495958/1‑76 (MQ=255)
                                                      aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATc                >  1:285751/1‑76 (MQ=255)
                                                         ttGAATACGGAAGAGATAGATGGTCGCGTCTACNGCGCCTTTATTCGTCTGACCGCNCGTGACCGAGNGCATCTgg             >  2:1278581/1‑76 (MQ=255)
                                                         ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg             >  1:835910/1‑76 (MQ=255)
                                                          tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGc            <  1:983826/76‑1 (MQ=255)
                                                           gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg           <  1:735331/76‑1 (MQ=255)
                                                           gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg           <  2:771604/76‑1 (MQ=255)
                                                           gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg           <  1:594855/76‑1 (MQ=255)
                                                            aaTACGGAAGAGATAGATGGTCGCGTCTACTGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa          >  2:66729/1‑76 (MQ=255)
                                                            aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa          >  2:581060/1‑76 (MQ=255)
                                                            aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa          >  2:522191/1‑76 (MQ=255)
                                                            aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa          >  2:289411/1‑76 (MQ=255)
                                                             aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGaa         >  2:811593/1‑76 (MQ=255)
                                                              tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat        <  1:127731/76‑1 (MQ=255)
                                                              tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat        >  2:1338144/1‑76 (MQ=255)
                                                                  gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa    <  2:731898/76‑1 (MQ=255)
                                                                   aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACGNCTNGTGACCGTGTGCATNTGGCGAATATCAt   >  2:1384508/1‑76 (MQ=255)
                                                                   aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt   >  2:524356/1‑76 (MQ=255)
                                                                    agagATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATg  >  2:849387/1‑76 (MQ=255)
                                                                         |                                                                      
GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG  >  REL606/3762668‑3762811

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: