Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,833 T→C V206A (GTT→GCT)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,8330TC98.0% 193.5 / ‑4.1 49V206A (GTT→GCT) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/19);  total (30/19)
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.47e-01

AAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTC  >  REL606/3928758‑3928905
                                                                           |                                                                        
aaGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGt                                                                          >  2:1686747/1‑76 (MQ=255)
   ccGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGctct                                                                       >  1:46537/1‑76 (MQ=255)
      gTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTctgct                                                                    <  1:596285/76‑1 (MQ=255)
      gTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTctgct                                                                    <  1:596297/76‑1 (MQ=255)
       tAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTctgctg                                                                   >  1:632138/1‑76 (MQ=255)
        aaaaCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGa                                                                  >  1:913870/1‑76 (MQ=255)
          aaCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATc                                                                >  1:651202/1‑76 (MQ=255)
           aCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATcg                                                               <  1:1441294/76‑1 (MQ=255)
              aTGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATcgacg                                                            <  1:1735032/76‑1 (MQ=255)
              aTGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATcgacg                                                            >  2:1430147/1‑76 (MQ=255)
                gctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGaa                                                          <  1:1500281/76‑1 (MQ=255)
                gctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGaa                                                          <  2:1693156/76‑1 (MQ=255)
                gctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGaa                                                          <  1:279205/76‑1 (MQ=255)
                 ctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAAc                                                         >  1:1413205/1‑76 (MQ=255)
                 ctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAGTGTGTGCTGATGGCTCTGCTGATCGACGAAc                                                         <  1:1476436/76‑1 (MQ=255)
                  tgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACg                                                        >  2:1528170/1‑76 (MQ=255)
                   gctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGt                                                       >  2:1802546/1‑76 (MQ=255)
                    ctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTc                                                      >  2:633475/1‑76 (MQ=255)
                       cAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCgg                                                   >  2:441691/1‑76 (MQ=255)
                          aaCATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGgaag                                                >  2:1698610/1‑76 (MQ=255)
                          aaCATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGgaag                                                <  2:562922/76‑1 (MQ=255)
                               tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTa                                           <  2:1413340/76‑1 (MQ=255)
                                gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTaa                                          >  1:685072/1‑76 (MQ=255)
                                   cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCg                                       >  1:808938/1‑76 (MQ=255)
                                    agagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCga                                      >  1:1621264/1‑76 (MQ=255)
                                     gagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCgag                                     <  2:318107/76‑1 (MQ=255)
                                     gagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCgag                                     >  2:1158489/1‑76 (MQ=255)
                                        cATTGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATg                                  <  2:1821536/76‑1 (MQ=255)
                                          ttGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCa                                >  1:1825252/1‑76 (MQ=255)
                                          ttGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCa                                >  1:1734511/1‑76 (MQ=255)
                                           tGCTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAg                               >  1:897780/1‑76 (MQ=255)
                                             cTTACAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCg                             >  1:678736/1‑76 (MQ=255)
                                                aCAACCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCt                          <  1:1131924/76‑1 (MQ=255)
                                                  aaCCACCCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTgg                        <  1:1078129/76‑1 (MQ=255)
                                                   accaccCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGt                       <  1:609709/76‑1 (MQ=255)
                                                   accaccCAGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGt                       <  1:609708/76‑1 (MQ=255)
                                                     caccCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaa                     <  1:684341/76‑1 (MQ=255)
                                                     caccCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaa                     >  1:515705/1‑76 (MQ=255)
                                                      accCGGATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGGaaa                    >  2:1012464/1‑76 (MQ=255)
                                                         cgtaTTGTGTGCTGATGACTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAgg                 <  2:431916/73‑1 (MQ=255)
                                                          ggATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGt                <  2:1471344/76‑1 (MQ=255)
                                                          ggATTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGt                <  1:136073/76‑1 (MQ=255)
                                                            aTTGTGTGCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGa              >  2:853267/1‑76 (MQ=255)
                                                                tgtgCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgtt          >  1:446922/1‑76 (MQ=255)
                                                                 gtgCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgttg         >  1:1406841/1‑76 (MQ=255)
                                                                  tgCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgttgt        >  2:978070/1‑76 (MQ=255)
                                                                  tgCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgttgt        >  2:1723289/1‑76 (MQ=255)
                                                                   gCTGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgttgtt       >  2:463379/1‑76 (MQ=255)
                                                                    cTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGGGAAGttgttg      >  1:386551/1‑76 (MQ=255)
                                                                     tGATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGc     >  2:90718/1‑76 (MQ=255)
                                                                      gATGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCt    >  1:149894/1‑76 (MQ=255)
                                                                        tGGCTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTc  >  2:1244189/1‑76 (MQ=255)
                                                                           |                                                                        
AAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTC  >  REL606/3928758‑3928905

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