Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,322,724 G→A pseudogene (385/408 nt) yciQ → hypothetical protein; b1268_1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,322,7240GA98.9% 743.1 / ‑1.4 187pseudogene (385/408 nt)yciQhypothetical protein; b1268_1
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (104/81);  total (105/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

Only 100 of 187 total aligned reads displayed.

TCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTA  >  REL606/1322651‑1322791
                                                                         |                                                                   
tCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAAt                                                                   >  1:2609271/1‑76 (MQ=255)
tCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAAt                                                                   >  1:2472117/1‑76 (MQ=255)
   tGACTGGCATATCGGGTTAATGATTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAAtgct                                                                >  1:1123339/1‑76 (MQ=255)
   tGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAAtgct                                                                >  1:5488576/1‑76 (MQ=255)
   tGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAAtgct                                                                >  1:4428192/1‑76 (MQ=255)
      cTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGcc                                                             >  2:4820533/1‑76 (MQ=255)
       tGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCa                                                            <  1:3791120/76‑1 (MQ=255)
       tGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCa                                                            >  2:5223511/1‑76 (MQ=255)
        ggCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCac                                                           >  2:1801053/1‑75 (MQ=255)
        ggCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCaa                                                           >  1:4269088/1‑76 (MQ=255)
        ggCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCaa                                                           >  2:3427299/1‑76 (MQ=255)
        ggCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCaa                                                           <  1:3450552/76‑1 (MQ=255)
         gCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAAc                                                          <  1:3261823/76‑1 (MQ=255)
          cagatCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACt                                                         <  1:1002330/73‑1 (MQ=255)
          cATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACt                                                         <  1:4692685/76‑1 (MQ=255)
          cATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGGAATGCTGCCAACt                                                         >  2:2030918/1‑76 (MQ=255)
          cATATCGGGTTAATGAGTATTGTCGTGGGCCATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACt                                                         >  1:3293979/1‑76 (MQ=255)
           atatCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACtt                                                        >  2:4176852/1‑76 (MQ=255)
           atatCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACtt                                                        >  2:2209227/1‑76 (MQ=255)
            tatCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTa                                                       <  2:3836436/76‑1 (MQ=255)
             atCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTAc                                                      >  2:294341/1‑76 (MQ=255)
             atCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTAc                                                      >  1:5931713/1‑76 (MQ=255)
             atCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTAc                                                      <  1:6198718/76‑1 (MQ=255)
              tCGGGTTAGTGAGTATTGGCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACt                                                     <  1:3325878/76‑1 (MQ=39)
               cGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTg                                                    <  1:5772586/76‑1 (MQ=255)
               cGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTg                                                    <  1:6124886/76‑1 (MQ=255)
               cGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTg                                                    >  1:4179286/1‑76 (MQ=255)
                 ggTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGAt                                                  >  1:2424636/1‑76 (MQ=255)
                 ggTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGAt                                                  >  2:3208498/1‑76 (MQ=255)
                 ggTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGAt                                                  >  2:38673/1‑76 (MQ=255)
                  gTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGAtt                                                 >  2:4658640/1‑76 (MQ=255)
                  gTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGAtt                                                 >  2:4951856/1‑76 (MQ=255)
                  gTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGAtt                                                 >  2:1139559/1‑76 (MQ=255)
                   ttAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGAttt                                                >  1:1730864/1‑76 (MQ=255)
                    tAATGAGTTTTGTCGTGTGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTa                                               <  2:4686327/76‑1 (MQ=38)
                    tAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTa                                               <  1:6259848/76‑1 (MQ=255)
                    tAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTa                                               <  1:2180719/76‑1 (MQ=255)
                     aaTGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAg                                              <  2:819508/76‑1 (MQ=255)
                      aTGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAgt                                             <  1:1095528/76‑1 (MQ=255)
                      aTGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAgt                                             >  1:4668574/1‑76 (MQ=255)
                       tGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAgtg                                            <  2:1042058/76‑1 (MQ=255)
                       tGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAgtg                                            >  2:4484551/1‑76 (MQ=255)
                        gAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAgtgt                                           <  2:2473606/76‑1 (MQ=255)
                        gAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAgtgt                                           >  1:1626120/1‑76 (MQ=255)
                          gTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTat                                         <  2:2327502/76‑1 (MQ=255)
                          gTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTat                                         <  1:5101301/76‑1 (MQ=255)
                           tATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTatg                                        >  1:5599347/1‑76 (MQ=255)
                           tATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTatg                                        <  2:5661278/76‑1 (MQ=255)
                           tATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTatg                                        >  2:5075822/1‑76 (MQ=255)
                             ttGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTatgat                                      <  1:4086397/76‑1 (MQ=255)
                             ttGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTatgat                                      >  2:1236331/1‑76 (MQ=255)
                              tGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTatgatg                                     >  2:532587/1‑76 (MQ=255)
                              tGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTatgatg                                     >  2:5124307/1‑76 (MQ=255)
                               gTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATgg                                    >  1:3175662/1‑76 (MQ=255)
                               gTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATgg                                    >  1:438902/1‑76 (MQ=255)
                                tCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGgt                                   <  1:3607415/76‑1 (MQ=255)
                                    gggCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttt                               <  2:3989758/76‑1 (MQ=255)
                                    gggCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttt                               <  2:888469/76‑1 (MQ=255)
                                    gggCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttt                               >  1:2025449/1‑76 (MQ=255)
                                    gggCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttt                               >  1:6025017/1‑76 (MQ=255)
                                    gggCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttt                               <  1:4630595/76‑1 (MQ=255)
                                     ggCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGtttt                              <  1:5883198/76‑1 (MQ=255)
                                      gCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttttt                             <  1:3480784/76‑1 (MQ=255)
                                      gCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttttt                             >  1:3994332/1‑76 (MQ=255)
                                       cGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                            <  2:1645859/76‑1 (MQ=255)
                                       cGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                            >  1:1563571/1‑76 (MQ=255)
                                       cGATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                            <  2:210337/76‑1 (MQ=255)
                                        gATAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGa                           <  1:4891248/76‑1 (MQ=255)
                                         aTAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAg                          >  2:171523/1‑76 (MQ=255)
                                          tAAACAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAgg                         >  2:3881975/1‑76 (MQ=255)
                                           aaaCAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                        >  1:3765402/1‑76 (MQ=255)
                                           aaaCAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                        >  1:1268672/1‑76 (MQ=255)
                                           aaaCAACGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                        >  1:1431741/1‑76 (MQ=255)
                                            aacaacGTTTCTTGCCTCAAGGCGACGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTg                       >  2:612135/1‑76 (MQ=255)
                                            aacaacGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTg                       <  2:2841830/76‑1 (MQ=255)
                                            aacaacGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTg                       >  2:3119487/1‑76 (MQ=255)
                                              caacGTTTCTTGCCTCAAGGCGACGGTACTGCTGCCAACTTACTGATTTAGTGTATGATGGTGCTTTTGAGGCGCt                     >  2:264436/1‑76 (MQ=18)
                                                acGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTcc                   >  1:9558/1‑76 (MQ=255)
                                                acGTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTcc                   <  2:4747216/76‑1 (MQ=255)
                                                  gTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAg                 <  2:746685/76‑1 (MQ=255)
                                                  gTTTCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAg                 >  2:4525472/1‑76 (MQ=255)
                                                   tttCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGGGTTTTTGAGGTGCTCCAGt                >  2:255433/1‑76 (MQ=255)
                                                    ttCTTGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTg               <  1:102231/76‑1 (MQ=255)
                                                       ttGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt            >  1:2314852/1‑76 (MQ=255)
                                                        tGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt           >  1:4919631/1‑76 (MQ=255)
                                                        tGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt           >  1:1937443/1‑76 (MQ=255)
                                                        tGCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt           <  1:461542/76‑1 (MQ=255)
                                                         gCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGGGGTGCTCCAGTGGCttt          >  1:140857/1‑75 (MQ=255)
                                                         gCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt            <  2:2678387/74‑1 (MQ=255)
                                                         gCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt            >  1:2678387/1‑74 (MQ=255)
                                                         gCCTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc          <  2:4185970/76‑1 (MQ=255)
                                                          ccTCAAGGTGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCt         <  1:2920466/76‑1 (MQ=255)
                                                           cTCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg        <  1:5962812/76‑1 (MQ=255)
                                                            tCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGt       >  1:1532767/1‑76 (MQ=255)
                                                            tCAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGt       >  2:5932842/1‑76 (MQ=255)
                                                             cAAGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGtt      >  2:3891597/1‑76 (MQ=255)
                                                              aaGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt     <  2:6234577/76‑1 (MQ=255)
                                                              aaGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt     <  2:2407882/76‑1 (MQ=255)
                                                              aaGGCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt     <  2:4572450/76‑1 (MQ=255)
                                                                 gCGACGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTa  <  1:3208498/76‑1 (MQ=17)
                                                                         |                                                                   
TCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTA  >  REL606/1322651‑1322791

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