Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,343 G→A A43T (GCA→ACA)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,3430GA99.6% 951.5 / ‑5.4 224A43T (GCA→ACA) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (110/113);  total (110/114)
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

Only 100 of 236 total aligned reads displayed.

CGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTT  >  REL606/3928270‑3928415
                                                                         |                                                                        
cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGcacaca                                                                        <  2:2482233/76‑1 (MQ=255)
cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGcacaca                                                                        <  2:1353931/76‑1 (MQ=255)
  aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACCCaaa                                                                      >  2:3740872/1‑76 (MQ=255)
  aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACaaa                                                                      <  2:4436388/76‑1 (MQ=255)
  aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACaaa                                                                      <  2:5409855/76‑1 (MQ=255)
  aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACaaa                                                                      <  2:2370997/76‑1 (MQ=255)
   aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAag                                                                     <  1:1121178/76‑1 (MQ=255)
   aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAag                                                                     >  2:6226419/1‑76 (MQ=255)
   aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAag                                                                     <  1:59620/76‑1 (MQ=255)
     atatGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAagag                                                                   >  1:3595093/1‑76 (MQ=255)
     atatGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAagag                                                                   >  2:1208639/1‑76 (MQ=255)
       atGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCCCACAAagaggg                                                                 >  1:2754526/1‑74 (MQ=255)
       atGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTg                                                                 <  1:4897469/76‑1 (MQ=255)
        tGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCCTCCTGAAGCAGCACACCAAGAgggg                                                                >  1:3321952/1‑76 (MQ=255)
        tGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTgg                                                                <  2:2654583/76‑1 (MQ=255)
         ggggCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGc                                                               <  1:3878329/76‑1 (MQ=255)
          gggCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCg                                                              <  1:1515731/76‑1 (MQ=255)
          gggCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCg                                                              <  2:4290645/76‑1 (MQ=255)
           ggCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGa                                                             <  1:4772815/76‑1 (MQ=255)
            gCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGaa                                                            <  2:1645461/76‑1 (MQ=255)
             cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAg                                                           <  1:4579578/76‑1 (MQ=255)
              tGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGa                                                          <  2:1108520/76‑1 (MQ=255)
              tGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGa                                                          >  2:6001729/1‑76 (MQ=255)
               ggAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGat                                                         <  1:3402613/76‑1 (MQ=255)
                 aaaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGatat                                                       <  2:5819169/76‑1 (MQ=255)
                 aaaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGatat                                                       <  1:5675695/76‑1 (MQ=255)
                   aaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCt                                                     <  1:138164/76‑1 (MQ=255)
                     ccTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCCTCCTGAAGCAGCACACAAAGAGTGGCGGAGATATCttt                                                   >  2:458510/1‑76 (MQ=255)
                       tGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTgg                                                 <  2:2757149/76‑1 (MQ=255)
                        ggCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTggg                                                >  2:3116325/1‑75 (MQ=255)
                        ggCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGt                                                <  1:2149537/76‑1 (MQ=255)
                        ggCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGAAATCttttgg                                                >  2:337311/1‑73 (MQ=255)
                         gCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTg                                               <  1:1816457/76‑1 (MQ=255)
                         gCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTg                                               >  2:1928287/1‑76 (MQ=255)
                          cTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGa                                              >  1:4514928/1‑76 (MQ=255)
                          cTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGa                                              >  2:3617458/1‑76 (MQ=255)
                          cTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGa                                              >  1:5323279/1‑76 (MQ=255)
                            cGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATg                                            >  2:1630874/1‑76 (MQ=255)
                             gTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATgg                                           >  1:608742/1‑76 (MQ=255)
                              tATGCGTAAGCAGGACATTATTTTTGTCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGc                                          >  2:6365291/1‑76 (MQ=255)
                              tATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGc                                          >  1:5180188/1‑76 (MQ=255)
                              tATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGc                                          <  2:2933938/76‑1 (MQ=255)
                              tATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGc                                          >  2:2157721/1‑76 (MQ=255)
                              tATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGc                                          <  1:2472848/76‑1 (MQ=255)
                              tATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGc                                          >  2:6895905/1‑76 (MQ=255)
                               aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                         >  1:2374294/1‑76 (MQ=255)
                               aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                         >  1:4327407/1‑76 (MQ=255)
                               aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                         <  1:2775310/76‑1 (MQ=255)
                               aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                         >  2:3063351/1‑76 (MQ=255)
                               aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                         <  2:2782654/76‑1 (MQ=255)
                                 gCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGTGTa                                       <  1:3350881/76‑1 (MQ=255)
                                 gCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTa                                       <  2:5304353/76‑1 (MQ=255)
                                  cGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGGACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTAc                                      >  1:6391227/1‑76 (MQ=255)
                                    tAAGCAGGACATTATTTTTGCCATCCTGAAGGAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTg                                    >  1:5805392/1‑76 (MQ=255)
                                    tAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTg                                    >  1:4665827/1‑76 (MQ=255)
                                     aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTgg                                   >  2:5001064/1‑76 (MQ=255)
                                     aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTgg                                   >  2:5939052/1‑76 (MQ=255)
                                     aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTgg                                   <  2:3982926/76‑1 (MQ=255)
                                       gCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgag                                 <  1:3865215/76‑1 (MQ=255)
                                        cAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgaga                                <  2:1755867/76‑1 (MQ=255)
                                        cAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgaga                                >  2:4279936/1‑76 (MQ=255)
                                         aGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGat                               >  2:4133138/1‑76 (MQ=255)
                                         aGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGat                               <  1:1682468/76‑1 (MQ=255)
                                         aGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGat                               <  2:2321603/76‑1 (MQ=255)
                                          ggACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGata                              >  2:3001188/1‑76 (MQ=255)
                                            aCATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATAtt                            >  2:4472531/1‑76 (MQ=255)
                                               ttattTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCa                         <  1:6513596/76‑1 (MQ=255)
                                                tattTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAg                        <  2:3263887/76‑1 (MQ=255)
                                                 attTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAgg                       <  2:1909393/76‑1 (MQ=255)
                                                  tttttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAgga                      <  2:729511/76‑1 (MQ=255)
                                                   ttttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggat                     <  2:3925381/76‑1 (MQ=255)
                                                   ttttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggat                     <  2:4819127/76‑1 (MQ=255)
                                                   ttttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggat                     >  2:240815/1‑76 (MQ=255)
                                                    tttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatg                    <  2:2248392/76‑1 (MQ=255)
                                                    tttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatg                    <  2:2541029/76‑1 (MQ=255)
                                                     ttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgg                   >  1:3089018/1‑76 (MQ=255)
                                                     ttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgg                   <  1:1764038/76‑1 (MQ=255)
                                                     ttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgg                   >  1:6760716/1‑76 (MQ=255)
                                                     ttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAAATATTGCAggatgg                   >  1:1453574/1‑76 (MQ=255)
                                                      tGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgga                  >  1:1076688/1‑76 (MQ=255)
                                                      tGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgga                  >  2:866086/1‑76 (MQ=255)
                                                      tGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgga                  <  1:1955324/76‑1 (MQ=255)
                                                      tGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgga                  >  1:3962204/1‑76 (MQ=255)
                                                       gCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatggat                 <  1:4038237/76‑1 (MQ=255)
                                                       gCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatggat                 >  1:4429135/1‑76 (MQ=255)
                                                        ccATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGAtt                >  1:3263856/1‑76 (MQ=255)
                                                          aTCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTg              <  2:922675/76‑1 (MQ=255)
                                                          aTCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTg              <  1:1970198/76‑1 (MQ=255)
                                                           tCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTgg             >  1:1982755/1‑76 (MQ=255)
                                                            ccTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGt            <  1:3766308/76‑1 (MQ=255)
                                                               gAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGgttt         >  1:5445837/1‑75 (MQ=255)
                                                               gAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGGTtc         >  2:787175/1‑76 (MQ=255)
                                                                aaGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTTGttttc        >  2:2832972/1‑74 (MQ=255)
                                                                aaGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGttttc        >  2:6195799/1‑74 (MQ=255)
                                                                 agcagcACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGGTTTGGTTtcct       >  1:4457224/1‑76 (MQ=255)
                                                                 agcagcACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcct       <  2:4871708/76‑1 (MQ=255)
                                                                 agcagcACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGGTtcct       >  1:1020525/1‑76 (MQ=255)
                                                                 agcagcACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGCAGATATTGCAGGATGGATTTGGTTtcct       <  2:3810866/76‑1 (MQ=255)
                                                                   cagcaCACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctcc     >  1:5039727/1‑76 (MQ=255)
                                                                    agcaCACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCg    >  1:6934215/1‑76 (MQ=255)
                                                                         |                                                                        
CGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTT  >  REL606/3928270‑3928415

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