Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,944 C→T A243V (GCG→GTG)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,9440CT98.2% 200.4 / ‑3.6 55A243V (GCG→GTG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (25/29);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

ATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCAT  >  REL606/3928871‑3929016
                                                                         |                                                                        
aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTgg                                                                        >  2:988066/1‑76 (MQ=255)
aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTgg                                                                        <  1:640681/76‑1 (MQ=255)
 tGCAGCGTCTGGTAAAAGGAGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGCTGTGGa                                                                       >  1:760788/1‑76 (MQ=255)
   cAGCGTCTGGTAANAGGTGAAGTTGTTGCTTCTACGTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGaaa                                                                     <  2:1305100/76‑1 (MQ=255)
      cGTCTGGTAAAAGGTGAAGTTGTTGTTTCTAGCTTTGACTNNCCCGCATCTCNNCGCGTTGAGGTTGTAGAAATgg                                                                  >  1:1283499/1‑76 (MQ=255)
       gTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGt                                                                 >  1:315067/1‑76 (MQ=255)
        tCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTg                                                                <  2:1309644/76‑1 (MQ=255)
        tCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTg                                                                <  2:1061557/76‑1 (MQ=255)
           ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATc                                                             <  2:676184/76‑1 (MQ=255)
           ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATc                                                             <  2:585289/76‑1 (MQ=255)
           ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTAGCCACGTTCAGGTTGTGGAAACGGTGATc                                                             >  1:72674/1‑76 (MQ=255)
            gTAAAATGTGAAATTGTTGCTACTGCCTATGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCg                                                            <  2:242632/76‑1 (MQ=255)
             tAAAAGGTGAAGTAGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCga                                                           >  1:101411/1‑76 (MQ=255)
              aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCgag                                                          >  1:1241360/1‑76 (MQ=255)
               aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCgaga                                                         <  2:288625/76‑1 (MQ=255)
                aaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGaa                                                        >  2:1230602/1‑76 (MQ=255)
                  ggTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGCTCAGGTTGTGGAAATGGTGATCGAGAAgg                                                      >  1:311099/1‑76 (MQ=255)
                   gTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGc                                                     <  1:539290/76‑1 (MQ=255)
                   gTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGc                                                     >  1:797522/1‑76 (MQ=255)
                     gAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCa                                                   <  2:1378822/76‑1 (MQ=255)
                      aaGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCaa                                                  <  1:623248/76‑1 (MQ=255)
                       aGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCaaa                                                 >  1:642109/1‑76 (MQ=255)
                         ttgttgCTTCTACCTTTGACGACCCCGCATCTCGCCACGTTCAGGTTGTGGAAAAGGTGATAGAGAAGGCCAAACg                                               >  2:213117/1‑76 (MQ=255)
                         ttgttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACg                                               <  1:216753/76‑1 (MQ=255)
                          tgttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGc                                              >  1:410700/1‑76 (MQ=255)
                           gttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGcc                                             >  1:581230/1‑76 (MQ=255)
                           gttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGcc                                             >  1:687544/1‑76 (MQ=255)
                            ttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCt                                            >  1:862432/1‑76 (MQ=255)
                             tgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTg                                           <  2:1135845/76‑1 (MQ=255)
                                ttCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt                                        >  1:462412/1‑76 (MQ=255)
                                ttCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt                                        <  2:942920/76‑1 (MQ=255)
                                 tCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGt                                         <  2:722939/74‑1 (MQ=255)
                                 tCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGt                                         >  1:722939/1‑74 (MQ=255)
                                      cTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcac                                  <  2:672780/76‑1 (MQ=255)
                                        ttGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaa                                <  2:1074890/76‑1 (MQ=255)
                                        ttGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaa                                <  1:213117/76‑1 (MQ=255)
                                        ttGACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaa                                <  2:872029/76‑1 (MQ=255)
                                         tGACGAACCCGCATCTCGCNACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaag                               >  1:1334997/1‑76 (MQ=255)
                                          gACGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                              <  1:978387/76‑1 (MQ=255)
                                          gACGAACCCGCATCTCGCCACGTACAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                              <  1:788766/76‑1 (MQ=255)
                                           aCGAACCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                             <  1:566495/76‑1 (MQ=255)
                                              aaCCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaaga                          <  1:1197196/76‑1 (MQ=255)
                                               aCCCGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGAc                         >  2:647489/1‑76 (MQ=255)
                                                  cGCATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGtt                      <  2:427633/76‑1 (MQ=255)
                                                    cATCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTAt                    <  1:240766/76‑1 (MQ=255)
                                                     aTCTCGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGATTGAGCACACGAAAGACGTTATc                   >  1:1229869/1‑76 (MQ=255)
                                                      tctcGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCa                  <  1:113646/76‑1 (MQ=255)
                                                      tctcGCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCa                  <  2:601509/76‑1 (MQ=255)
                                                          gCCACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCt              <  2:1378542/76‑1 (MQ=255)
                                                          cccACGTTCAGGTTGTGGAAATGGTGATGGAGAATGCCAAACGCCTGNTTAAGCACAAGAAAGACGTTATCATTCt              >  1:124079/2‑76 (MQ=255)
                                                            cACGTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGc            >  2:619913/1‑76 (MQ=255)
                                                               gTTCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCg         <  2:6599/76‑1 (MQ=255)
                                                                 tCAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGAc       >  1:1051845/1‑76 (MQ=255)
                                                                  cAGGTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACt      >  1:1081383/1‑76 (MQ=255)
                                                                     gTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCa   <  2:760788/76‑1 (MQ=255)
                                                                     gTTGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCa   >  2:58992/1‑76 (MQ=255)
                                                                      ttGTGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCAt  >  1:583505/1‑76 (MQ=255)
                                                                         |                                                                        
ATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCAT  >  REL606/3928871‑3929016

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: