Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,665,639 A→T 92.4% L100Q (CTG→CAG)  rplS ← 50S ribosomal protein L19

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,665,6390A→T92.4% 18.5 ‑1.2 ‑0.0132L100Q (CTG→CAG) rplS50S ribosomal protein L19
Reads supporting (aligned to +/- strand):  new base (71/51):  ref base (8/2):  total (79/53)
Fisher's exact test strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.67e-01

Only 100 of 134 total aligned reads displayed.

TTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCAC  >  REL606/2665605‑2665672
                                  |                                 
tttGATACGAGCAGCCTTACCAGTACGCTCACGCa                                   >  1:5515046/1‑35 (MQ=255)
  tGATACGAGCAGCCTTACCAGTACGCTCACGCTgg                                 >  1:13936008/1‑35 (MQ=255)
  tGATACGAGCAGCCTTACCAGTACGCTCACGCTgg                                 >  1:275721/1‑35 (MQ=255)
  tGATACGAGCAGCCTTACCAGTACGCTCACGCTgg                                 >  1:7186687/1‑35 (MQ=255)
  tGATACGAGCAGCCTTACCAGTACGCTCACGCTgg                                 <  1:7319290/35‑1 (MQ=255)
  tGATACGACCAGCCTTACCAGTACGCTCACGCTgg                                 <  1:4665737/35‑1 (MQ=37)
   gATACGAGCAGCCTTACCAGTACGCTCACGCTGgt                                <  1:3792791/35‑1 (MQ=255)
   gATACGAGCAGCCTTACCAGTACGCTCACGCTGgt                                >  1:5473548/1‑35 (MQ=255)
   gATACGAGCAGCCTTACCAGTACGCTCACGCTGgt                                >  1:7624979/1‑35 (MQ=255)
    aTACGAGCAGCCTTACCAGTACGCTCACGCTGgta                               >  1:11776932/1‑35 (MQ=255)
    aTACGAGCAGCCTTACCAGTACGCTCACGCTGgta                               >  1:7344232/1‑35 (MQ=255)
    aTACGAGCAGCCTTACCAGTACGCTCACGCTGgta                               >  1:3116700/1‑35 (MQ=255)
    aTACGAGCAGCCTTACCAGTACGCTCACGCTGgta                               >  1:4295322/1‑35 (MQ=255)
    aTACGAGCAGCCTTACCAGTACGCTCACGCTGgta                               >  1:14113331/1‑35 (MQ=255)
     tACGAGCAGCCTTACCAGTACGCTCACGCTGgtag                              <  1:2969940/35‑1 (MQ=255)
     tACGAGCAGCCTTACCAGTACGCTCACGCTGgtag                              >  1:2035697/1‑35 (MQ=255)
     tACGAGCAGCCTTACCAGTACGCTCACGCAGgtag                              >  1:16504179/1‑35 (MQ=255)
      aCGAGCAGCCTTACCAGTACGCTCACGCTGgtagt                             >  1:20589012/1‑35 (MQ=255)
      aCGAGCAGCCTTACCAGTACGCTCACGCTGgtagt                             >  1:11387645/1‑35 (MQ=255)
      aCGAGCAGCCTTACCAGTACGCTCACGCTGgtagt                             >  1:823277/1‑35 (MQ=255)
      aCGAGCAGCCTTACCAGTACGCTCACGCTGgtagt                             >  1:8353292/1‑35 (MQ=255)
      aCGAGCAGCCTTACCAGTACGCTCACGCTGgtagt                             <  1:14713774/35‑1 (MQ=255)
      aCGAGCAGCCTTACCAGTACGCTCACGCTGgtagt                             >  1:3376153/1‑35 (MQ=255)
      aCGAGCAGCCTTACCAGTACGCTCACGCAGgtagt                             <  1:18140914/35‑1 (MQ=255)
       cGAGCAGCCTTACCAGTACGCTCACGCTGgtagta                            >  1:4709377/1‑35 (MQ=255)
        gAGCAGCCTTACCAGTACGCTCACTCTGGTAGTAc                           >  1:18977453/1‑35 (MQ=37)
        gAGCAGCCTTACCAGTACGCTCACGCTGGTAGTAc                           <  1:15042608/35‑1 (MQ=255)
        gAGCAGCCTTACCAGTACGCTCACGCTGGTAGTAc                           <  1:6427937/35‑1 (MQ=255)
        gAGCAGCCTTACCAGTACGCTCACGCTGGTAGTAc                           <  1:11410283/35‑1 (MQ=255)
         agcagcCTTACCAGTACGCTCACGCTGGTAGTACa                          >  1:12173694/1‑35 (MQ=255)
         agcagcCTTACCAGTACGCTCACGCTGGTAGTACa                          >  1:2531972/1‑35 (MQ=255)
         agcagcCTTACCAGTACGCTCACGCTGGTAGTACa                          <  1:10352679/35‑1 (MQ=255)
         agcagcCTTACCAGTACGCTCACGCTGGTAGTACa                          <  1:10004302/35‑1 (MQ=255)
          gcagcCTTACCAGTACGCTCACGCTGGTAGTACAg                         >  1:4032236/1‑35 (MQ=255)
          gcagcCTTACCAGTACGCTCACGCTGGTAGTACAg                         >  1:10383004/1‑35 (MQ=255)
           cagcCTTACCAGTACGCTCACGCTGGTAGTACAgc                        <  1:14702703/35‑2 (MQ=255)
           cagcCTTACCAGTACGCTCACGCTGGTAGTACAGt                        >  1:4741505/1‑35 (MQ=255)
           cagcCTTACCAGTACGCTCACGCTGGTAGTACAGt                        <  1:16820350/35‑1 (MQ=255)
            agcCTTACCAGTACGCTCACGCTGGTAGTACAGtt                       >  1:1947993/1‑35 (MQ=255)
            agcCTTACCAGTACGCTCACGCTGGTAGTACAGtt                       <  1:19206820/35‑1 (MQ=255)
            agcCTTACCAGTACGCTCACGCTGGTAGTACAGtt                       <  1:16950681/35‑1 (MQ=255)
            agcCTTACCAGTACGCTCACGCTGGTAGGACCGtt                       >  1:6007680/1‑35 (MQ=25)
            agcCTTACCAGTACGCTCACGCAGGTAGTACAGtt                       >  1:12473550/1‑35 (MQ=255)
             gcCTTACCAGTACGCTCACGCTGGTAGTACAGttt                      >  1:15220580/1‑35 (MQ=255)
             gcCTTACCAGTACGCTCACGCTGGTAGTACAGttt                      >  1:19071216/1‑35 (MQ=255)
              ccTTACCAGTACGCTCACGCTGGTAGTACAGTTTa                     >  1:11870616/1‑35 (MQ=255)
               cTTACCAGTACGCTCACGCTGGTAGTACAGTTTAg                    <  1:3317085/35‑1 (MQ=255)
               cTTACCAGTACGCTCACGCTGGTAGTACAGTTTAg                    >  1:12703485/1‑35 (MQ=255)
               cTTACCAGTACGCTCACGCTGGTAGTACAGTTTAg                    >  1:16210276/1‑35 (MQ=255)
               cTTACCAGTACGCTCACGCTGGTAGTACAGTTTAg                    >  1:2073263/1‑35 (MQ=255)
               cTTACCAGTACGCTCACGCTGGTAGTACAGTTTAg                    <  1:6858286/35‑1 (MQ=255)
               cTTACCAGTACGCTCACGCTGGTAGTACAGTTTAg                    >  1:9761550/1‑35 (MQ=255)
                ttACCAGTACGCTCACGCTGGTAGTACAGTTTAGc                   <  1:18601326/35‑1 (MQ=255)
                ttACCAGTACGCTCACGCTGGTAGTACAGTTTAGc                   <  1:6571734/35‑1 (MQ=255)
                ttACCAGTACGCTCACGCTGGTAGTACAGTTTAGc                   >  1:19178026/1‑35 (MQ=255)
                 tACCAGTACGCTCACGCTGGTAGTACAGTTTAGCt                  <  1:11279379/35‑1 (MQ=255)
                 tACCAGTACGCTCACGCTGGTAGTACAGTTTAGCt                  <  1:12327200/35‑1 (MQ=255)
                 tACCAGTACGCTCACGCTGGTAGTACAGTTTAGCt                  <  1:20322399/35‑1 (MQ=255)
                  aCCAGTACGCTCACGCTGGTAGTACAGTTTAGCtt                 <  1:16210867/35‑1 (MQ=255)
                  aCCAGTACGCTCACGCTGGTAGTACAGTTTAGCtt                 <  1:7016237/35‑1 (MQ=255)
                  aCCAGTACGCTCACGCTGGTAGTACAGTTTAGCtt                 >  1:5880398/1‑35 (MQ=255)
                  aCCAGTACGATCACGCTGGTAGTACAGTTTAGCtt                 <  1:10769071/35‑1 (MQ=37)
                   ccAGTACGCTCACGCTGGTAGTACAGTTTAGCttt                >  1:6088406/1‑35 (MQ=255)
                   ccAGTACGCTCACGCTGGTAGTACAGGTTAGCttt                >  1:19763183/1‑35 (MQ=37)
                   ccAGTACGCTCACGCTGGTAGGACAGTTTAGCttt                >  1:17888678/1‑35 (MQ=37)
                   ccAGTACGCTCACGCTCGTAGTACAGTTTAGCttt                >  1:3300996/1‑35 (MQ=37)
                    cAGTACGCTCACGCTGGTAGTACAGTTTAGCTTTa               >  1:13263099/1‑35 (MQ=255)
                    cAGTACGCTCACGCTGGTAGTACAGTTTAGCTTTa               <  1:13017303/35‑1 (MQ=255)
                    cAGTACGCTCACGCTGGTAGTACAGTTTAGCTTTa               <  1:3080493/35‑1 (MQ=255)
                     attACGCTCACGCTGGTAGTACAGTTTAGCTTTAc              >  1:7824470/3‑35 (MQ=38)
                     aGTACGCTCACGCTGGTAGTACAGTTTAGCTTTAc              >  1:5458709/1‑35 (MQ=255)
                     aGTACGCTCACGCTGGTAGTACAGTTTAGCTTTAc              <  1:17810108/35‑1 (MQ=255)
                     aGTACGCTCACGCTGGTAGTACAGTTTAGCTTTAc              <  1:2186630/35‑1 (MQ=255)
                      gTACGCTCACGCTTGTAGTACAGTTTAGCTTTACg             >  1:139228/1‑35 (MQ=38)
                      gTACGCTCACGCTGGTAGTACAGTTTAGCTTTACg             <  1:3213128/35‑1 (MQ=255)
                      gTACGCTCACGCTGGTAGTACAGTTTAGCTTTACg             >  1:13628570/1‑35 (MQ=255)
                       tACGCTCACGCTGGTAGTACAGTTTAGCTTTACGa            >  1:17246029/1‑35 (MQ=255)
                       tACGCTCACGCAGGTAGTACAGTTTAGCTTTACGa            >  1:11023201/1‑35 (MQ=255)
                        aCGCTCACGCTGGTAGTACAGTTTAGCTTTACGaa           <  1:8486635/35‑1 (MQ=255)
                        aCGCTCACGCTGGTAGTACAGTTTACCTTTACGaa           >  1:16466273/1‑35 (MQ=37)
                          gCTCACGCTGGTAGTACAGTTTAGCTTTACGAACa         <  1:9656235/35‑1 (MQ=255)
                          gCTCACGCTGGTAGTACAGTTTAGCTTTACGAACa         >  1:14334172/1‑35 (MQ=255)
                          gCTCACGCTGGTAGTACAGTTTAGCTTTACGAACa         >  1:18870179/1‑35 (MQ=255)
                          gCTCACGCTGGTAGTACAGTTTAGCTTTACGAACa         >  1:6613687/1‑35 (MQ=255)
                          gCTCACGCTGGTAGTACAGTTTAGCTTTACGAACa         >  1:16218465/1‑35 (MQ=255)
                          gCTCACGCTGGTAGGACAGTTTAGCTTTACGAACa         >  1:3117236/1‑35 (MQ=37)
                           cTCACGCTGGTAGTACAGTTTAGCTTTACGAACAg        <  1:1198612/35‑1 (MQ=255)
                           cTCACGCTGGTAGTACAGTTTAGCTTTACGAACAg        >  1:7379897/1‑35 (MQ=255)
                           cTCACGCTGGTAGTACAGTTTAGCTTTACGAACAg        <  1:18728168/35‑1 (MQ=255)
                            tCACGCTGGTAGTACAGTTTAGCTTTACGAACAGc       >  1:1662816/1‑35 (MQ=255)
                             cACGCTGGTAGTACAGTTTAGCTTTACGAACAGca      >  1:3193754/1‑35 (MQ=255)
                             cACGCTGGTAGTACAGTTTAGCTTTACGAACAGca      <  1:1456811/35‑1 (MQ=255)
                             cACGCTGGTAGTACAGTTTAGCTTTACGAACAGca      >  1:5037194/1‑35 (MQ=255)
                             cACGCTGGTAGTACAGTTTAGCTTTACGAACAGca      <  1:15709246/35‑1 (MQ=255)
                             cACGCAGGTAGTACAGTTTAGCTTTACGAACAGca      >  1:4733925/1‑35 (MQ=255)
                               cGCTGGTAGTACAGTTTAGCTTTACGAACAGcacc    >  1:8078493/1‑35 (MQ=255)
                               cGCTGGTAGTACAGTTTAGCTTTACGAACAGcacc    <  1:1404176/35‑1 (MQ=255)
                               cGCAGGTAGTACAGTTTAGCTTTACGAACAGcacc    >  1:2416179/1‑35 (MQ=255)
                                 cAGGTAGTACAGTTTAGCTTTACGAACANnaccan  >  1:10738184/1‑35 (MQ=255)
                                 cAGGTAGTACAGTTTAGCTTTACGAACAGcaccac  >  1:17503299/1‑35 (MQ=255)
                                  |                                 
TTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCAC  >  REL606/2665605‑2665672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ 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.