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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,5930AT93.6% 160.4 / 0.1 47N126I (AAC→ATC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (24/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGG  >  REL606/3928521‑3928665
                                                                        |                                                                        
gcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGc                                                                       >  1:139173/1‑76 (MQ=255)
 cccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGcc                                                                      >  1:587338/1‑76 (MQ=255)
    gccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGc                                                                   <  1:657926/76‑1 (MQ=255)
      cgAAAGAAGGTAAACGCTATTTTGCGCTGCTGAAAGGTAACGAAGTTAACTTCGACAAACAAAAAATCGCCCGcaa                                                                 <  2:1366485/76‑1 (MQ=255)
       gAAAGAAGGTGAACGTTATTTTGCGCTGCTAAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGcaac                                                                <  1:1400705/76‑1 (MQ=255)
            aagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAAt                                                           <  1:234555/76‑1 (MQ=255)
             agGTGAACGCTATTTTGCGCTGCTGAAAGTTAACAAATTAACCTTCGACAAACCTGAAATCGCCCGCAACAAAATc                                                          <  2:932633/76‑1 (MQ=255)
                   aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCttt                                                    <  2:1316398/76‑1 (MQ=255)
                    cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTg                                                   <  1:1366599/76‑1 (MQ=255)
                    cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTg                                                   >  1:1556822/1‑76 (MQ=255)
                      cTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTgag                                                 >  1:262021/1‑76 (MQ=255)
                      cTATTTTGCGCTGCTAAAAGTTAACGAATTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTgag                                                 <  1:380663/76‑1 (MQ=255)
                       tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTgaga                                                >  1:1232215/1‑76 (MQ=255)
                        aTTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGaa                                               >  2:344805/1‑76 (MQ=255)
                        aTTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGaa                                               >  2:517059/1‑76 (MQ=255)
                        aTTTTGCGCTGATGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGaa                                               >  1:1307249/1‑76 (MQ=255)
                        aTTTTGAGCTGCTGAAAGTTAACGAGGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGaa                                               <  1:103603/76‑1 (MQ=255)
                         ttttGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAAc                                              >  2:985344/1‑76 (MQ=255)
                           ttGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACtt                                            <  2:1020826/76‑1 (MQ=255)
                             gcgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTaa                                          <  1:432133/76‑1 (MQ=255)
                              cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAAc                                         >  1:524800/1‑76 (MQ=255)
                              cgcTGCTGAAAGTAAACGAATTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAAc                                         <  2:759948/76‑1 (MQ=255)
                               gctgctGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAAcc                                        >  1:1688926/1‑76 (MQ=255)
                               gctgctCAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAAcc                                        >  2:908686/1‑76 (MQ=255)
                                ctgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAAccc                                       >  1:803132/1‑76 (MQ=255)
                                 tgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAAcccc                                      <  2:1714402/76‑1 (MQ=255)
                                  gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                     >  1:909576/1‑76 (MQ=255)
                                  gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                     >  1:909595/1‑76 (MQ=255)
                                  gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                     >  1:131582/1‑76 (MQ=255)
                                  gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                     >  2:1624107/1‑76 (MQ=255)
                                    tgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCt                                   >  2:422826/1‑76 (MQ=255)
                                     gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTg                                  >  2:202719/1‑76 (MQ=255)
                                           tAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCa                            >  2:971673/1‑76 (MQ=255)
                                           tAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCa                            >  1:1745897/1‑76 (MQ=255)
                                              cGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAAc                         >  1:1558641/1‑76 (MQ=255)
                                                aaGTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActc                       <  2:275650/76‑1 (MQ=255)
                                                aaGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActc                       >  2:430257/1‑76 (MQ=255)
                                                  gTTAACTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                     <  1:76160/76‑1 (MQ=255)
                                                  gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                     >  2:1702843/1‑76 (MQ=255)
                                                       cTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTg                <  1:1310433/76‑1 (MQ=255)
                                                       cTTCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTg                <  2:1053783/76‑1 (MQ=255)
                                                        ttCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGc               <  2:1073674/76‑1 (MQ=255)
                                                        ttCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGc               >  2:1052025/1‑76 (MQ=255)
                                                         tCGACAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCg              <  1:440310/76‑1 (MQ=255)
                                                             cAAACCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATg          <  1:1654576/76‑1 (MQ=255)
                                                                aCCTGAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGaa       <  1:908686/76‑1 (MQ=255)
                                                                    gAAATCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTg   >  2:572234/1‑76 (MQ=255)
                                                                     aaatCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg  >  2:59527/1‑76 (MQ=255)
                                                                        |                                                                        
GCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGG  >  REL606/3928521‑3928665

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