Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,472,392 G→T 6.4% pseudogene (3325/5622 nt) ebh ← very large surface anchored protein (pseudogene)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,472,3920GT 6.4% 1246.9 / 31.7 317pseudogene (3325/5622 nt)ebhvery large surface anchored protein (pseudogene)
Reads supporting (aligned to +/- strand):  ref base G (161/132);  new base T (4/20);  total (165/152)
Fisher's exact test for biased strand distribution p-value = 4.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

CTTGTTTCGCATTTTCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGATTTCGATC  >  BX571857/1472296‑1472488
                                                                                                |                                                                                                
cTTGTTTCGCATTTTCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAg                                                                                                  <  2:2491877/97‑1 (MQ=255)
   gTTTCGCATTTTCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAgcgc                                                                                               >  1:3119733/1‑97 (MQ=255)
      tCGCATTTTCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCAtt                                                                                            >  1:2111066/1‑97 (MQ=255)
      tCGCATTTTCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCAtt                                                                                            <  1:1541738/97‑1 (MQ=255)
          aTTTTCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTa                                                                                        >  2:2817277/1‑97 (MQ=255)
            tttCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAAt                                                                                      >  2:2025620/1‑97 (MQ=255)
               cTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTc                                                                                   <  2:4286579/97‑1 (MQ=255)
                 ttAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGt                                                                                 <  1:2373006/97‑1 (MQ=255)
                   aGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTa                                                                               <  1:2489591/97‑1 (MQ=255)
                   aGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTa                                                                               >  1:1659470/1‑97 (MQ=255)
                     cTTTCACTAAATCTTTAATGTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATg                                                                             <  2:3284402/97‑1 (MQ=255)
                      tttCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGa                                                                            >  2:1328378/1‑97 (MQ=255)
                      tttCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGa                                                                            >  1:3356002/1‑97 (MQ=255)
                       ttCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGAc                                                                           >  1:112334/1‑97 (MQ=255)
                         cACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCt                                                                         >  1:3516695/1‑97 (MQ=255)
                          aCTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCtt                                                                        >  1:4260835/1‑97 (MQ=255)
                           cTAAATCTTTGATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCttg                                                                       <  2:2056202/97‑1 (MQ=37)
                           cTAAATCTTTAATTTCTTGTAATGCTTGCGCAGGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCttg                                                                       <  2:957351/97‑1 (MQ=255)
                            tAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCttgt                                                                      >  2:1805810/1‑97 (MQ=255)
                               aTCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGa                                                                   >  1:3737907/1‑97 (MQ=255)
                                tCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGaa                                                                  >  2:3843828/1‑97 (MQ=255)
                                tCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGaa                                                                  >  1:719480/1‑97 (MQ=255)
                                tCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGaa                                                                  >  2:1175186/1‑97 (MQ=255)
                                 cTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAg                                                                 >  1:2059628/1‑97 (MQ=255)
                                  tttAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGt                                                                >  2:149774/1‑97 (MQ=255)
                                   ttAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTa                                                               >  2:3942043/1‑97 (MQ=255)
                                    tgaTTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTAt                                                              <  1:888391/95‑1 (MQ=255)
                                      aTTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTAttt                                                            >  1:1926453/1‑97 (MQ=255)
                                      aTTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTAttt                                                            >  1:4073155/1‑97 (MQ=255)
                                        ttCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGa                                                          >  1:569261/1‑97 (MQ=255)
                                         tCTTGTCATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGAt                                                         >  2:305691/1‑97 (MQ=255)
                                         tCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGAt                                                         >  1:1515125/1‑97 (MQ=255)
                                         tCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGAt                                                         >  1:2008048/1‑97 (MQ=255)
                                           ttGTAATGCTTGCGCAAGTTGTGCTTTGGCGTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTa                                                       >  1:1245901/1‑97 (MQ=255)
                                            tGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTaa                                                      <  1:2234104/97‑1 (MQ=255)
                                            tGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTaa                                                      <  2:1729262/97‑1 (MQ=255)
                                            tGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTaa                                                      <  1:2946106/97‑1 (MQ=255)
                                               aaTGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTa                                                               <  1:4106538/85‑1 (MQ=39)
                                                aTGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACg                                                  >  2:2570116/1‑97 (MQ=255)
                                                aTGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACg                                                  >  1:2803425/1‑97 (MQ=255)
                                                aTGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACg                                                  <  1:2774004/97‑1 (MQ=255)
                                                aTGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACg                                                  <  1:3174530/97‑1 (MQ=255)
                                                 tGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGa                                                 <  2:2719359/97‑1 (MQ=255)
                                                 tGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGa                                                 >  2:317710/1‑97 (MQ=255)
                                                  gCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGAt                                                >  1:1268498/1‑97 (MQ=255)
                                                  gCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGAt                                                >  2:2950662/1‑97 (MQ=255)
                                                   cTTGTGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATc                                               <  2:4113189/97‑1 (MQ=255)
                                                   cTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATc                                               >  1:104887/1‑97 (MQ=255)
                                                   cTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATc                                               >  2:3463194/1‑97 (MQ=255)
                                                   cTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATc                                               >  1:3600346/1‑97 (MQ=255)
                                                   cTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATAAGATc                                               >  1:1162938/1‑97 (MQ=255)
                                                     tGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCtt                                             >  2:225217/1‑97 (MQ=255)
                                                      gcgcAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGAcc                                                                          <  1:686577/67‑1 (MQ=40)
                                                      gcgcAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCttt                                            >  1:1856111/1‑97 (MQ=255)
                                                      gcgcAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCttt                                            >  1:3799197/1‑97 (MQ=255)
                                                       cgcAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTa                                           <  1:3991939/97‑1 (MQ=255)
                                                       cgcAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTa                                           >  1:1343654/1‑97 (MQ=255)
                                                        gcAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTaa                                          <  2:3520312/97‑1 (MQ=255)
                                                          aaGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCATCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGt                                        <  2:2367763/97‑1 (MQ=25)
                                                          aaGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGt                                        <  2:210677/97‑1 (MQ=255)
                                                           aGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTg                                       <  1:808094/97‑1 (MQ=255)
                                                            gTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGc                                      <  1:1484722/97‑1 (MQ=40)
                                                             ttGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCt                                     <  1:2784714/97‑1 (MQ=255)
                                                             ttGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCt                                     >  2:450944/1‑97 (MQ=255)
                                                             ttGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCt                                     >  2:106588/1‑97 (MQ=255)
                                                             ttGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCt                                     <  1:4357264/97‑1 (MQ=255)
                                                             ttGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCt                                     >  1:3439932/1‑97 (MQ=255)
                                                              tgtgCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCtt                                    >  2:3213025/1‑97 (MQ=255)
                                                               gtgCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTAt                                                              <  1:448504/70‑1 (MQ=35)
                                                                tgCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGt                                  >  2:830831/1‑97 (MQ=255)
                                                                  cTTTGGCTTGTTCAATTTCTTCTTTAGTCATCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGttt                                <  2:772512/97‑1 (MQ=21)
                                                                  cTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGttt                                >  1:2284087/1‑97 (MQ=255)
                                                                   tttGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGtttt                               >  2:786038/1‑97 (MQ=40)
                                                                    ttGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGttttt                              >  2:474248/1‑97 (MQ=255)
                                                                     tGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTaa                                                                                       >  2:4619508/1‑39 (MQ=255)
                                                                     tGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTaa                                                                                       <  1:4619508/39‑1 (MQ=255)
                                                                     tGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTc                             >  2:3666533/1‑97 (MQ=255)
                                                                       gCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCt                           <  1:3272004/97‑1 (MQ=255)
                                                                        cTTGTTCTATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCtt                          <  2:4429633/97‑1 (MQ=255)
                                                                         ttGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTa                         >  1:783234/1‑97 (MQ=255)
                                                                          tGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAAt                                                                                      <  2:4321819/35‑1 (MQ=255)
                                                                           gTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATc                       <  2:340053/97‑1 (MQ=40)
                                                                            ttCAATTTCTTCTTTAGTCATCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGt                                        <  2:1628952/79‑1 (MQ=255)
                                                                             tCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTg                     <  1:790801/97‑1 (MQ=255)
                                                                              cAATTTCTTCTTTAGTCATCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGAt                                                >  2:4035123/1‑69 (MQ=255)
                                                                              cAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGt                    >  1:4583132/1‑97 (MQ=255)
                                                                              cAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGt                    >  1:3196771/1‑97 (MQ=255)
                                                                                 tttCTTCTTTAGTCATCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAg                 <  2:2682108/97‑1 (MQ=21)
                                                                                 tttCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAg                 >  2:679793/1‑97 (MQ=255)
                                                                                     ttcttTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATGCGATCTTTAAGTGCTTGTTTTTCCTTATCt                      <  1:3541790/88‑1 (MQ=255)
                                                                                      tcttTAGTCATCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTg            <  2:3737104/97‑1 (MQ=21)
                                                                                       cttTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGt                    <  1:1400581/88‑1 (MQ=255)
                                                                                       cttTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTgg           <  2:1950977/97‑1 (MQ=255)
                                                                                        tttAGTCATCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGa          <  2:3026075/97‑1 (MQ=18)
                                                                                         ttAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGAt         <  1:1872623/97‑1 (MQ=255)
                                                                                          tAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGAtt        >  1:2701458/1‑97 (MQ=255)
                                                                                            gTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGATTTc      <  2:130147/97‑1 (MQ=255)
                                                                                             tCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGATTTCg     <  1:3035549/97‑1 (MQ=255)
                                                                                              cAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGATTTCGa    >  2:3152811/1‑97 (MQ=255)
                                                                                               aGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGATTTCGAt   >  1:1132793/1‑97 (MQ=255)
                                                                                                |                                                                                                
CTTGTTTCGCATTTTCTTTAGCTTTCACTAAATCTTTAATTTCTTGTAATGCTTGCGCAAGTTGTGCTTTGGCTTGTTCAATTTCTTCTTTAGTCAGCGCATTGTTAATGTCGTTATGACCTTGTTGAAGTATTTGATTAATACGATCTTTAAGTGCTTGTTTTTCCTTATCTGTTAGATTTGGATTTCGATC  >  BX571857/1472296‑1472488

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