Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,595 G→A A127T (GCC→ACC)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,5950GA98.9% 1041.2 / NA 277A127T (GCC→ACC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (137/136);  minor base C (0/3);  total (137/140)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 279 total aligned reads displayed.

CCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTA  >  REL606/3928522‑3928667
                                                                         |                                                                        
cccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACAcc                                                                        >  1:4313616/1‑76 (MQ=255)
  cgccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCg                                                                      <  1:4581128/76‑1 (MQ=255)
   gccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGc                                                                     <  2:4308352/76‑1 (MQ=255)
       aaaGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGcaaca                                                                 >  2:4623266/1‑76 (MQ=255)
        aagaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGcaacaa                                                                >  2:5444508/1‑76 (MQ=255)
        aagaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGcaacaa                                                                <  1:4897166/76‑1 (MQ=255)
        aagaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGcaacaa                                                                <  1:4069232/76‑1 (MQ=255)
         agaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACaaa                                                               >  1:4198169/1‑76 (MQ=255)
         agaagGAGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACaaa                                                               >  2:5182765/1‑76 (MQ=255)
          gaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACaaaa                                                              >  2:5305162/1‑76 (MQ=255)
          gaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACaaaa                                                              <  1:3039217/76‑1 (MQ=255)
           aagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAAt                                                             <  1:1510587/76‑1 (MQ=255)
              gTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCct                                                          <  2:2910637/76‑1 (MQ=255)
              gTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCct                                                          <  2:3893499/76‑1 (MQ=255)
                gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCctct                                                        <  2:79198/76‑1 (MQ=255)
                gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCctct                                                        >  2:6579642/1‑76 (MQ=255)
                gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCctct                                                        >  1:747138/1‑76 (MQ=255)
                   cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTg                                                     <  1:1185060/76‑1 (MQ=255)
                   cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTg                                                     <  2:1538853/76‑1 (MQ=255)
                   cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTg                                                     <  2:1261368/76‑1 (MQ=255)
                   cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTg                                                     <  2:892259/76‑1 (MQ=255)
                    gCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTga                                                    <  1:982801/76‑1 (MQ=255)
                     cTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTgag                                                   <  1:5617501/76‑1 (MQ=255)
                     cTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTgag                                                   <  1:1577293/76‑1 (MQ=255)
                      tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAGAAACCTGAAAACACCCGCAACAAAATCCTCTTTgaga                                                  <  1:1053599/76‑1 (MQ=255)
                      tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTgaga                                                  >  2:3342160/1‑76 (MQ=255)
                      tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTgaga                                                  >  1:5681478/1‑76 (MQ=255)
                      tATTTTGCGCTGCTGAAAGTTAACGAAGTAAACTTCGACAAACCAGAAAACACCCGCAACAAAATCCTCTTTgaga                                                  <  2:3172193/76‑1 (MQ=255)
                       aTTTTNCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGaa                                                 >  2:2840412/1‑76 (MQ=255)
                       aTTTTGCGCTGCTGAATGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGaa                                                 <  1:3551903/76‑1 (MQ=255)
                       aTTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGaa                                                 <  1:6888798/76‑1 (MQ=255)
                          ttGCGCTGCTGAAAGTTAACGAAGTTAGCTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACtt                                              <  2:931068/76‑1 (MQ=255)
                          ttGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACtt                                              <  1:5837850/76‑1 (MQ=255)
                          ttGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACtt                                              >  2:661349/1‑76 (MQ=255)
                          ttGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACtt                                              >  1:3936767/1‑76 (MQ=255)
                           tGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACCCCCGCAACAAAATCCTCTTTGAGAACTTa                                             <  2:2487634/76‑1 (MQ=255)
                           tGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTa                                             >  2:7117376/1‑76 (MQ=255)
                            gcgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTaa                                            <  1:596738/76‑1 (MQ=255)
                             cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAAc                                           <  1:173415/76‑1 (MQ=255)
                             cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAAc                                           <  1:6532902/76‑1 (MQ=255)
                             cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAAc                                           <  2:395034/76‑1 (MQ=255)
                             cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAAc                                           <  2:1847602/76‑1 (MQ=255)
                              gctgctGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAAcc                                          <  2:4871640/76‑1 (MQ=255)
                                 gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                       >  2:735735/1‑76 (MQ=255)
                                 gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                       <  2:1773407/76‑1 (MQ=255)
                                  ctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGc                                      >  2:1607599/1‑76 (MQ=255)
                                    gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTg                                    >  2:3245516/1‑76 (MQ=255)
                                      aaGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCa                                  >  2:5014876/1‑76 (MQ=255)
                                      aaGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCa                                  <  2:3448573/76‑1 (MQ=255)
                                      aaGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCa                                  >  1:2049814/1‑76 (MQ=255)
                                       aGTTAACGTAGTCAACTTCGACAAACCTGAAAGCACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCAc                                 <  2:4501309/76‑1 (MQ=255)
                                       aGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCAc                                 <  1:5094990/76‑1 (MQ=255)
                                          tAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCa                              >  1:800513/1‑76 (MQ=255)
                                          tAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCa                              >  2:1934102/1‑76 (MQ=255)
                                           aaCGAATTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaa                             <  1:4838071/76‑1 (MQ=255)
                                           aaCGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaa                             >  1:2558934/1‑76 (MQ=255)
                                           aaCGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaa                             >  2:5790914/1‑76 (MQ=255)
                                           aaCGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaa                             >  1:3744725/1‑76 (MQ=255)
                                           aaCGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaa                             <  1:5208427/76‑1 (MQ=255)
                                            aCGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaaa                            <  2:5413181/76‑1 (MQ=255)
                                            aCGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaaa                            <  2:1636281/76‑1 (MQ=255)
                                             cGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAAc                           <  1:4438408/76‑1 (MQ=255)
                                              gAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAAct                          >  1:2827858/1‑76 (MQ=255)
                                              gAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAAct                          >  2:4011710/1‑76 (MQ=255)
                                                aGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActct                        >  1:7042649/1‑76 (MQ=255)
                                                 gTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                       <  2:7042749/76‑1 (MQ=255)
                                                  ttAACTTCTACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCg                      >  1:1867167/1‑76 (MQ=255)
                                                  ttAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCg                      >  1:6443982/1‑76 (MQ=255)
                                                  ttAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCg                      >  2:418220/1‑76 (MQ=255)
                                                  ttAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCg                      >  2:5601625/1‑76 (MQ=255)
                                                  ttAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCg                      >  1:3405925/1‑76 (MQ=255)
                                                   tAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGt                     <  2:5695774/76‑1 (MQ=255)
                                                   tAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGt                     >  2:6317006/1‑76 (MQ=255)
                                                     aCTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                   <  1:6029308/76‑1 (MQ=255)
                                                     aCTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                   <  1:2167525/76‑1 (MQ=255)
                                                     aCTTCCACAAACCTGAAACCACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                   <  1:6322194/76‑1 (MQ=255)
                                                      cTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTg                  <  1:3180226/76‑1 (MQ=255)
                                                       ttCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGc                 >  1:5428027/1‑76 (MQ=255)
                                                         cGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGt               >  1:720721/1‑76 (MQ=255)
                                                         cGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGt               <  2:7467669/76‑1 (MQ=255)
                                                          gACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTa              <  2:2752768/76‑1 (MQ=255)
                                                          gACAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTa              >  2:2775845/1‑76 (MQ=255)
                                                           aCAAACCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTAt             >  1:4996996/1‑76 (MQ=255)
                                                            cAAACCTGAAAACACCCGCAACAAAATCTTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATg            >  2:1368624/1‑76 (MQ=255)
                                                             aaaCCTGAAAACCCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATgg           <  1:927177/76‑1 (MQ=255)
                                                             aaaCCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATgg           >  2:3736734/1‑76 (MQ=255)
                                                             aaaCCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATgg           <  2:2370122/76‑1 (MQ=255)
                                                              aaCCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGa          <  2:6700511/76‑1 (MQ=255)
                                                              aaCCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGa          <  2:2206721/76‑1 (MQ=255)
                                                               aCCTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGaa         >  1:4379071/1‑76 (MQ=255)
                                                                ccTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAc        >  1:266196/1‑76 (MQ=255)
                                                                 cTGAAAACCCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACg       <  1:626247/76‑1 (MQ=255)
                                                                 cTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACg       <  2:6184487/76‑1 (MQ=255)
                                                                 cTGAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACg       >  2:2684706/1‑76 (MQ=255)
                                                                  taaaaaCACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGt      <  1:5858939/74‑1 (MQ=255)
                                                                   gAAAACACCCGCAACAAAATCCTCTTTGAGAACTTAATCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTg     <  2:567408/76‑1 (MQ=255)
                                                                     aaaCACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt   >  1:1418055/1‑76 (MQ=255)
                                                                     aaaCACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt   <  2:3851866/76‑1 (MQ=255)
                                                                     aaaCACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt   <  1:668775/76‑1 (MQ=255)
                                                                      aaCACCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTa  >  1:4712480/1‑76 (MQ=255)
                                                                         |                                                                        
CCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTA  >  REL606/3928522‑3928667

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