Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 953,823 T→G 17.6% K203T (AAG→ACG)  selU ← tRNA 2‑selenouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122953,8230TG17.6% 132.5 / 4.7 72K203T (AAG→ACG) selUtRNA 2‑selenouridine synthase
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base G (13/2);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 8.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTGGTGAAGCTGTCTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTCCAGCTGGTGCAGC  >  ACB122/953688‑953957
                                                                                                                                       |                                                                                                                                      
cgGTTGGTGAAGCTGTCTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCt                                                                                                                                        >  7:212060/1‑136 (MQ=255)
  gTTGGTGAAGCTGTCTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCtt                                                                                                                                       <  1:79193/135‑1 (MQ=255)
         aaGCTGTCTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCacc                                                                                                                                >  6:4062/1‑133 (MQ=255)
            cTGTCTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGt                                                                                                                            <  8:102543/136‑1 (MQ=255)
             tGTCTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGt                                                                                                                            <  6:304085/135‑1 (MQ=255)
                cTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGCCAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGa                                                                                                                         >  2:170095/1‑135 (MQ=255)
                cTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGGCGa                                                                                                                         >  6:46795/1‑135 (MQ=255)
                            ccAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGt                                                                                                             >  4:261327/1‑135 (MQ=255)
                            ccAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGt                                                                                                             >  6:72758/1‑135 (MQ=255)
                                   gTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGAATGCCAGCTGGTTTTCaaa                                                                                                      >  6:208222/1‑135 (MQ=255)
                                         ctgcATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGAt                                                                                                <  7:116678/135‑1 (MQ=255)
                                          tgcATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGCACTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGGCGATTGCCAGCTGGTTTTGCAAGTCGctct                                                                                              >  8:8980/1‑136 (MQ=255)
                                            cATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGGGTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTAGATCTg                                                                                             >  6:32022/1‑135 (MQ=255)
                                               gCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGCTCTgcgc                                                                                          >  1:59078/1‑135 (MQ=255)
                                               gCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTgcgc                                                                                          >  2:250353/1‑135 (MQ=255)
                                               gCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTGAGCACGGCGCTTGGCAGCGGGGTTTGAAAGTCGCTCTgggc                                                                                          >  5:272175/1‑135 (MQ=255)
                                               gCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGTCGATTGCCAGCTGGGTTTCAAAGTCGATCTGCGcc                                                                                         >  8:383102/1‑136 (MQ=255)
                                                    ggTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCCCGTCGAGTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTg                                                                                    >  8:290638/1‑136 (MQ=255)
                                                       aCAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAACTCGATCTGCGCCGGTTgcg                                                                                  >  3:167485/1‑135 (MQ=255)
                                                         aaCTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTGATCACGTCGATTGGCAGCTGGTTTTGTAAGTCGTTCTGCGCCGGTTgcggc                                                                                >  1:139444/1‑133 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGGGTTTTAAGCACGTCGCTTTGCAGCAGGGGTTTCAAGTAGGTATGGGCGCGTTGCGCCGGGgc                                                                           >  5:219219/1‑135 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTGAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGgcg                                                                          >  8:19457/1‑136 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCCCGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGgc                                                                           >  4:322741/1‑135 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGgcg                                                                          >  7:45489/1‑136 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGgcg                                                                          >  8:203176/1‑136 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGgcg                                                                          >  8:264569/1‑136 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTAAGCACGGCGATTTCCAGCTGGTTTTCAAATTCGATATGCGCCGGTTGCGCCGTGgc                                                                           >  4:263809/1‑135 (MQ=255)
                                                              cAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGTCGATTGCCAGCTGGTTTTGAAAGTCGATCTGCGCCGGTTGCGCCGTGgc                                                                           >  3:233549/1‑135 (MQ=255)
                                                               aGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTTGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGgc                                                                           <  6:336207/134‑1 (MQ=255)
                                                                   gCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGACCGAACTGCCCGATACAGCGGGCGCGGTTTTTCAGCCCGTCGATTGCCAGCTGGTTTTCAAACTCCATCAGCGCCCGTTGCGCCCTGGCCGGCtt                                                                     >  7:111254/1‑136 (MQ=255)
                                                                   gCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCCGCACGTCGATTGCCCGCTGGTTTTCAAAGTGGATCTGCGCCGGTTGCGCCGTGGCGCGCtt                                                                     >  7:61478/1‑136 (MQ=255)
                                                                   gCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCtt                                                                     <  8:86021/136‑1 (MQ=255)
                                                                   gCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCt                                                                      >  5:155994/1‑135 (MQ=255)
                                                                        gTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCGTGcc                                                                  >  1:45089/1‑134 (MQ=255)
                                                                        gTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGGTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCgc                                                                >  8:206929/1‑135 (MQ=255)
                                                                        gTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGGTTTCAGCACGGCGAGTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCGTGCCg                                                                 >  1:61810/1‑135 (MQ=255)
                                                                        gTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCCCGTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGTGCGCTTGCCg                                                                 >  7:5113/1‑135 (MQ=255)
                                                                          gCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGaa                                                               <  3:188462/135‑1 (MQ=255)
                                                                           cAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTGAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGaaa                                                              >  2:63144/1‑135 (MQ=255)
                                                                           cAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTGAGCACGGCGATTGGCAGCGGGTTTTCAAAGTCGAGCTGCGCCGGTTGCGCCGTGGCGCGCGTGCCGaaa                                                              >  2:92585/1‑135 (MQ=255)
                                                                           cAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGTCGATTGCCCGCTGGGTGTCAAAGTCGCGCAGAGGCGGTTGCGCCGTGGCGCGCGTGCCgcc                                                               >  6:96914/1‑132 (MQ=255)
                                                                           cAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGTAGAATGCCAGCGGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGaaa                                                              >  7:177889/1‑135 (MQ=255)
                                                                              cTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTACGCACGTCGATTGCCAGCTGGTTTTCAAAGTAGAGCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTg                                                           >  3:166113/1‑135 (MQ=255)
                                                                                  cGACAATGCGGCCTTCATCTTCCAGAACTAACTGCCAGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCTGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGc                                                       <  6:217268/135‑1 (MQ=255)
                                                                                   gACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGcc                                                      <  1:89314/135‑1 (MQ=255)
                                                                                           ggCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGGGTTTTCAGCACCGCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGGGTGCCGACACTGGCGCCGCGCTggg                                              >  5:367531/1‑134 (MQ=255)
                                                                                           ggCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCATTTTAAGCCCGTCGATTTCCAGCAGGGTTTCAAAGTCGATCTGCGCCCGTTGCGCCGTGGCGCGCTTGCCCAAACTGGAGCCGCCAAGGt                                              >  6:283792/1‑135 (MQ=255)
                                                                                              cTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCCCGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTgg                                           >  1:1005/1‑135 (MQ=255)
                                                                                              cTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTgg                                           <  6:54659/135‑1 (MQ=255)
                                                                                              cTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGc                                          <  8:70816/136‑1 (MQ=255)
                                                                                               ttCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGc                                          <  3:202149/135‑1 (MQ=255)
                                                                                               ttCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGc                                          <  3:263809/135‑1 (MQ=255)
                                                                                               ttCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGc                                          <  5:72758/135‑1 (MQ=255)
                                                                                               ttCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGc                                          <  2:384845/135‑1 (MQ=255)
                                                                                                  aTCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGt                                        <  1:250353/134‑1 (MQ=255)
                                                                                                       ccAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCt                                   <  8:134559/134‑1 (MQ=255)
                                                                                                        cAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCt                                   <  6:155994/133‑1 (MQ=255)
                                                                                                        cAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCt                                   <  5:352738/133‑1 (MQ=255)
                                                                                                        cAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCt                                   <  5:32022/133‑1 (MQ=255)
                                                                                                        cAGAAAGACCTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCAAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGATTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCt                                   <  1:175221/133‑1 (MQ=255)
                                                                                                          gaacgaacTGCCCGATACCGCGGGCGCGCGTTTTCAGCACGTAGATTGCCAGCTGGTTTTCAAAGTCGAACTTCGCCCGTTTCGCCGTGGCGCGCATGCCGAAACTGGAGCCGCGATGGTTTGCGTTGGCttaca                               >  2:104513/1‑132 (MQ=255)
                                                                                                            acgaacTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAgg                             <  1:28762/135‑1 (MQ=255)
                                                                                                               aacTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTcc                          <  5:4062/135‑1 (MQ=255)
                                                                                                                  tGCCCGATACCGCGGGCGCGCTTTTTCCGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGc                       <  5:329587/135‑1 (MQ=255)
                                                                                                                     ccGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCAGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCa                      <  5:208222/133‑1 (MQ=255)
                                                                                                                           ccGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTcc              <  7:203176/135‑1 (MQ=255)
                                                                                                                           ccGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTcc              <  2:139610/135‑1 (MQ=255)
                                                                                                                            cgcgGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGGTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTcc              <  4:363183/134‑1 (MQ=255)
                                                                                                                                 gcgcgcTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCtt                                                                     <  1:223417/74‑1 (MQ=255)
                                                                                                                                 gcgcgcTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCtt                                                                     >  2:223417/1‑74 (MQ=255)
                                                                                                                                 gcgcgcTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATATGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGCGCCGCGCTGGTTGGCGTGGCCGTCAAGGTCCCGCACGTTGTACAGCTgg        >  4:159844/1‑135 (MQ=255)
                                                                                                                                  cgcgcTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTCCAGCTGGt       <  3:296438/135‑1 (MQ=255)
                                                                                                                                  cgcgcTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGGTGGCGTGGCCTTCAAGGTCCAGCACGTTGTCCAGCTGGt       <  6:219219/135‑1 (MQ=255)
                                                                                                                                      cTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTCCAGCTGGTgcagc  <  7:264569/136‑1 (MQ=255)
                                                                                                                                       tttttCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGAGCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTCCAGCTGGTgcagc  <  2:115972/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
CGGTTGGTGAAGCTGTCTTCAAGCCACACCAGCGGGTAATGCTGCATGCCTTGGTACAACTCCAGCGGCACCGTGCAGCTGCCGACAATGCGGCCTTCATCTTCCAGAACGAACTGCCCGATACCGCGGGCGCGCTTTTTCAGCACGTCGATTGCCAGCTGGTTTTCAAAGTCGATCTGCGCCGGTTGCGCCGTGGCGCGCTTGCCGAAACTGGAGCCGCGATGGTTGGCGTGGCCTTCAAGGTCCAGCACGTTGTCCAGCTGGTGCAGC  >  ACB122/953688‑953957

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