Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,841 A→T I209F (ATC→TTC)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,8410AT97.4% 295.5 / ‑2.3 76I209F (ATC→TTC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (35/39);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTG  >  REL606/3928771‑3928913
                                                                      |                                                                        
cATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTcgac                                                                     <  2:565747/76‑1 (MQ=255)
 aTGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTcgacg                                                                    >  1:1912712/1‑76 (MQ=255)
   gctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGaa                                                                  >  1:965265/1‑76 (MQ=255)
    ctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAAc                                                                 >  2:881613/1‑76 (MQ=255)
     tgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACg                                                                <  2:1971261/76‑1 (MQ=255)
     tgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACg                                                                <  1:991212/76‑1 (MQ=255)
     tgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACg                                                                >  1:666893/1‑76 (MQ=255)
       ctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTc                                                              >  1:1120197/1‑76 (MQ=255)
       ctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTc                                                              <  1:1920484/76‑1 (MQ=255)
         gcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCg                                                            >  2:78348/1‑76 (MQ=255)
           aGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGga                                                          <  1:1175760/76‑1 (MQ=255)
              aCATTGCTCAGAGCATTGCTTACAACCCCCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGgaaga                                                       <  2:1291113/76‑1 (MQ=255)
              aCATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGgaaga                                                       >  1:1759557/1‑76 (MQ=255)
              aCATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGgaaga                                                       <  2:632559/76‑1 (MQ=255)
              aCATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGgaaga                                                       >  2:701769/1‑76 (MQ=255)
                aTTGCTCAGAGCCTTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGaagaag                                                     <  1:1873285/76‑1 (MQ=255)
                 ttGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGt                                                    <  1:463984/76‑1 (MQ=255)
                 ttGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGt                                                    <  1:955203/76‑1 (MQ=255)
                  tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTa                                                   <  2:1458878/76‑1 (MQ=255)
                   gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTaa                                                  <  2:1072543/76‑1 (MQ=255)
                   ccTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTaa                                                  >  1:2101567/2‑76 (MQ=255)
                    cTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAAc                                                 >  2:1527911/1‑76 (MQ=255)
                     tCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAAcc                                                <  2:84234/76‑1 (MQ=255)
                      cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCg                                               >  1:2077725/1‑76 (MQ=255)
                      cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCg                                               <  2:1566565/76‑1 (MQ=255)
                        gagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCgag                                             <  1:1625432/76‑1 (MQ=255)
                         agCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCgaga                                            <  2:1912712/76‑1 (MQ=255)
                           cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATg                                          >  2:303647/1‑76 (MQ=255)
                            aTTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGc                                         <  1:1777298/76‑1 (MQ=255)
                             ttGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCa                                        >  1:2099362/1‑76 (MQ=255)
                               gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGc                                      >  1:884974/1‑76 (MQ=255)
                               gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGc                                      <  2:965265/76‑1 (MQ=255)
                                cTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCg                                     >  1:453394/1‑76 (MQ=255)
                                cTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCg                                     <  2:1135259/76‑1 (MQ=255)
                                  tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTc                                   <  1:195234/76‑1 (MQ=255)
                                    cAACCCCCCGGATTGTGTGCTGATTGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTg                                 <  1:1270746/76‑1 (MQ=255)
                                    cAACCACCCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTg                                 <  1:445366/76‑1 (MQ=255)
                                      accaccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGt                               <  2:40173/76‑1 (MQ=255)
                                      accaccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGt                               <  2:307749/76‑1 (MQ=255)
                                       ccaccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTa                              >  2:99342/1‑76 (MQ=255)
                                       ccaccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTa                              <  2:2037200/76‑1 (MQ=255)
                                        caccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaa                             <  2:1741422/76‑1 (MQ=255)
                                        caccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaa                             <  2:545153/76‑1 (MQ=255)
                                        caccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaa                             <  2:728676/76‑1 (MQ=255)
                                        caccCGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaa                             <  2:1743744/76‑1 (MQ=255)
                                           ccGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAg                          >  1:197659/1‑76 (MQ=255)
                                           ccGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAg                          <  2:1501004/76‑1 (MQ=255)
                                           ccGGATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAg                          >  2:787043/1‑76 (MQ=255)
                                              gATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTg                       >  2:1548953/1‑76 (MQ=255)
                                              gATTGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTg                       >  2:1021749/1‑76 (MQ=255)
                                                ttGTGTGCTTATGGTTCTGCTGATCCACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGaa                     <  2:1847559/76‑1 (MQ=255)
                                                ttGTGTGCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGaa                     <  1:665906/76‑1 (MQ=255)
                                                 tgtgtgCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAg                    >  1:317296/1‑76 (MQ=255)
                                                 tgtgtgCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAg                    <  1:268554/76‑1 (MQ=255)
                                                 tgtgtgCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAg                    >  1:502833/1‑76 (MQ=255)
                                                 tgtgtgCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAg                    >  1:15773/1‑76 (MQ=255)
                                                     tgCTGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgttgt                >  2:442176/1‑76 (MQ=255)
                                                        tGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGc             <  2:371228/76‑1 (MQ=255)
                                                        tGATGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGCGAAGTTGTTGc             >  1:365284/1‑76 (MQ=255)
                                                          aTGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCtt           >  1:320217/1‑76 (MQ=255)
                                                          aTGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCtn           <  1:815438/76‑2 (MQ=255)
                                                           ttgTTCTTCTGTTCGACGAACGGCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTc          <  2:994786/74‑1 (MQ=255)
                                                           tGGTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTc          >  1:1733213/1‑76 (MQ=255)
                                                            ggTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCt         <  1:239477/76‑1 (MQ=255)
                                                            ggTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCt         <  2:838055/76‑1 (MQ=255)
                                                            ggTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCt         >  2:1754968/1‑76 (MQ=255)
                                                             gTTCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTa        >  2:181597/1‑76 (MQ=255)
                                                              ttCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTAc       >  2:887939/1‑76 (MQ=255)
                                                              ttCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTAc       >  2:966672/1‑76 (MQ=255)
                                                               tCTGCTGTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTAcc      >  1:1494798/1‑76 (MQ=255)
                                                                ctgctgTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCt     <  1:1479695/76‑1 (MQ=255)
                                                                ctgctgTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCt     <  1:1479683/76‑1 (MQ=255)
                                                                  gctgTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCttt   >  2:71635/1‑76 (MQ=255)
                                                                   ctgTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTCCCTTTg  >  2:1949667/1‑76 (MQ=255)
                                                                   ctgTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTg  <  2:676825/76‑1 (MQ=255)
                                                                   ctgTTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTg  >  2:134255/1‑76 (MQ=255)
                                                                      |                                                                        
CATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTG  >  REL606/3928771‑3928913

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