Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,570,953 G→T D52Y (GAT→TAT)  uspA → universal stress global response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,570,9530GT97.5% 299.0 / 1.5 79D52Y (GAT→TAT) uspAuniversal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (38/39);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTCCACTA  >  REL606/3570880‑3571025
                                                                         |                                                                        
tCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCttatttat                                                                        >  2:1252751/1‑76 (MQ=255)
 cGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTAtg                                                                       <  1:639638/76‑1 (MQ=255)
 cGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTAtg                                                                       >  2:722091/1‑76 (MQ=255)
    cccTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGa                                                                    >  2:1001666/1‑76 (MQ=255)
      cTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAAt                                                                  >  2:1576322/1‑76 (MQ=255)
       tACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATc                                                                 <  1:433434/76‑1 (MQ=255)
        aCAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCt                                                                <  1:1710983/76‑1 (MQ=255)
         cAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTg                                                               <  2:1856629/76‑1 (MQ=255)
         cAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTg                                                               >  2:1018159/1‑76 (MQ=255)
          aaTGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTgg                                                              >  1:257262/1‑76 (MQ=255)
            tGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGt                                                            <  2:1159226/76‑1 (MQ=255)
                aaaGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGata                                                        <  1:1924353/76‑1 (MQ=255)
                 aaGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCATATTTATGTGAATCTGGGTGatat                                                       >  1:1873918/1‑76 (MQ=255)
                   tgttCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGc                                                     <  1:1773648/74‑1 (MQ=255)
                    tttCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCa                                                    <  2:1296166/76‑1 (MQ=255)
                    tttCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATATGAATCTGGGTGATATGCa                                                    <  2:1120205/76‑1 (MQ=255)
                     ttCTCTTATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAg                                                   <  1:1245280/76‑1 (MQ=255)
                     ttCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAg                                                   <  1:1833748/76‑1 (MQ=255)
                        tctGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGaaa                                                <  1:74920/76‑1 (MQ=255)
                         ctGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAAc                                               >  1:1126797/1‑76 (MQ=255)
                           gATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGc                                             >  1:826168/1‑76 (MQ=255)
                            aTCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCa                                            >  1:36507/1‑76 (MQ=255)
                            aTCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCa                                            >  1:1352147/1‑76 (MQ=255)
                            aTCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCa                                            >  1:1409899/1‑76 (MQ=255)
                             tCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCAt                                           >  1:1084180/1‑76 (MQ=255)
                              ccACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCAtc                                          <  2:475211/76‑1 (MQ=255)
                               cACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCAtct                                         >  2:1807198/1‑76 (MQ=255)
                               cACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCAtct                                         >  2:435257/1‑76 (MQ=255)
                                 cGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCa                                       <  1:1812654/76‑1 (MQ=255)
                                  gTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCag                                      <  2:1924019/76‑1 (MQ=255)
                                  gTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCag                                      >  1:675205/1‑76 (MQ=255)
                                   tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCaga                                     <  1:259759/76‑1 (MQ=255)
                                   tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCaga                                     >  1:113593/1‑76 (MQ=255)
                                   tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCaga                                     >  1:1575862/1‑76 (MQ=255)
                                    aGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCagaa                                    >  1:678280/1‑76 (MQ=255)
                                    aGATGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCagaa                                    >  1:88754/1‑76 (MQ=255)
                                       tGTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAagag                                 <  1:123014/76‑1 (MQ=255)
                                        gTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAgaga                                <  1:1001667/76‑1 (MQ=255)
                                        gTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAgaga                                >  1:1123529/1‑76 (MQ=255)
                                        gTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAgaga                                <  1:544953/76‑1 (MQ=255)
                                        gTAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAgaga                                >  1:657082/1‑76 (MQ=255)
                                         tAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGac                               >  1:1226225/1‑76 (MQ=255)
                                         tAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGac                               >  1:513408/1‑76 (MQ=255)
                                         tAAACTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGac                               >  2:1730906/1‑76 (MQ=255)
                                          aaaCTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGaca                              >  2:234239/1‑76 (MQ=255)
                                          aaaCTACTCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGaca                              >  1:1526836/1‑76 (MQ=255)
                                           aaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGacac                             <  2:840183/76‑1 (MQ=255)
                                               actCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACAcatca                         <  1:1121960/76‑1 (MQ=255)
                                               actCTGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACAcatca                         >  1:1492101/1‑76 (MQ=255)
                                                ctctGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACAcatcat                        >  2:568321/1‑76 (MQ=255)
                                                ctctGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACAcatcat                        <  2:973215/76‑1 (MQ=255)
                                                ctctGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACAcatcat                        <  1:143362/76‑1 (MQ=255)
                                                  ctGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATgc                      <  2:1315592/76‑1 (MQ=255)
                                                   tGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATgcg                     <  1:795128/76‑1 (MQ=255)
                                                   tGACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATgcg                     <  2:1462303/76‑1 (MQ=255)
                                                    gACCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATgcgc                    >  1:1753194/1‑76 (MQ=255)
                                                     aCCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCt                   >  1:18148/1‑76 (MQ=255)
                                                     aCCTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCt                   <  2:1443509/76‑1 (MQ=255)
                                                     aCCTATACACCGGGCTTATTTATGTGAATCTGAGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCt                   <  2:1443515/76‑1 (MQ=255)
                                                     aCCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCt                   >  2:551490/1‑76 (MQ=255)
                                                       cTATACACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGa                 <  1:1883399/76‑1 (MQ=255)
                                                         ataCACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGAcc               <  1:1618922/76‑1 (MQ=255)
                                                          taCACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCg              >  1:34559/1‑76 (MQ=255)
                                                          taCACCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCg              <  1:188333/76‑1 (MQ=255)
                                                           acacCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGa             <  2:1918694/76‑1 (MQ=255)
                                                           acacCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGa             >  2:1821468/1‑76 (MQ=255)
                                                           acacCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGa             >  2:1319480/1‑76 (MQ=255)
                                                            cacCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAg            <  1:1573549/76‑1 (MQ=255)
                                                             acCGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGc           >  1:176004/1‑76 (MQ=255)
                                                              ccGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCt          >  1:1325143/1‑76 (MQ=255)
                                                              ccGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCt          <  1:104582/76‑1 (MQ=255)
                                                               cGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCtt         >  1:600095/1‑76 (MQ=255)
                                                               cGGGCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCtt         <  1:907514/76‑1 (MQ=255)
                                                                gggCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCttt        <  1:1756109/76‑1 (MQ=255)
                                                                 ggCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTc       >  1:1555611/1‑76 (MQ=255)
                                                                  gCTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTcc      <  1:372900/76‑1 (MQ=255)
                                                                   cTTATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTCCa     >  1:1862122/1‑76 (MQ=255)
                                                                    ttatttatGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTCCAc    >  1:1268856/1‑76 (MQ=255)
                                                                    ttatttatGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTCCAc    >  2:957226/1‑76 (MQ=255)
                                                                     tATTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTCCACt   <  2:1422653/76‑1 (MQ=255)
                                                                      aTTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTCCACTa  <  1:1950107/76‑1 (MQ=255)
                                                                      aTTTATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACAGCATGCGCTGACCGAGCTTTCCACtg  <  1:547747/76‑2 (MQ=255)
                                                                         |                                                                        
TCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCATGCGCTGACCGAGCTTTCCACTA  >  REL606/3570880‑3571025

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