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.0% 181.5 / ‑2.9 50I15N (ATC→AAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (18/30);  total (19/31)
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.98e-01

TACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928189‑3928332
                                                                       |                                                                        
tACCCCAAGTTTAAGAACTCGCACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTAAGCTGAACAct                                                                      >  2:134564/1‑76 (MQ=255)
tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                      >  2:166713/1‑76 (MQ=255)
 aCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActc                                                                     >  2:778582/1‑76 (MQ=255)
  ccccAAGTTTAAGAACTCACGCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCActct                                                                    >  2:132714/1‑76 (MQ=255)
   cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   <  2:426218/76‑1 (MQ=255)
   cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   <  2:426238/76‑1 (MQ=255)
   cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   >  1:547191/1‑76 (MQ=255)
   cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   <  1:1435949/76‑1 (MQ=255)
      aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGc                                                                <  1:573946/76‑1 (MQ=255)
      aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGc                                                                <  1:1193297/76‑1 (MQ=255)
       aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACANTCTCGGCg                                                               <  1:1650975/76‑1 (MQ=255)
       aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               <  2:517199/76‑1 (MQ=255)
       aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               <  1:1168126/76‑1 (MQ=255)
       aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               <  1:953431/76‑1 (MQ=255)
       aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               <  1:227925/76‑1 (MQ=255)
         tttNAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             >  1:1537090/1‑76 (MQ=255)
         tttAAGAACTCTCACCACTATGAATCTTACCGAATTAACGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             >  1:100463/1‑76 (MQ=255)
         tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             <  2:1079572/76‑1 (MQ=255)
           tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           >  1:220900/1‑76 (MQ=255)
             aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         >  2:1330604/1‑76 (MQ=255)
                  tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg                                                    >  2:424031/1‑76 (MQ=255)
                  tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg                                                    >  2:818560/1‑76 (MQ=255)
                   cacaCCANCATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGc                                                   >  2:1390885/1‑76 (MQ=255)
                    acacCACTATGAATCTTACCGAGTTAANGANTGNGCTGGTTNCTTAGCAGAACACTCTCGGCGAAAATATGGGGCt                                                  <  1:1478716/76‑1 (MQ=255)
                    acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt                                                  >  1:456453/1‑76 (MQ=255)
                     caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                 <  1:1332272/76‑1 (MQ=255)
                      accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                <  1:599341/76‑1 (MQ=255)
                      accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATANGGGGCTgg                                                <  1:1642481/76‑1 (MQ=255)
                       ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGa                                               >  1:590854/1‑76 (MQ=255)
                         acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa                                             <  2:925724/76‑1 (MQ=255)
                          cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaaa                                            <  2:644040/76‑1 (MQ=255)
                            aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          <  2:1510004/76‑1 (MQ=255)
                            aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  1:172823/1‑76 (MQ=255)
                                     aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                                 <  1:744783/76‑1 (MQ=255)
                                      ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  2:122988/76‑1 (MQ=255)
                                        gAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc                              <  1:410046/76‑1 (MQ=255)
                                         aaTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg                             <  2:751780/76‑1 (MQ=255)
                                           ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           <  2:1321978/76‑1 (MQ=255)
                                           ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           >  2:843194/1‑76 (MQ=255)
                                            tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  2:732770/76‑1 (MQ=255)
                                                   tACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTGGTATGCGTAAGCAGGAc                   <  1:1704266/76‑1 (MQ=255)
                                                     cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 >  1:1047998/1‑76 (MQ=255)
                                                     cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 <  1:1160693/76‑1 (MQ=255)
                                                     cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 <  1:330035/76‑1 (MQ=255)
                                                     cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 <  1:465322/76‑1 (MQ=255)
                                                     cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 <  1:725398/76‑1 (MQ=255)
                                                          gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt            <  1:1516696/76‑1 (MQ=255)
                                                           tttCTGAGCTTAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt           <  2:365121/76‑1 (MQ=255)
                                                               tGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGcc       >  1:347366/1‑76 (MQ=255)
                                                                    tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  1:298880/1‑76 (MQ=255)
                                                                       |                                                                        
TACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928189‑3928332

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