Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 621,623 G→A 69.4% A42A (GCG→GCA SAS0539 → putative uracil‑DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857621,6230GA69.4% 429.5 / 251.4 232A42A (GCG→GCASAS0539putative uracil‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (42/29);  new base A (97/64);  total (139/93)
Fisher's exact test for biased strand distribution p-value = 8.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 233 total aligned reads displayed.

AACAACGAAACATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCATTTTCAGT  >  BX571857/621527‑621718
                                                                                                |                                                                                               
aacaacGAAACATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCg                                                                                                 <  1:1940999/97‑1 (MQ=255)
   aacGAAACATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGttt                                                                                              >  2:509572/1‑97 (MQ=255)
   aacGAAACATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCAttt                                                                                              <  1:987593/97‑1 (MQ=255)
      gAAACATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGAt                                                                                           >  2:3220699/1‑97 (MQ=255)
      gAAACATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGAt                                                                                           >  1:3127505/1‑97 (MQ=255)
       aaaCATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGAtt                                                                                          <  1:2554826/97‑1 (MQ=255)
       aaaCATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGAtt                                                                                          >  1:3242345/1‑97 (MQ=255)
        aaCATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGAttt                                                                                         <  2:28128/97‑1 (MQ=255)
         aCATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTa                                                                                        >  2:543458/1‑97 (MQ=255)
         aCATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTa                                                                                        >  2:370480/1‑97 (MQ=255)
          cATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTaa                                                                                       >  1:1866342/1‑97 (MQ=255)
          cATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTaa                                                                                       >  1:1147652/1‑97 (MQ=255)
          cATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTaa                                                                                       >  1:2150078/1‑97 (MQ=255)
           aTGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAac                                                                                      <  2:1679907/97‑1 (MQ=255)
           aTGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGAATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAac                                                                                      >  1:1744060/1‑97 (MQ=255)
              aCTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACAcc                                                                                   >  1:2371410/1‑97 (MQ=255)
                tttAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGt                                                                                 >  1:1985977/1‑97 (MQ=255)
                    aaGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGa                                                                             <  2:471258/97‑1 (MQ=255)
                      gCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGaaa                                                                           <  1:2030327/97‑1 (MQ=255)
                       cTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGaaaa                                                                          >  2:39300/1‑97 (MQ=255)
                        tATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAat                                                                         <  2:1767367/97‑1 (MQ=255)
                         aTGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAata                                                                        <  2:3354620/97‑1 (MQ=255)
                          tGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAatat                                                                       >  1:2255112/1‑97 (MQ=255)
                            cATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCa                                                                     >  1:2405463/1‑97 (MQ=255)
                            cATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCa                                                                     >  1:1294723/1‑97 (MQ=255)
                            cATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCa                                                                     >  2:149088/1‑97 (MQ=255)
                             aTGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCaa                                                                    <  1:2617592/97‑1 (MQ=255)
                             aTGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCaa                                                                    <  2:3255117/97‑1 (MQ=255)
                              tGATTTTTTAGAAAAAGAACATTAGACTGCAATCGTATACCCGGAGAGGGCAAATATATATCAAGCAGATGATTTAACACCGTTTGAAACTATCaaa                                                                   >  2:35806/1‑97 (MQ=255)
                               gATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAg                                                                  <  2:2089008/97‑1 (MQ=255)
                                aTTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAgt                                                                 >  1:1376518/1‑97 (MQ=255)
                                 tttttaAGAAAAAGAATATTCGGCTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGGTTGAAAAAATCAAAgtt                                                                <  2:1440306/97‑1 (MQ=255)
                                  tttttAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAgttg                                                               <  1:2329453/97‑1 (MQ=255)
                                   ttttAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAgttgt                                                              <  2:2255252/97‑1 (MQ=255)
                                   ttttAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAgttgt                                                              >  2:156102/1‑97 (MQ=255)
                                   ttttAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAgttgt                                                              >  1:1658436/1‑97 (MQ=255)
                                     ttAGAATAAGAATATTCGACGGCAATCGTATACCCTGATAGGGAAAATATATATCACGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTta                                                            <  2:2069823/97‑1 (MQ=255)
                                     ttAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTta                                                            >  1:2867391/1‑97 (MQ=255)
                                     ttAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTta                                                            >  2:1605615/1‑97 (MQ=255)
                                         aaaaaGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATAtt                                                        <  1:694802/97‑1 (MQ=255)
                                           aaaGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAg                                                      <  2:2535158/97‑1 (MQ=255)
                                            aaGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAgg                                                     >  2:2104891/1‑97 (MQ=255)
                                             aGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGa                                                    >  1:248357/1‑97 (MQ=255)
                                             aGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGa                                                    >  2:1906248/1‑97 (MQ=255)
                                              gAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGAc                                                   >  2:3314969/1‑97 (MQ=255)
                                              gAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGAc                                                   >  1:2634299/1‑97 (MQ=255)
                                               aaTATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACa                                                  >  2:592616/1‑97 (MQ=255)
                                                atatTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGATGTTATATTAGGACaa                                                 >  1:1010593/1‑97 (MQ=255)
                                                 tatTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAg                                                >  2:223476/1‑97 (MQ=255)
                                                 tatTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAg                                                >  1:358094/1‑97 (MQ=255)
                                                 tatTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAg                                                >  2:2659080/1‑97 (MQ=255)
                                                 tatTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAg                                                <  2:1173888/97‑1 (MQ=255)
                                                  atTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGa                                               >  1:65201/1‑97 (MQ=255)
                                                   ttCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAc                                              <  1:1979160/97‑1 (MQ=255)
                                                   ttCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAc                                              >  1:2074831/1‑97 (MQ=255)
                                                    tCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAcc                                             <  1:1634875/97‑1 (MQ=255)
                                                    tCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAcc                                             >  2:1227838/1‑97 (MQ=255)
                                                    tCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAcc                                             <  1:1375848/97‑1 (MQ=255)
                                                     cGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAccc                                            >  1:1042154/1‑97 (MQ=255)
                                                     cGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGCTTCATTCAACACCGTCTGAAAATATCAAAGTCGTTATATTAGGACAAGAccc                                            >  2:892582/1‑97 (MQ=255)
                                                     cGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAccc                                            >  2:2396380/1‑97 (MQ=255)
                                                     cGACTACAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGAccc                                            >  2:2452091/1‑97 (MQ=255)
                                                      gACTGCAATCGTAGACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCg                                           <  2:1965955/97‑1 (MQ=255)
                                                       aCTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGt                                          <  1:1752430/97‑1 (MQ=255)
                                                       aCTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGt                                          <  1:1076981/97‑1 (MQ=255)
                                                        cTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTa                                         <  1:539864/97‑1 (MQ=255)
                                                         tGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTAt                                        >  2:2953912/1‑97 (MQ=255)
                                                          gCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATc                                       >  2:806572/1‑97 (MQ=255)
                                                          gCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATc                                       <  1:3177503/97‑1 (MQ=255)
                                                          gCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATc                                       <  1:794276/97‑1 (MQ=255)
                                                          gCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATc                                       >  1:383567/1‑97 (MQ=255)
                                                          gCAATCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATc                                       >  1:2582438/1‑97 (MQ=255)
                                                            aaTCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCAt                                     <  2:440845/97‑1 (MQ=255)
                                                              tCGTATACCCTGATAGGGAACATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATgg                                   >  2:3149623/1‑97 (MQ=255)
                                                              tCGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATgg                                   <  2:649239/97‑1 (MQ=255)
                                                               cGTATACCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGt                                  >  1:2212857/1‑97 (MQ=255)
                                                                   taCCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCaa                              >  2:1351258/1‑97 (MQ=255)
                                                                   taCCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCaa                              >  1:3013447/1‑97 (MQ=255)
                                                                   taCCCTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCaa                              >  1:2198618/1‑97 (MQ=255)
                                                                     cccTGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAAc                            >  1:3298015/1‑97 (MQ=255)
                                                                        tGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCaa                         <  2:2883431/97‑1 (MQ=255)
                                                                        tGATAGGGAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCaa                         >  1:3326074/1‑97 (MQ=255)
                                                                             gggAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGcaca                    <  2:3206721/97‑1 (MQ=255)
                                                                             gggAAAATATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGcaca                    >  1:844275/1‑97 (MQ=255)
                                                                               gAAAATATATATCAAGCATTTGATTTAACTCCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATg                  <  1:1263634/97‑1 (MQ=255)
                                                                                aaaaTATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGGTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATgg                 <  2:2716817/97‑1 (MQ=255)
                                                                                  aaTATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTAACATGGTCCAAACCAAGCACATGGAt               >  1:1338175/1‑97 (MQ=255)
                                                                                  aaTATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGAt               <  1:3187923/97‑1 (MQ=255)
                                                                                  aaTATATATCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGAt               >  2:2771186/1‑97 (MQ=255)
                                                                                   atatatatCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGAtt              >  2:1257294/1‑97 (MQ=255)
                                                                                    tatatatCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACCTGGATTa             >  2:2878851/1‑97 (MQ=255)
                                                                                     atatatCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAg            >  1:1563588/1‑97 (MQ=255)
                                                                                       atatCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCa          <  2:73958/97‑1 (MQ=255)
                                                                                       atatCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCa          <  1:3009226/97‑1 (MQ=255)
                                                                                        tatCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCAt         <  2:2388883/97‑1 (MQ=255)
                                                                                         atCAAGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCAtt        >  1:1232271/1‑97 (MQ=255)
                                                                                            aaGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCATTTTc     >  1:87502/1‑97 (MQ=255)
                                                                                            aaGCATTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCATTTTc     <  2:1325211/97‑1 (MQ=255)
                                                                                              gCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCATTTTctg   >  2:257248/1‑95 (MQ=255)
                                                                                               cGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCATTTTCAGt  >  1:577101/1‑97 (MQ=255)
                                                                                                |                                                                                               
AACAACGAAACATGACTTTAAAGCTATGCATGATTTTTTAGAAAAAGAATATTCGACTGCAATCGTATACCCTGATAGGGAAAATATATATCAAGCGTTTGATTTAACACCGTTTGAAAATATCAAAGTTGTTATATTAGGACAAGACCCGTATCATGGTCCAAACCAAGCACATGGATTAGCATTTTCAGT  >  BX571857/621527‑621718

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