Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,259 A→T I15F (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,2590AT95.1% 231.5 / 3.2 61I15F (ATC→TTC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (23/35);  total (24/37)
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.05e-01

TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC  >  REL606/3928188‑3928327
                                                                       |                                                                    
ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAc                                                                  >  1:1054711/1‑76 (MQ=255)
 tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAct                                                                 >  2:308390/1‑76 (MQ=255)
 tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAct                                                                 <  1:512815/76‑1 (MQ=255)
    cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACtctc                                                              <  2:1330166/76‑1 (MQ=255)
    cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACtctc                                                              <  2:1615375/76‑1 (MQ=255)
      cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCgg                                                            >  1:804928/1‑76 (MQ=255)
       aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGc                                                           <  2:493606/76‑1 (MQ=255)
        aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCg                                                          <  2:1239620/76‑1 (MQ=255)
          tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGaa                                                        <  1:354244/76‑1 (MQ=255)
            tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaaaa                                                      >  1:803855/1‑76 (MQ=255)
             aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAat                                                     <  1:227678/76‑1 (MQ=255)
              aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAata                                                    >  2:1365004/1‑76 (MQ=255)
              aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAata                                                    <  1:974870/76‑1 (MQ=255)
               taaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAatat                                                   >  2:190630/2‑76 (MQ=255)
                aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATg                                                  >  1:386779/1‑76 (MQ=255)
                aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATg                                                  >  1:386798/1‑76 (MQ=255)
                   tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATgggg                                               <  1:950473/76‑1 (MQ=255)
                    cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGc                                              >  2:876990/1‑76 (MQ=255)
                    cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGc                                              >  2:200546/1‑76 (MQ=255)
                    cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGATTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGc                                              >  2:527725/1‑76 (MQ=255)
                     acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCt                                             <  1:668242/76‑1 (MQ=255)
                      caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTg                                            <  1:1239038/76‑1 (MQ=255)
                        ccacTATGAATCTTACCAAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGa                                          <  1:360440/76‑1 (MQ=255)
                           cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGaaaa                                       <  1:97172/76‑1 (MQ=255)
                                  tCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCt                                <  2:1647518/76‑1 (MQ=255)
                                     tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                             <  1:1008864/76‑1 (MQ=255)
                                      aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                            <  2:139196/76‑1 (MQ=255)
                                       ccGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                           <  2:1100724/76‑1 (MQ=255)
                                        cGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATg                          <  2:804927/76‑1 (MQ=255)
                                         gAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc                         <  1:1705251/76‑1 (MQ=255)
                                          aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg                        <  1:785439/76‑1 (MQ=255)
                                            ttAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                      <  2:758289/76‑1 (MQ=255)
                                            ttAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                      <  1:18902/76‑1 (MQ=255)
                                            ttAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                      <  1:308391/76‑1 (MQ=255)
                                             tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGtga                     <  2:1575077/76‑3 (MQ=255)
                                             tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                     <  1:785653/76‑1 (MQ=255)
                                             tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                     <  2:1614569/76‑1 (MQ=255)
                                             tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                     >  1:807673/1‑76 (MQ=255)
                                              aaaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                    >  1:1263558/1‑76 (MQ=255)
                                               aaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                   <  1:335157/76‑1 (MQ=255)
                                               aaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                   <  2:511921/76‑1 (MQ=255)
                                                aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                  <  2:368875/76‑1 (MQ=255)
                                                aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                  >  2:1271410/1‑76 (MQ=255)
                                                aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                  <  2:368861/76‑1 (MQ=255)
                                                aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                  >  1:1032411/1‑76 (MQ=255)
                                                 gAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                 >  1:475415/1‑76 (MQ=255)
                                                  aaTACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                             <  1:319212/63‑1 (MQ=255)
                                                  aaTACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                             >  2:319211/1‑63 (MQ=255)
                                                  aaTACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg                <  2:930296/76‑1 (MQ=255)
                                                    tACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc              <  1:567306/76‑1 (MQ=255)
                                                    tACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc              >  2:681456/1‑76 (MQ=255)
                                                      cGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat            <  1:1374018/76‑1 (MQ=255)
                                                         cGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat         <  1:1088032/76‑1 (MQ=255)
                                                          ggTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt        >  1:1012908/1‑76 (MQ=255)
                                                            tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt      <  2:1576687/76‑1 (MQ=255)
                                                            tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt      <  2:1576682/76‑1 (MQ=255)
                                                            tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt      <  2:1038447/76‑1 (MQ=255)
                                                            tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt      >  1:329900/1‑76 (MQ=255)
                                                            nttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt      >  1:173989/2‑76 (MQ=255)
                                                             ttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttttt     >  1:66523/1‑76 (MQ=255)
                                                               cTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGc   >  2:957769/1‑76 (MQ=255)
                                                                tGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGcc  >  1:266267/1‑76 (MQ=255)
                                                                       |                                                                    
TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC  >  REL606/3928188‑3928327

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