Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 200,540 C→T 5.9% L385F (CTT→TTT)  yaeL → zinc metallopeptidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606200,5400C→T 5.9% 11.3 0.0 ‑0.0153L385F (CTT→TTT) yaeLzinc metallopeptidase
Reads supporting (aligned to +/- strand):  new base (4/5):  ref base (60/84):  total (64/89)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.61e-01

Only 100 of 172 total aligned reads displayed.

ACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTT  >  REL606/200506‑200574
                                  |                                  
 cTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCt                                   >  1:4310794/1‑35 (MQ=255)
 cTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCt                                   >  1:182494/1‑35 (MQ=255)
  tCGGGGTTGTTTATTCCCTGCCGTTTCTTGCGTtt                                  <  1:17445357/35‑1 (MQ=37)
  tCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCtt                                  >  1:4562648/1‑35 (MQ=255)
   ggggtTTGTTTATTCCCTGCCGTTTCTTGCGCtta                                 <  1:6568817/34‑1 (MQ=25)
   cGGGTTTTTTTATTCCCTGCCGTTTCTTGCGCtta                                 <  1:18700411/35‑1 (MQ=25)
   cGGGGTTGTTTATTCCCTGCCGTTTCTTGCGCtta                                 <  1:8959641/35‑1 (MQ=255)
   cGGGGTTGTTTATTACCTGCCGTTTCTTGCGCtta                                 >  1:3111840/1‑35 (MQ=255)
    ggggTTGTTTATTACCTGCCGTTTCTTGCGCttat                                >  1:329602/1‑35 (MQ=255)
     ggtttgttTATTCCCTGCCGTTTCTTGCGCttatt                               <  1:7667097/32‑1 (MQ=37)
     ggtttgttTATTACCTGCCGTTTCTTGCGCttatt                               >  1:8823015/4‑35 (MQ=255)
     gggTTGTTTATTACCTGCCGTTTCTTGCGCttatt                               >  1:10581711/1‑35 (MQ=255)
     gggTTGTTTATTACCTGCCGTTTCTTGCGCttatt                               >  1:494595/1‑35 (MQ=255)
     gggTTGTTTATTACCTGCCGTTTCTTGCGCttatt                               <  1:6776674/35‑1 (MQ=255)
      ggTTGTTTATTACCTGCCGGTTCTTGCGCttttta                              >  1:19470836/1‑35 (MQ=37)
       tttgttTATTACCTGCCGTTTCTTGCGCTTATTAg                             <  1:14841305/34‑1 (MQ=255)
       tttgttTATTACCTGCCGTTTCTTGCGCTTATTAg                             <  1:13863930/34‑1 (MQ=255)
       gttgttTATTACCTGCCGTTTCTTGCGCTTATTAg                             >  1:2967075/1‑35 (MQ=255)
       gttgttTATTACCTGCCGTTTCTTGCGCTTATTAg                             >  1:16911207/1‑35 (MQ=255)
        ttgttTATTCCCTGCCGTTTCTTGCGCTTATTAGc                            <  1:7172478/35‑1 (MQ=255)
        ttgttTATTCCCTGCCGTTTCTTGCGCTTATTAGc                            <  1:3692733/35‑1 (MQ=255)
        ttgttTATTACCTGCCGTTTCTTGCGCTTATTAGc                            >  1:20341222/1‑35 (MQ=255)
         tgttTATTCCCTGCCGTTTCTTGCGCTTATTAGCg                           <  1:915901/35‑1 (MQ=255)
         tgttTATTACCTGCCGTTTCTTGCGTTTATTAGCg                           <  1:16453471/35‑1 (MQ=255)
         tgttTATTACCTGCCGTTTCTTGCGCTTATTAGCg                           >  1:17795633/1‑35 (MQ=255)
          ttttATCCCCTGCCGTTTCTTGCGCTTATTAGCGt                          <  1:2159550/34‑1 (MQ=25)
          gttTATTCCCTGCCGTTTCTTGCGCTTATTAGCGt                          <  1:14351112/35‑1 (MQ=255)
          gttTATTCCCTGCCGTTTCTTGCGCTTATTAGCGt                          <  1:5287180/35‑1 (MQ=255)
          gttTATTCCCTGCCGTTTCTTGCGCTTATTAGCGt                          <  1:6258/35‑1 (MQ=255)
          gttTATTACCTGCCGTTTCTTGCGCTTATTAGCGt                          >  1:1034937/1‑35 (MQ=255)
           tttATTCCCTGCCGTTTCTTGCGCTTATTAGCGTg                         <  1:2811160/35‑1 (MQ=255)
           tttATTACCTGCCGTTTCTTGCGCTTATTAGCGTg                         >  1:18333046/1‑35 (MQ=255)
            ttattgCCTGCCGTTTCTTGCGCTTATTAGCGTGa                        <  1:20219361/35‑1 (MQ=255)
            ttattaCCTGCCGTTTCTTGCGCTTATTAGCGTGa                        >  1:13720420/1‑35 (MQ=255)
             tattaCCTGCCGTTTCTTGCGCTTATTAGCGTGaa                       >  1:8440493/1‑35 (MQ=255)
              attcCCTGCCGTTTCTTGCGCTTATTAGCGTGGAc                      <  1:14833818/35‑1 (MQ=37)
              attcCCTGCCGTTTCTTGCGCTTATTAGCGTGAAc                      <  1:404180/35‑1 (MQ=255)
              attcCCTGCCGTTTCTTGCGCTTATTAGCGTGAAc                      <  1:14933968/35‑1 (MQ=255)
               ttaCCTGCCGTTTCTTGCGCTTATTAGCGTGAACt                     >  1:14263645/1‑35 (MQ=255)
                tcccTGCCGTTTCTTGCGCTTATTAGCGTGAACtt                    <  1:17225312/33‑1 (MQ=255)
                tcccTGCCGTTTCTTGCGCTTATTAGCGTGAACtt                    <  1:8328741/33‑1 (MQ=255)
                taCCTGCCGTTTCTTGCGCTTATTAGAGTGAACtt                    <  1:19062146/35‑1 (MQ=255)
                 cccTGCCGTTTCTTGCGCTTATTAGCGTGAACTTa                   <  1:16665852/34‑1 (MQ=255)
                 cccTGCCGTTTCTTGCGCTTATTAGCGTGAACTTa                   <  1:15264053/34‑1 (MQ=255)
                 aCCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTa                   >  1:14692554/1‑35 (MQ=255)
                  ccTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAg                  <  1:11931058/35‑1 (MQ=255)
                  ccTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAg                  >  1:5880311/1‑35 (MQ=255)
                  ccTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAg                  >  1:10701733/1‑35 (MQ=255)
                  ccTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAg                  >  1:13266352/1‑35 (MQ=255)
                  ccTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAg                  >  1:19198160/1‑35 (MQ=255)
                   cTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAgg                 <  1:10037848/35‑1 (MQ=255)
                   cTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAgg                 <  1:13432450/35‑1 (MQ=255)
                    tGCCGTTTCTTGCGCTTATTAGCGTGAACTTAggg                <  1:11527405/35‑1 (MQ=255)
                    tGCCGTTTCTTGCGCTTATTAGCGTGAACTTAggg                <  1:1549718/35‑1 (MQ=255)
                     gCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGa               >  1:8634716/1‑35 (MQ=255)
                     gCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGa               <  1:311601/35‑1 (MQ=255)
                     gCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGa               >  1:15423011/1‑35 (MQ=255)
                      ccGTTTCTTGCGTTTATTAGCGTGAACTTAGGGAt              >  1:18955868/1‑35 (MQ=255)
                       cGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATa             >  1:9878254/1‑35 (MQ=255)
                       cGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATa             >  1:10797717/1‑35 (MQ=255)
                       cGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATa             >  1:12256527/1‑35 (MQ=255)
                       cGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATa             >  1:208946/1‑35 (MQ=255)
                       cGTTTCTTGCGCTTATTAGCGTGAACTTAGCGATa             >  1:13182437/1‑35 (MQ=255)
                        gTTTCTTGCGTTTATTAGCGTGAACTTAGGGATaa            <  1:7757001/35‑1 (MQ=255)
                        gTTTCTTGCGCTTATTAGCGTGAACTTAGGGATaa            >  1:7481256/1‑35 (MQ=255)
                        gTTTCTTGCGCTTATTAGCGTGAACTTAGGGATaa            >  1:16445862/1‑35 (MQ=255)
                        gTTTCTTGCGCTTATTAGCGTGAACTTAGGGATaa            >  1:14575109/1‑35 (MQ=255)
                        gTTTCTTGCGCTTATTAGCGTGAACTTAGGGATaa            >  1:7004092/1‑35 (MQ=255)
                         tttCTTGCGCTTATTAGCGTGAACTTANGGATAAt           <  1:5713630/35‑1 (MQ=255)
                         tttCTTGCGCTTATTAGCGTGAACTTAGGGATAAt           >  1:2669298/1‑35 (MQ=255)
                         tttCTTGCGCTTATTAGCGTGAACTTAGGGATAAt           <  1:12261567/35‑1 (MQ=255)
                          ttCTTGCGCTTATTAGCGTGAACTTAGGGATAAtt          <  1:2204933/35‑1 (MQ=255)
                          ttCTTGCGCTTATTAGCGTGAACTTAGGGATAAtt          <  1:20382676/35‑1 (MQ=255)
                          ttCTTGCGCTTATTAGCGTGAACTTAGGGATAAtt          <  1:3922877/35‑1 (MQ=255)
                          ttCTTGCGCTTATTAGCGTGAACTTAGGGATAAtt          <  1:13406729/35‑1 (MQ=255)
                          ttCTTGCGCTTATTAGAGTGAACTTAGGGATAAtt          <  1:19088911/35‑1 (MQ=255)
                           tCTTGCGTTTATTAGCGTGAACTTAGGGATAATTa         <  1:256838/35‑1 (MQ=255)
                           tCTTGCGCTTATTAGCGTGAACTTAGGGATAATTa         <  1:5095913/35‑1 (MQ=255)
                           tCTTGCGCTTATTAGCGTGAACTTAGGGATAATTa         >  1:5540806/1‑35 (MQ=255)
                           tCTTGCGCTTATTAGCGTGAACTTAGGGATAATTa         <  1:13130608/35‑1 (MQ=255)
                           tCTTGCGCTTATTAGCGTGAACTTAGGGATAATTa         >  1:18768811/1‑35 (MQ=255)
                           tCTTGCGCTTATTAGCGTGAACTTAGGGATAATTa         <  1:11300148/35‑1 (MQ=255)
                           tCTTGCGCTTATTAGCGTGAACTTAGGGATAATTa         <  1:200884/35‑1 (MQ=255)
                            cTTGCGCTTATTAGCGTGAACTTAGGGATAATTaa        >  1:13796221/1‑35 (MQ=255)
                             ttGCGCTTATTAGCGTGAACTTAGGGATAATTAAc       >  1:2371477/1‑35 (MQ=255)
                             ttGCGCTTATTAGCGTGAACTTAGGGATAATTAAc       >  1:10167284/1‑35 (MQ=255)
                             ttGCGCTTATTAGCGTGAACTTAGGGATAATTAAc       >  1:8502549/1‑35 (MQ=255)
                             ttGCGCTTATTAGCGTGAACTTAGGGATAATTAAc       <  1:7848494/35‑1 (MQ=255)
                             ttGCGCTTATTAGCGTGAACTTAGGGATAATTAAc       >  1:16413632/1‑35 (MQ=255)
                             ttGCGCTTATTAGCGTGAACTTAGGGATAATTAAc       <  1:18070353/35‑1 (MQ=255)
                             ttGCGCTTATTAGCGTGAACTTAGGGATAATTAAc       >  1:12641723/1‑35 (MQ=255)
                              tGCGCTTATTGGCGTGAACTTAGGGATAATTAAcc      <  1:889699/35‑1 (MQ=255)
                              tGCGCTTATTAGCGTGAACTTAGGGATAATTAAcc      >  1:11815698/1‑35 (MQ=255)
                               gcgcTTATTAGCGTGAACTTAGGGATAATTAACCt     <  1:14335513/35‑1 (MQ=255)
                               gcgcTTATTAGCGTGAACTTAGGGATAATTAACCt     <  1:19759108/35‑1 (MQ=255)
                               gcgcTTATTAGCGTGAACTTAGGGATAATTAACCt     <  1:18713172/35‑1 (MQ=255)
                               gcgcTTATTAGCGTGAACTTAGGGATAATTAACCt     <  1:6475976/35‑1 (MQ=255)
                                cgcTTATTAGCGTGAACTTAGGGATAATTAACCTg    <  1:1260417/35‑1 (MQ=255)
                                 gcTTATTAGCGTGAACTTAGGGATAATTAACCTGt   <  1:14553579/35‑1 (MQ=255)
                                 gcTTATTAGCGTGAACTTAGGGATAATTAACCTGt   >  1:1724743/1‑35 (MQ=255)
                                  |                                  
ACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTT  >  REL606/200506‑200574

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.