Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 2,407,915 G→A 41.4% E35K (GAA→AAA)  SAS2252 → response regulator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718572,407,9150GA41.4% 145.3 / inf 222E35K (GAA→AAA) SAS2252response regulator protein
Reads supporting (aligned to +/- strand):  ref base G (56/74);  new base A (47/45);  total (103/119)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

CAATGTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGTAATACAT  >  BX571857/2407819‑2408011
                                                                                                |                                                                                                
cAATGTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAgtggtg                                                                                                  <  1:1400131/97‑1 (MQ=255)
 aaTGTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGa                                                                                                 <  2:2319457/97‑1 (MQ=255)
  aTGTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGaa                                                                                                >  1:1219484/1‑97 (MQ=255)
  aTGTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGaa                                                                                                >  1:1982114/1‑97 (MQ=255)
   tGTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAg                                                                                               <  1:795609/97‑1 (MQ=255)
    gTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGc                                                                                              >  2:2378941/1‑97 (MQ=255)
    gTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGc                                                                                              >  1:340681/1‑97 (MQ=255)
    gTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGc                                                                                              >  1:679238/1‑97 (MQ=255)
    gTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGc                                                                                              <  1:2401779/97‑1 (MQ=255)
     tCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACACTGATGCATACACACAACCAAGTGGTGAAGCt                                                                                             >  2:182238/1‑97 (MQ=255)
      cTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTg                                                                                            <  2:796682/97‑1 (MQ=255)
      cTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTg                                                                                            <  1:3667962/97‑1 (MQ=255)
       ttgttgTCGATGACGATCCTCGAGTCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGc                                                                                           <  2:3560770/97‑1 (MQ=255)
        tgttgtCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACATACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCa                                                                                          <  2:2835319/97‑1 (MQ=255)
        tgttgtCGATGACGATCCTCGAATCCTAAATTATAAAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCa                                                                                          <  2:1419291/97‑1 (MQ=255)
            gtCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTaa                                                                                      >  1:1021527/1‑97 (MQ=255)
              cGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTaaaa                                                                                    >  2:245562/1‑97 (MQ=255)
                aTGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACt                                                                                  <  2:3069428/97‑1 (MQ=255)
                  gACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACttt                                                                                <  2:292776/97‑1 (MQ=255)
                  gACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACttt                                                                                <  2:3510248/97‑1 (MQ=255)
                     gATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAg                                                                             >  2:3317561/1‑97 (MQ=255)
                      aTCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGa                                                                            >  1:360725/1‑97 (MQ=255)
                        ccTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGaaa                                                                          >  1:2310856/1‑97 (MQ=255)
                         cTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGaaaa                                                                         >  1:2565438/1‑97 (MQ=255)
                         cTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGaaaa                                                                         <  1:800003/97‑1 (MQ=255)
                           cGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAAc                                                                       >  1:2514357/1‑97 (MQ=255)
                               tCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCg                                                                   >  1:3580504/1‑97 (MQ=255)
                                ccTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCgt                                                                  <  2:1682808/97‑1 (MQ=255)
                                ccTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCgt                                                                  <  1:1169488/97‑1 (MQ=255)
                                ccTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCgc                                                                  >  2:189595/1‑96 (MQ=255)
                                 cTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCgtg                                                                 <  1:667540/97‑1 (MQ=255)
                                 cTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCgtg                                                                 >  1:266836/1‑97 (MQ=255)
                                     aTTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGa                                                             <  2:2062643/97‑1 (MQ=255)
                                       tatataGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGata                                                           >  2:92305/1‑97 (MQ=255)
                                          ataGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTg                                                        <  2:640356/97‑1 (MQ=255)
                                          ataGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTg                                                        <  1:3087713/97‑1 (MQ=255)
                                           tagctagcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGc                                                       >  2:199642/1‑97 (MQ=255)
                                           tagctagcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGc                                                       <  2:3208748/97‑1 (MQ=255)
                                              ctcgcCATTTACAAACAGAGCACATTGATGCGTACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGt                                                    <  2:2408532/94‑1 (MQ=255)
                                              ctagcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGt                                                    >  2:2283542/1‑97 (MQ=255)
                                              ctagcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGt                                                    >  2:164491/1‑97 (MQ=255)
                                               tagcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTg                                                   >  1:1663979/1‑97 (MQ=255)
                                               tagcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTg                                                   >  1:1329868/1‑97 (MQ=255)
                                               tagcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTg                                                   >  1:602750/1‑97 (MQ=255)
                                                agcCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTgg                                                  <  2:610448/97‑1 (MQ=255)
                                                   cATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAg                                               <  1:2750032/97‑1 (MQ=255)
                                                   cATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAg                                               >  2:1150289/1‑97 (MQ=255)
                                                    aTTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGGGGTAGa                                              >  2:2756103/1‑97 (MQ=255)
                                                     tttACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGat                                             <  2:3478812/97‑1 (MQ=255)
                                                     tttACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGat                                             <  1:2976928/97‑1 (MQ=255)
                                                       tACAAACAGAGCACATTGATGCATCCACACAACCAAGTGGTGAAGCTGCATTAAAACTCTTAGGAAAACAGCGTGTCGATATTGCAGTGGTAGatat                                           >  2:822358/1‑97 (MQ=255)
                                                       tACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGatat                                           >  1:927177/1‑97 (MQ=255)
                                                       tACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGatat                                           >  2:122/1‑97 (MQ=255)
                                                        aCAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGAtatt                                          <  2:299654/97‑1 (MQ=255)
                                                          aaaCAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGAtattat                                        <  1:1855797/97‑1 (MQ=255)
                                                           aaCAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTatg                                       >  1:2752871/1‑97 (MQ=255)
                                                           aaCAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGCAGAtattacg                                       >  2:2224647/1‑95 (MQ=255)
                                                             cAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTatgaa                                     >  1:311564/1‑96 (MQ=255)
                                                             cAGAGCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTatgat                                     >  1:3107800/1‑97 (MQ=255)
                                                              agagCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTatgatg                                    <  2:881962/97‑1 (MQ=255)
                                                              agagCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTatgatg                                    <  1:2930264/97‑1 (MQ=255)
                                                               gagCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATgatgg                                   <  1:3192006/97‑1 (MQ=255)
                                                                agCACATTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATgatgga                                  <  1:1598504/97‑1 (MQ=255)
                                                                  cacaTTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATgatggatg                                <  1:1731270/97‑1 (MQ=255)
                                                                  cacaTTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATgatggatg                                <  2:301447/97‑1 (MQ=255)
                                                                   acaTTGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGatggatgg                               <  1:501546/97‑1 (MQ=255)
                                                                       tGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATg                           >  2:3459524/1‑97 (MQ=255)
                                                                       tGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTGTAGAAAAACAGCGTGTCGACATTGCTGTGGTAGATATTATGATGGATGGTATg                           >  2:1359037/1‑97 (MQ=255)
                                                                       tGATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATg                           <  2:502846/97‑1 (MQ=255)
                                                                        gATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATgg                          >  1:3565336/1‑97 (MQ=255)
                                                                        gATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATgg                          >  2:2529309/1‑97 (MQ=255)
                                                                        gATGCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATgg                          <  2:756946/97‑1 (MQ=255)
                                                                           gCATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACg                       <  1:272407/97‑1 (MQ=255)
                                                                            cATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACgg                      >  2:700686/1‑97 (MQ=255)
                                                                            cATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACgg                      >  1:1094581/1‑97 (MQ=255)
                                                                            cATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACgg                      >  1:22592/1‑97 (MQ=255)
                                                                            cATACACACAACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACgg                      <  1:462365/97‑1 (MQ=255)
                                                                                cacacaACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCttt                  >  1:2077127/1‑97 (MQ=255)
                                                                                 acacaACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTc                 <  2:2290370/97‑1 (MQ=255)
                                                                                 acacaACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTc                 >  1:327812/1‑97 (MQ=255)
                                                                                 acacaACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTc                 <  1:2821547/97‑1 (MQ=255)
                                                                                 acacaACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTc                 >  1:1072596/1‑97 (MQ=255)
                                                                                 acacaACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTc                 <  1:2540416/97‑1 (MQ=255)
                                                                                  cacaACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCa                <  1:216013/97‑1 (MQ=255)
                                                                                   acaACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCaa               <  1:2667152/97‑1 (MQ=255)
                                                                                   acaACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCaa               <  1:1124315/97‑1 (MQ=255)
                                                                                   acaACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCaa               <  2:464480/97‑1 (MQ=255)
                                                                                   acaACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCaa               <  1:2598308/97‑1 (MQ=255)
                                                                                    cgaCCAAGTGGTGGGGCTGCATTAAAACTTTTAGGAGAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAAt              <  2:2145586/95‑1 (MQ=255)
                                                                                    caACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAAt              >  2:1465395/1‑97 (MQ=255)
                                                                                    caACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAAt              <  2:555861/97‑1 (MQ=255)
                                                                                    caACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAAt              <  1:223186/97‑1 (MQ=255)
                                                                                    caACCAAGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAAt              <  1:3464260/97‑1 (MQ=255)
                                                                                         aaGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGt         <  1:1937584/97‑1 (MQ=255)
                                                                                         aaGTGGTAAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGt         <  2:3616451/97‑1 (MQ=255)
                                                                                           gtggtgAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGTaa       >  1:2888348/1‑97 (MQ=255)
                                                                                           gtggtaAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGCCGGCTTTCAATTATGTaa       >  2:1099948/1‑97 (MQ=255)
                                                                                           gtggtaAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGTaa       <  2:1555856/97‑1 (MQ=255)
                                                                                             ggtgAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGTAATa     <  1:1517714/97‑1 (MQ=255)
                                                                                               tgAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGTAATACa   <  2:3377925/97‑1 (MQ=255)
                                                                                                |                                                                                                
CAATGTCTTGTTGTCGATGACGATCCTCGAATCCTAAATTATATAGCTAGCCATTTACAAACAGAGCACATTGATGCATACACACAACCAAGTGGTGAAGCTGCATTAAAACTTTTAGAAAAACAGCGTGTCGATATTGCAGTGGTAGATATTATGATGGATGGTATGGACGGCTTTCAATTATGTAATACAT  >  BX571857/2407819‑2408011

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