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,2600TA94.5% 206.4 / 5.2 55I15N (ATC→AAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (27/25);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-02

TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928187‑3928332
                                                                         |                                                                        
tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACa                                                                        <  2:1501134/76‑1 (MQ=255)
 ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       >  1:3503/1‑76 (MQ=255)
   aCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActc                                                                     >  1:1443817/1‑76 (MQ=255)
    ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct                                                                    <  2:1310859/76‑1 (MQ=255)
    ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct                                                                    <  2:685654/76‑1 (MQ=255)
     cccAAGTTTAAGAACTCGCACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   <  1:1731333/76‑1 (MQ=255)
      ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg                                                                  >  2:466815/1‑76 (MQ=255)
      ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg                                                                  <  2:473285/76‑1 (MQ=255)
       cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCgg                                                                 <  1:438056/76‑1 (MQ=255)
         aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               >  1:1770794/1‑76 (MQ=255)
           tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             <  2:1772321/76‑1 (MQ=255)
           tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             >  1:1549524/1‑76 (MQ=255)
            ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa                                                            >  1:60219/1‑76 (MQ=255)
             tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           >  2:1874697/1‑76 (MQ=255)
                gAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAatat                                                        <  2:723258/76‑1 (MQ=255)
                 aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATg                                                       >  1:809123/1‑76 (MQ=255)
                 aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATg                                                       <  1:630870/76‑1 (MQ=255)
                   cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAATATgggg                                                    <  2:872592/76‑1 (MQ=255)
                    tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg                                                    >  2:1955134/1‑76 (MQ=255)
                      acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt                                                  <  2:1807902/76‑1 (MQ=255)
                       caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                 <  1:375946/76‑1 (MQ=255)
                        accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                <  2:1344098/76‑1 (MQ=255)
                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa                                              <  1:205372/76‑1 (MQ=255)
                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa                                              <  1:77245/76‑1 (MQ=255)
                           acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa                                             <  1:1933328/76‑1 (MQ=255)
                            cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaaa                                            >  2:1306411/1‑76 (MQ=255)
                             tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAc                                           <  2:1193492/76‑1 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  2:1368485/1‑76 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  1:16274/1‑76 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  2:83070/1‑76 (MQ=255)
                                  aTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGc                                      >  1:681383/1‑76 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    <  2:1907849/76‑1 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    >  1:486740/1‑76 (MQ=255)
                                     ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                                   <  2:1487282/76‑1 (MQ=255)
                                     ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                                   >  2:1987332/1‑76 (MQ=255)
                                      tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                                  <  1:1488833/76‑1 (MQ=255)
                                        ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  2:45024/76‑1 (MQ=255)
                                           aaTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg                             >  1:328154/1‑76 (MQ=255)
                                             ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           >  1:973874/1‑76 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCCCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  1:456057/76‑1 (MQ=255)
                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         >  2:1665566/1‑76 (MQ=255)
                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         >  2:1944470/1‑76 (MQ=255)
                                                aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                        <  1:261027/76‑1 (MQ=255)
                                                  gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      >  1:991884/1‑76 (MQ=255)
                                                  gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      <  1:328693/76‑1 (MQ=255)
                                                   aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg                     <  1:304032/76‑1 (MQ=255)
                                                       cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTGGTATGCGTAAGCAGGACat                 <  1:1245537/76‑1 (MQ=255)
                                                       cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 <  1:1951183/76‑1 (MQ=255)
                                                         ccGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta               <  2:812945/76‑1 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTCAGCAGGCCattat              >  2:719595/1‑76 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              >  2:372847/1‑76 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              >  2:719592/1‑76 (MQ=255)
                                                           ngTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt             >  1:345390/2‑76 (MQ=255)
                                                            gTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            <  1:1222418/76‑1 (MQ=255)
                                                                  gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa      >  1:558913/1‑76 (MQ=255)
                                                                   aGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAt     >  1:7230/1‑76 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  1:970654/1‑76 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  2:1125389/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: