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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,2600TA97.2% 128.4 / ‑3.0 36I15N (ATC→AAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (18/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928188‑3928332
                                                                        |                                                                        
ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       >  2:1140584/1‑76 (MQ=255)
 tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                      >  2:975496/1‑76 (MQ=255)
    cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGNACACtctc                                                                   <  1:1210568/76‑1 (MQ=255)
      cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCNGGTTGCTCAGCTGAACACTCGCgg                                                                 <  1:106282/76‑1 (MQ=255)
      cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCgg                                                                 <  1:888588/76‑1 (MQ=255)
       aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGc                                                                >  2:633333/1‑76 (MQ=255)
        aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACGCTCGGCg                                                               <  2:487396/76‑1 (MQ=255)
              aGAACTCATACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTTAACACTCTCGGCGAAAata                                                         >  2:18446/1‑76 (MQ=255)
              aGAACTCACACCACTATGAATCTTACCGACTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         <  2:34249/76‑1 (MQ=255)
                  cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATggg                                                     >  1:559430/1‑76 (MQ=255)
                   tCACACCACTATGANTCTTATCGAACTAAAGAATCCGCCAGTTCCTGAGCAGAACACTCTCGGCGAATATATgggg                                                    <  2:28161/76‑1 (MQ=255)
                     acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt                                                  <  1:985922/76‑1 (MQ=255)
                       accaccCTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                >  1:1188431/1‑76 (MQ=255)
                            tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGANAAc                                           <  1:1140585/76‑1 (MQ=255)
                            tATGAATCTTAACGAATTAAAGAATACGCCGGTCTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAc                                           >  1:1163031/1‑76 (MQ=255)
                             aTGAATCTTTCCGAATTAAAGAATACGCCGGNGTCTAAGGTGAACACTCTCGGCGAAAATATGCNGCTGGAAAAcc                                          >  1:149979/1‑76 (MQ=255)
                               gAATCTTACCGAATCAAAGAATACGCCGGTGTCTGAGCTCAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTg                                        <  1:18446/76‑1 (MQ=255)
                                  tCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCt                                     <  2:761946/76‑1 (MQ=255)
                                   cTTACCGAATTAGAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    <  1:1017059/76‑1 (MQ=255)
                                         gAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc                              <  1:855734/76‑1 (MQ=255)
                                             tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          >  1:572297/1‑76 (MQ=255)
                                              aaaGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         <  1:344258/76‑1 (MQ=255)
                                               aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                        <  2:567187/76‑1 (MQ=255)
                                                 gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      <  1:557982/76‑1 (MQ=255)
                                                        ccGGTTCCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta               >  1:1185167/1‑76 (MQ=255)
                                                         cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              >  2:516268/1‑76 (MQ=255)
                                                          ggTTTCTGAGCTGAACACTCTCGGCGAAAATGTGNGGCTTGANCANCTNGCTCTCATGCGCAAGCAGGACattatt             >  1:1105333/1‑76 (MQ=255)
                                                           gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            <  2:501440/76‑1 (MQ=255)
                                                           gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            <  2:147370/76‑1 (MQ=255)
                                                           gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            >  2:574898/1‑76 (MQ=255)
                                                           gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            >  1:501140/1‑76 (MQ=255)
                                                           gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            >  2:911905/1‑76 (MQ=255)
                                                            tttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt           >  2:939817/1‑76 (MQ=255)
                                                                 gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa      >  2:930786/1‑76 (MQ=255)
                                                                   gCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATc    <  2:1203839/76‑1 (MQ=255)
                                                                    cTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATcc   >  1:1024082/1‑76 (MQ=255)
                                                                     tGAACACTCTCGTCGAAAAAATGGGGCTGAAAATCCTGACTCGACTGCGTAAGCAGGACATTATTTTATCCATCCt  >  1:1222655/1‑76 (MQ=255)
                                                                        |                                                                        
TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928188‑3928332

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