Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,570,803 Δ1 bp coding (4/435 nt) uspA → universal stress global response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,570,8020G.98.4% 314.0 / ‑0.9 64coding (3/435 nt)uspAuniversal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (30/33);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

GACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTAT  >  REL606/3570727‑3570876
                                                                           |                                                                          
gACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATg                                                                            <  2:1090159/76‑1 (MQ=255)
gACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATg                                                                            >  1:1700476/1‑76 (MQ=255)
      ggTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTtata                                                                     <  2:2103268/76‑1 (MQ=255)
      ggTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTtata                                                                     <  1:37610/76‑1 (MQ=255)
        tGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATaaa                                                                   <  1:1473831/76‑1 (MQ=255)
         gACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAac                                                                  <  2:465779/76‑1 (MQ=255)
          aCGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAaca                                                                 >  1:486949/1‑76 (MQ=255)
           cGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGGCTTATAAaca                                                                 <  1:1810562/76‑1 (MQ=255)
           cGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAacac                                                                >  2:54446/1‑76 (MQ=255)
            gTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAAcaca                                                               >  1:1489847/1‑76 (MQ=255)
            gTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAAcaca                                                               <  1:2094654/76‑1 (MQ=255)
             tAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACAt                                                              <  2:1156414/76‑1 (MQ=255)
              aaCGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACAtt                                                             <  2:59662/76‑1 (MQ=255)
              aaCGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACAtt                                                             >  2:1886364/1‑76 (MQ=255)
                cGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATtct                                                           <  2:377663/76‑1 (MQ=255)
                cGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATtct                                                           <  2:278037/76‑1 (MQ=255)
                         aaaCGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGt                                                  >  1:1016913/1‑76 (MQ=255)
                          aaCGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTc                                                 <  2:588379/76‑1 (MQ=255)
                             gCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAt                                              <  1:1980108/76‑1 (MQ=255)
                               tAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtct                                            >  1:1417090/1‑76 (MQ=255)
                               tAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtct                                            >  1:1350772/1‑76 (MQ=255)
                                aGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctc                                           >  1:121590/1‑76 (MQ=255)
                                aGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctc                                           <  1:61442/76‑1 (MQ=255)
                                aGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctc                                           >  1:1773911/1‑76 (MQ=255)
                                aGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctc                                           >  1:1364997/1‑76 (MQ=255)
                                 gCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctct                                          <  1:604682/76‑1 (MQ=255)
                                 gCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctct                                          <  2:195386/76‑1 (MQ=255)
                                  cTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctctc                                         >  1:687447/1‑76 (MQ=255)
                                  cTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGAtctctc                                         >  2:1719100/1‑76 (MQ=255)
                                    ggCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTccc                                       <  2:549796/76‑1 (MQ=255)
                                     gCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTcccc                                      <  2:1530059/76‑1 (MQ=255)
                                     gCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTcccc                                      >  1:1881926/1‑76 (MQ=255)
                                      ccAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCg                                     <  1:76317/76‑1 (MQ=255)
                                      ccAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCg                                     <  1:927025/76‑1 (MQ=255)
                                      ccAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCg                                     <  2:2050173/76‑1 (MQ=255)
                                       cAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGTGGTCGCTCTCTCCCCgg                                    >  2:1425292/1‑76 (MQ=255)
                                       cAGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCgg                                    >  2:1208063/1‑76 (MQ=255)
                                        aGTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGa                                   <  2:1700476/76‑1 (MQ=255)
                                         gTCATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGaa                                  <  1:911642/76‑1 (MQ=255)
                                           cATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAg                                <  2:1719499/76‑1 (MQ=255)
                                           cATCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAg                                <  2:548977/76‑1 (MQ=255)
                                            aTCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGc                               >  2:810187/1‑76 (MQ=255)
                                            aTCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGATCGATCTCTCCCCGGAAAGc                               <  1:627263/76‑1 (MQ=255)
                                             tCGACAACTTTATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCa                              <  1:1675563/76‑1 (MQ=255)
                                                     tttATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTg                      >  1:635490/1‑76 (MQ=255)
                                                      ttATGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTgg                     <  2:993290/76‑1 (MQ=255)
                                                        aTGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTa                   >  1:745122/1‑76 (MQ=255)
                                                        aTGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTa                   >  1:476235/1‑76 (MQ=255)
                                                         tGGAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTag                  >  2:1143586/1‑76 (MQ=255)
                                                          ggAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCGCCGGAAAGCAAAGTTCTGGTaga                 >  2:1873768/1‑76 (MQ=255)
                                                          ggAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTaga                 >  1:2069477/1‑76 (MQ=255)
                                                          ggAAGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTaga                 <  2:189329/76‑1 (MQ=255)
                                                            aaGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAgaga               <  1:1282562/76‑1 (MQ=255)
                                                             aGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGaa              <  1:530302/76‑1 (MQ=255)
                                                             aGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGaa              >  1:380534/1‑76 (MQ=255)
                                                             aGGAGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGaa              >  2:267730/1‑76 (MQ=255)
                                                                aGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGc           <  1:881903/76‑1 (MQ=255)
                                                                aGTAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGc           <  1:1246933/76‑1 (MQ=255)
                                                                  tAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAg         <  2:629215/76‑1 (MQ=255)
                                                                  tAACACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAg         >  2:1654510/1‑76 (MQ=255)
                                                                   aaCACTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGt        >  1:1823129/1‑76 (MQ=255)
                                                                    acacTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGtc       >  2:1691574/1‑76 (MQ=255)
                                                                     cacTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGtct      <  2:1440876/76‑1 (MQ=255)
                                                                     cacTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGtct      >  2:2050690/1‑76 (MQ=255)
                                                                     cacTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGtct      >  2:801032/1‑76 (MQ=255)
                                                                      acTATGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGtctc     >  2:382375/1‑76 (MQ=255)
                                                                          tGGCTTATAAACACACTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTAt  <  2:1469155/76‑1 (MQ=255)
                                                                           |                                                                          
GACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGGCTTATAAACACATTCTCATCGCGGTCGATCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTAT  >  REL606/3570727‑3570876

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