Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,434,365 G→A 13.2% G144S (GGT→AGT)  SAS1338 → putative tetrahydrodipicolinate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,434,3650GA13.2% 679.4 / 29.9 219G144S (GGT→AGT) SAS1338putative tetrahydrodipicolinate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (87/101);  new base A (29/2);  total (116/103)
Fisher's exact test for biased strand distribution p-value = 2.92e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

TGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATCGG  >  BX571857/1434268‑1434462
                                                                                                 |                                                                                                 
tGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTg                                                                                                   >  1:352111/1‑98 (MQ=255)
tGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACAACTg                                                                                                   >  2:1336391/1‑98 (MQ=255)
 gTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACCATTgg                                                                                                  >  2:506724/1‑98 (MQ=255)
  tCGTTGGGATGGTCGCAGCTGTTAATGTTGGCGCGGTCGGTGGCGAAGGTACAATGATTGATATGGATGCTGCTCGCGGTGGTCGTGCTACAACTGGt                                                                                                 <  2:1360796/98‑1 (MQ=255)
     ttATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                                              <  1:3114787/98‑1 (MQ=255)
      tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                                             >  1:3503600/1‑98 (MQ=255)
       atgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                                            <  1:384656/98‑1 (MQ=255)
       atgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                                            <  1:1896740/98‑1 (MQ=255)
       atgatgCGCGTAACGATGAAGGTTGGCGCGGTCGGTGGCGAAGGTACAATGATTGATATGAATGCTTCTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                                            <  2:677680/98‑1 (MQ=255)
         gatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                          >  1:2212115/1‑98 (MQ=255)
           tgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTa                                                                                        <  2:792998/98‑1 (MQ=255)
            gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                                       >  1:3053314/1‑98 (MQ=255)
            gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCATTAGAAAAAAAGTAc                                                                                       >  2:2390605/1‑98 (MQ=255)
              gcgcAACAATTAATGTTGGTGCTGCTGTTGGCGAAAGTACAACGATTGATATGAATGCTGCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACAt                                                                                     <  2:124491/98‑1 (MQ=255)
              gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTGGTAAAAATGTACAt                                                                                     >  2:198749/1‑98 (MQ=255)
              gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACCACTGGAAAAAATGTACAt                                                                                     >  2:2603905/1‑98 (MQ=255)
               cgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                                    >  1:2870505/1‑98 (MQ=255)
                  aacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTAGAAAAAAAATACATGTAg                                                                                 >  2:3361188/1‑98 (MQ=255)
                   acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgg                                                                                <  1:3275204/98‑1 (MQ=255)
                        tAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACATCTGGTAAAAATGTACATGTAGGGGCTg                                                                           <  1:1752915/98‑1 (MQ=255)
                          atatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgc                                                                         <  1:816092/98‑1 (MQ=255)
                            atTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                                       >  1:1134453/1‑98 (MQ=255)
                                gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                                   <  1:3583138/98‑1 (MQ=255)
                                    aGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTATCACTGAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAg                                                               >  2:1526189/1‑98 (MQ=255)
                                     gacGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCGCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGc                                                              <  2:1162294/96‑1 (MQ=255)
                                     gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGc                                                              <  1:1095216/98‑1 (MQ=255)
                                      tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                             <  1:2214724/98‑1 (MQ=255)
                                      tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                             >  1:2515309/1‑98 (MQ=255)
                                      tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                             >  1:3609572/1‑98 (MQ=255)
                                      tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                             >  1:1515478/1‑98 (MQ=255)
                                      tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTGGAAAAAATATACATGTAGGGGCTGGCGCAGCATTAGCa                                                             >  2:748986/1‑98 (MQ=255)
                                      tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTAGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                             >  2:1855362/1‑98 (MQ=255)
                                       ggTTGGCGATGGTAGAATGATTGGTAGGAATGCTGCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAg                                                            <  2:833525/97‑1 (MQ=255)
                                       cGTTGGCGAAGGTACAATGAGTGATATGAATGCTGCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAg                                                            <  2:262399/98‑1 (MQ=255)
                                         ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                          >  1:1324521/1‑98 (MQ=255)
                                          tGGCGGAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                                         <  2:1141885/98‑1 (MQ=255)
                                             cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                      <  1:1482324/98‑1 (MQ=255)
                                             cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                      <  1:703861/98‑1 (MQ=255)
                                              gAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                                     <  1:1224583/98‑1 (MQ=255)
                                              gAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTAGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                                     >  2:232726/1‑98 (MQ=255)
                                               aaGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTGAAAAAAAAATACAAGTAGGGGCAGGCGCAGTATTAGCAGGTGTGAtt                                                    >  2:1338358/1‑98 (MQ=255)
                                                aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                                   >  1:2475796/1‑98 (MQ=255)
                                                aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                                   >  1:1701047/1‑98 (MQ=255)
                                                aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGCTGCGATTg                                                   >  2:1639947/1‑98 (MQ=255)
                                                 ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                                  >  1:1065467/1‑98 (MQ=255)
                                                 ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                                  <  1:3554510/98‑1 (MQ=255)
                                                 ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                                  <  1:1922352/98‑1 (MQ=255)
                                                 ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                                  >  2:2616461/1‑98 (MQ=255)
                                                  gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                                 >  1:1503436/1‑98 (MQ=255)
                                                  gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                                 <  2:2745424/98‑1 (MQ=255)
                                                  gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                                 <  1:210640/98‑1 (MQ=255)
                                                  gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGGACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                                 <  2:648935/98‑1 (MQ=255)
                                                     cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccc                                              <  1:1804311/98‑1 (MQ=255)
                                                      aaTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                             <  1:3217304/98‑1 (MQ=255)
                                                      aaTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                             >  1:2591574/1‑98 (MQ=255)
                                                       aTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccccc                                            >  1:2978105/1‑98 (MQ=255)
                                                         gattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTa                                          <  1:1168531/98‑1 (MQ=255)
                                                           ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                                        >  1:3447060/1‑98 (MQ=255)
                                                                atGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATc                                   <  1:479680/98‑1 (MQ=255)
                                                                  gAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCAc                                 >  1:548388/1‑98 (MQ=255)
                                                                  gAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCAc                                 <  2:2093732/98‑1 (MQ=255)
                                                                   aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAGTGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCAcc                                <  2:2965959/98‑1 (MQ=255)
                                                                   aaTGCTACTCTCGGGGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCAct                                <  2:2474350/98‑2 (MQ=255)
                                                                   aaTGCGCCTCTCGGGGGTCGGGCGACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCAcc                                <  2:2720211/98‑1 (MQ=255)
                                                                    aTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCg                               <  1:2750592/98‑1 (MQ=255)
                                                                    aTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCg                               <  1:742528/98‑1 (MQ=255)
                                                                     tGCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCgg                              >  2:3140703/1‑98 (MQ=255)
                                                                      gCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAACCCCTAGTGCATCACCGGt                             >  2:1878807/1‑98 (MQ=255)
                                                                       ctcctctCGGTGGTCGTGCTACAGCTGGTAAAAATGTACATGTAGGGGCTGGTGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGtt                            <  2:3432401/95‑1 (MQ=255)
                                                                       cTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGtt                            <  2:1294952/98‑1 (MQ=255)
                                                                           tctcGGTGTTCGTGCAACAACTAGTAAAAATGTACATGTAGGGGCGGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATaa                        <  2:1065920/98‑1 (MQ=255)
                                                                           tctcGGTGGTCGTGGTACAAGTGGTAAGAGTGTACATGGAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATaa                        <  2:1600352/98‑1 (MQ=255)
                                                                           tctcGGTGGTCGTGCTACAGCTGGTAAAAATGTACATGGAGGGGCTGGCGAAGTATTAGCAGGTGGGATTGAACCCCCTAGTGCATCACCGGTTATaa                        <  2:2314514/98‑1 (MQ=255)
                                                                           tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATaa                        <  2:339165/98‑1 (MQ=255)
                                                                            ctcGGTGGTTGCGCTAAGACTGGGAAAAATGTACATGTAGGGGCTGGCGTAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTtatact                       <  2:3032853/98‑3 (MQ=255)
                                                                            ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAAt                       <  1:1027798/98‑1 (MQ=255)
                                                                                gtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAg                   >  1:932266/1‑98 (MQ=255)
                                                                                 tggtCGTGCTACCACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGCACCCCCTAGTGCATCACCGGTTATAATCGAgg                  >  2:2572704/1‑98 (MQ=255)
                                                                                 tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAgg                  >  2:3058805/1‑98 (MQ=255)
                                                                                 tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAgg                  >  2:1473766/1‑98 (MQ=255)
                                                                                 tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAgg                  <  1:2341414/98‑1 (MQ=255)
                                                                                 tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAgg                  >  1:3307963/1‑98 (MQ=255)
                                                                                  ggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGga                 >  2:3122456/1‑98 (MQ=255)
                                                                                   gtCGTGCTACAACTGGTGAAGGTGTACATGTAGGGGCGGGCGCAGTATTAGAAGGTGGGATTGAACCCCCTGGTGCATCACCGGTTATAATCGAGgat                <  2:3590131/98‑1 (MQ=255)
                                                                                    tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGgatg               >  1:395202/1‑98 (MQ=255)
                                                                                     cGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGgatga              <  1:1429213/98‑1 (MQ=255)
                                                                                      gTGCTACAACTGGTGAAAATGTACATGTAGGGGCTGGGGGAGTATTAGCAGGTGTGGTTGAACCCGCTAGCGCATCACCGGTTATAATCGAGgatgat             <  2:2770761/98‑1 (MQ=255)
                                                                                        gCTACAACTGGTAAGAAGGTGGATGTAGGGGCTGGCGGAGCATAAGCAGTACTAATTGAGCCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGt           <  2:2944104/98‑1 (MQ=255)
                                                                                         cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGACCCCCCTAGTGCATCACCGGTTATAATCGAGGatgatgca          <  2:3528382/98‑3 (MQ=255)
                                                                                         cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTa          >  2:211857/1‑98 (MQ=255)
                                                                                         cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTa          <  2:2204415/98‑1 (MQ=255)
                                                                                         cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTa          <  2:1421408/98‑1 (MQ=255)
                                                                                             aaCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTaa      <  1:2250496/98‑1 (MQ=255)
                                                                                               cTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATc    >  1:1267900/1‑98 (MQ=255)
                                                                                               cTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATc    <  1:2837916/98‑1 (MQ=255)
                                                                                                tGGTAAAAATGTACATGTAGGGGCTGGGGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATCg   <  1:1367782/98‑1 (MQ=255)
                                                                                                 ggTAAAGATGTTCATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATCgg  <  1:1829843/98‑1 (MQ=255)
                                                                                                 ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCGGGTGTGATTGAACCCCCTGGTGCGTCACCGGTTATAATCGAGGATGATGTATTAATggg  <  2:1486395/98‑1 (MQ=255)
                                                                                                 ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATCgg  >  2:3508877/1‑98 (MQ=255)
                                                                                                 ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATCgg  <  1:3364709/98‑1 (MQ=255)
                                                                                                 |                                                                                                 
TGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGTATTAATCGG  >  BX571857/1434268‑1434462

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