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,2600TA96.7% 240.2 / ‑0.3 60I15N (ATC→AAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (30/28);  total (31/29)
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.72e-01

TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928187‑3928332
                                                                         |                                                                        
tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACa                                                                        <  1:1051558/76‑1 (MQ=255)
tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACa                                                                        <  1:83626/76‑1 (MQ=255)
 ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       >  2:1746823/1‑76 (MQ=255)
     cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   <  2:75361/76‑1 (MQ=255)
     cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   <  2:687487/76‑1 (MQ=255)
      ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCg                                                                  >  1:271159/1‑76 (MQ=255)
       cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCgg                                                                 <  2:1865130/76‑1 (MQ=255)
            ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa                                                            <  2:1510430/76‑1 (MQ=255)
             tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           <  2:1921381/76‑1 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          >  1:630507/1‑76 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          <  1:401865/76‑1 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          >  2:1098434/1‑76 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          >  1:1351272/1‑76 (MQ=255)
               aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         <  2:1415488/76‑1 (MQ=255)
               aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         >  1:669966/1‑76 (MQ=255)
               aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         <  2:669966/76‑1 (MQ=255)
                gAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAatat                                                        >  2:1389343/1‑76 (MQ=255)
                gAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAatat                                                        >  2:89311/1‑76 (MQ=255)
                  aCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgg                                                      <  2:322575/76‑1 (MQ=255)
                    tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg                                                    >  2:73825/1‑76 (MQ=255)
                        accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                <  1:1211262/76‑1 (MQ=255)
                         ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGa                                               >  2:620474/1‑76 (MQ=255)
                         ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGa                                               >  1:1941797/1‑76 (MQ=255)
                            cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGaaaa                                            <  1:77082/76‑1 (MQ=255)
                            cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaaa                                            <  2:1152811/76‑1 (MQ=255)
                             tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAc                                           >  1:1318386/1‑76 (MQ=255)
                             tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAc                                           >  1:140930/1‑76 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  1:545100/1‑76 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  2:1438299/1‑76 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATAAGGGGCTGGAAAAcc                                          >  1:1152943/1‑76 (MQ=255)
                               tGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCt                                         <  1:1341839/76‑1 (MQ=255)
                                 aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                       >  2:675030/1‑76 (MQ=255)
                                 aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                       >  2:1131452/1‑76 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    >  1:529046/1‑76 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    <  2:553242/76‑1 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    >  2:1860724/1‑76 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    >  2:7167/1‑76 (MQ=255)
                                     ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                                   <  2:625004/76‑1 (MQ=255)
                                     ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                                   >  2:1238026/1‑76 (MQ=255)
                                       aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                                 <  2:1274162/76‑1 (MQ=255)
                                       aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                                 >  1:829673/1‑76 (MQ=255)
                                        ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  2:1849497/76‑1 (MQ=255)
                                        ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  1:150307/76‑1 (MQ=255)
                                          gAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc                              <  2:1791450/76‑1 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          >  1:1998509/1‑76 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  1:1751668/76‑1 (MQ=255)
                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         <  1:309287/76‑1 (MQ=255)
                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         <  2:1938378/76‑1 (MQ=255)
                                                 aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                       <  2:1996333/76‑1 (MQ=255)
                                                 aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                       <  1:1948049/76‑1 (MQ=255)
                                                 aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                       >  2:1768976/1‑76 (MQ=255)
                                                 aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                       >  1:1471086/1‑76 (MQ=255)
                                                   aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg                     >  1:1440973/1‑76 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              <  2:396452/76‑1 (MQ=255)
                                                            gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            <  2:957268/76‑1 (MQ=255)
                                                              ttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttttt          >  2:327688/1‑76 (MQ=255)
                                                                cTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGc        >  1:292570/1‑76 (MQ=255)
                                                                 tGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGcc       >  1:1103920/1‑76 (MQ=255)
                                                                   aGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAt     <  1:675030/76‑1 (MQ=255)
                                                                    gCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATc    >  2:83684/1‑76 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  2:1736073/1‑76 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  1:1225085/1‑76 (MQ=255)
                                                                         |                                                                        
TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928187‑3928332

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