Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 868,946 G→A 22.3% A315T (GCT→ACT)  SAS0786 → putative selenocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857868,9460GA22.3% 494.6 / 152.5 202A315T (GCT→ACT) SAS0786putative selenocysteine lyase
Reads supporting (aligned to +/- strand):  ref base G (80/77);  new base A (26/19);  total (106/96)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 206 total aligned reads displayed.

ATTGCTCAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACAAGATGTAC  >  BX571857/868850‑869042
                                                                                                |                                                                                                
aTTGCTCAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGCTGCAATTCATAAATATGAACAAGAATTAACGATATATg                                                                                                  >  2:601556/1‑97 (MQ=255)
  tGCTCAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCt                                                                                                >  1:15650/1‑97 (MQ=255)
   gCTCAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCtt                                                                                               >  1:1466582/1‑97 (MQ=255)
   gCTCAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCtt                                                                                               >  1:1388497/1‑97 (MQ=255)
     tCAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTAt                                                                                             <  1:1783530/97‑1 (MQ=255)
      cAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATg                                                                                            <  2:337366/97‑1 (MQ=255)
      cAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATg                                                                                            <  1:949687/97‑1 (MQ=255)
      cAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATg                                                                                            >  1:1590269/1‑97 (MQ=255)
      cAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATg                                                                                            <  1:126536/97‑1 (MQ=255)
        aGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAg                                                                                          >  2:190446/1‑97 (MQ=255)
         gCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGc                                                                                         <  2:1762313/97‑1 (MQ=255)
          cAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCa                                                                                        <  2:871530/97‑1 (MQ=255)
          cAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCa                                                                                        <  2:272789/97‑1 (MQ=255)
           aaTTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCaa                                                                                       <  2:1810216/97‑1 (MQ=255)
             ttGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAAt                                                                                     >  2:1440837/1‑97 (MQ=255)
              tGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATg                                                                                    <  1:2072570/97‑1 (MQ=255)
              tGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATg                                                                                    <  1:1983578/97‑1 (MQ=255)
              ggggCTTGCAGAAGCTCTTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATg                                                                                    <  1:2117856/96‑1 (MQ=255)
               gggCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGATTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGt                                                                                   >  2:2218354/1‑97 (MQ=255)
                 gCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCt                                                                                 <  2:1491738/97‑1 (MQ=255)
                 gCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCt                                                                                 >  1:2416482/1‑97 (MQ=255)
                   ttGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGc                                                                               <  1:1182988/97‑1 (MQ=255)
                     gCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCaa                                                                             >  1:1429049/1‑97 (MQ=255)
                       aGAAGCTATTCGGGATTGAGGACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAAtt                                                                           <  2:1608262/97‑1 (MQ=255)
                         aaGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGa                                                                         <  1:1971729/97‑1 (MQ=255)
                          aGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGaa                                                                        >  2:955317/1‑97 (MQ=255)
                           gCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAg                                                                       <  2:287619/97‑1 (MQ=255)
                            cTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAgg                                                                      >  2:675454/1‑97 (MQ=255)
                             tATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGa                                                                     >  2:831111/1‑97 (MQ=255)
                                  gCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGa                                                                <  2:896229/97‑1 (MQ=255)
                                   cTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGaa                                                               <  1:1500541/97‑1 (MQ=255)
                                   cTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGaa                                                               <  1:2533566/97‑1 (MQ=255)
                                      tttAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAAtt                                                            >  1:862644/1‑97 (MQ=255)
                                        tAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTa                                                          >  2:236899/1‑97 (MQ=255)
                                        tAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTa                                                          <  1:1642478/97‑1 (MQ=255)
                                        tAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTa                                                          >  2:264751/1‑97 (MQ=255)
                                                 tAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGCTATATACTTATGAGCAAATGTCTGCAATTGAAGGAATGGAAATTTATGGCCCGcc                                                 >  2:1227081/1‑97 (MQ=255)
                                                 tAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGcc                                                 >  1:1439459/1‑97 (MQ=255)
                                                  aGGGTTTGAGGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCa                                                <  2:1879657/97‑1 (MQ=255)
                                                   ggTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCaa                                               >  2:323081/1‑97 (MQ=255)
                                                   ggTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCaa                                               >  1:433451/1‑97 (MQ=255)
                                                     ttttGATGCAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAg                                             >  1:88298/1‑97 (MQ=255)
                                                      tttGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAgg                                            <  1:112338/97‑1 (MQ=255)
                                                       ttGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCACTTGAAGGAATTGAAATTTATGGCCCGCCAAAGGa                                           >  2:1966127/1‑97 (MQ=255)
                                                        tGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAt                                          <  2:504774/97‑1 (MQ=255)
                                                         gATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtc                                         <  1:1241399/97‑1 (MQ=255)
                                                         gATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtc                                         <  2:1391673/97‑1 (MQ=255)
                                                           tGCAATTGATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtcgt                                       <  2:1063546/97‑1 (MQ=255)
                                                           tGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtcgt                                       <  2:2510599/97‑1 (MQ=255)
                                                            gCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtcgtc                                      >  1:738951/1‑97 (MQ=255)
                                                            gCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtcgtc                                      >  1:1993088/1‑97 (MQ=255)
                                                             cAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtcgtcg                                     >  2:342784/1‑97 (MQ=255)
                                                             cAATTCATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGAtcgtcg                                     >  2:1051258/1‑97 (MQ=255)
                                                              aaTTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATcgtcgt                                    >  1:676022/1‑97 (MQ=255)
                                                              aaTTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAACTGAAGGAAATGAAATTTATGGCCCGCCAAAGGATcgtcgt                                    >  1:759048/1‑97 (MQ=255)
                                                               aTTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTg                                   <  2:1309346/97‑1 (MQ=255)
                                                                  cATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAg                                >  2:1252175/1‑97 (MQ=255)
                                                                  cATAAATATGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAg                                >  2:1528769/1‑97 (MQ=255)
                                                                   ggaaaTATGACGAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAgg                               <  2:2179393/95‑1 (MQ=255)
                                                                   aTAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAgg                               >  1:2480067/1‑97 (MQ=255)
                                                                   aTAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAgg                               <  1:341857/97‑1 (MQ=255)
                                                                   aTAAAGAGGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAgg                               <  2:745410/97‑1 (MQ=255)
                                                                    tAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGgt                              <  2:1239852/97‑1 (MQ=255)
                                                                     aaaTATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGgtg                             >  1:1415329/1‑97 (MQ=255)
                                                                     aaaTATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGgtg                             >  1:1875493/1‑97 (MQ=255)
                                                                     aaaTATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGgtg                             >  1:363100/1‑97 (MQ=255)
                                                                      aaTATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGgtgt                            <  1:379846/97‑1 (MQ=255)
                                                                      aaTATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGgtgt                            >  1:128249/1‑97 (MQ=255)
                                                                        tatGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGtaa                          >  1:1204987/1‑97 (MQ=255)
                                                                         atGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGtaat                         <  1:2434379/97‑1 (MQ=255)
                                                                         atGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGtaat                         >  1:2470921/1‑97 (MQ=255)
                                                                         atGAACAAGAATTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGtaat                         >  1:158017/1‑97 (MQ=255)
                                                                         atGAACAAGAACTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGtaat                         <  2:1280958/97‑1 (MQ=255)
                                                                          tGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGtaata                        >  1:2095319/1‑97 (MQ=255)
                                                                           gAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGtaataa                       >  1:472971/1‑97 (MQ=255)
                                                                            aaCAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAAc                      <  1:2221206/97‑1 (MQ=255)
                                                                            aaCAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAAc                      <  2:2055376/97‑1 (MQ=255)
                                                                              cAAGAATTAACGATATATGCTTATGAGCACATGACTGCAATTGAAGGACTTGAACTTTATGGCCCGCCAAAGGACCGTCGGGCCGGCGTAATAACGt                    >  1:620425/1‑97 (MQ=255)
                                                                              cAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGt                    >  2:2403775/1‑97 (MQ=255)
                                                                              cAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGt                    >  1:275663/1‑97 (MQ=255)
                                                                               aaGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGtt                   <  1:527020/97‑1 (MQ=255)
                                                                               aaGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGtt                   >  1:1297891/1‑97 (MQ=255)
                                                                                aGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGttt                  >  1:1953450/1‑97 (MQ=255)
                                                                                 gAATTAGCGATATATGCTTATGGGCAAATGGCTGCAATGGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTa                 <  2:2458794/97‑1 (MQ=255)
                                                                                   aTTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAAt               >  1:1526404/1‑97 (MQ=255)
                                                                                   aTTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAAt               >  1:132227/1‑97 (MQ=255)
                                                                                   aTTAACGATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAAt               >  1:2382450/1‑97 (MQ=255)
                                                                                    ttAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAAtt              >  2:546695/1‑97 (MQ=255)
                                                                                    ttAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAAtt              >  2:980001/1‑97 (MQ=255)
                                                                                    ttAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAAtt              >  1:820414/1‑97 (MQ=255)
                                                                                      aaCGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTa            <  2:2185026/97‑1 (MQ=255)
                                                                                      aaCGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTa            >  1:2241752/1‑97 (MQ=255)
                                                                                         gATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACaa         >  1:504145/1‑97 (MQ=255)
                                                                                         gATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACaa         >  1:2309884/1‑97 (MQ=255)
                                                                                         gATATATACTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACaa         <  2:2332099/97‑1 (MQ=255)
                                                                                           tatatGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACAAgg       >  2:461351/1‑96 (MQ=255)
                                                                                            atatGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACAAGAt      >  1:1037374/1‑97 (MQ=255)
                                                                                               tGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACAAGATGTa   <  2:1815497/97‑1 (MQ=255)
                                                                                               tGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGAGTAATAACGTTTAATTTACAAGATGTa   >  2:2274862/1‑97 (MQ=255)
                                                                                                gCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACAAGATGTAc  >  2:785902/1‑97 (MQ=255)
                                                                                                |                                                                                                
ATTGCTCAAGCAATTGGGCTTGCAGAAGCTATTCGCTATTTAGAACGCATAGGTTTTGATGCAATTCATAAATATGAACAAGAATTAACGATATATGCTTATGAGCAAATGTCTGCAATTGAAGGAATTGAAATTTATGGCCCGCCAAAGGATCGTCGTGCAGGTGTAATAACGTTTAATTTACAAGATGTAC  >  BX571857/868850‑869042

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