Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 2,602,946 G→T 77.3% G116G (GGC→GGA SAS2426 ← putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718572,602,9460GT77.3% 482.2 / 227.7 229G116G (GGC→GGASAS2426putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/29);  new base T (97/80);  total (120/109)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-02

Only 100 of 233 total aligned reads displayed.

ATCGTATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCTTTAATTTTC  >  BX571857/2602851‑2603042
                                                                                               |                                                                                                
  cGTATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCa                                                                                               >  1:1337821/1‑97 (MQ=255)
   gTATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCaa                                                                                              >  1:121928/1‑97 (MQ=255)
   gTATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCaa                                                                                              >  2:2645512/1‑97 (MQ=255)
   gTATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCaa                                                                                              >  1:473622/1‑97 (MQ=255)
    tATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCaaa                                                                                             <  1:878906/97‑1 (MQ=255)
    tATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCaaa                                                                                             >  1:1848340/1‑97 (MQ=255)
     aTCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAAt                                                                                            >  1:3217640/1‑97 (MQ=255)
      tCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATg                                                                                           <  1:148546/97‑1 (MQ=255)
       cTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGc                                                                                          <  1:502451/97‑1 (MQ=255)
       cTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGc                                                                                          <  1:1900284/97‑1 (MQ=255)
        tGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCt                                                                                         >  2:1226350/1‑97 (MQ=255)
        tGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCt                                                                                         >  2:7224/1‑97 (MQ=255)
           accaccTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTa                                                                                      <  2:2732181/97‑1 (MQ=255)
            ccaccTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTaa                                                                                     <  2:629858/97‑1 (MQ=255)
            ccaccTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTaa                                                                                     >  1:2573933/1‑97 (MQ=255)
            ccaccTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTaa                                                                                     >  2:1884803/1‑97 (MQ=255)
             caccTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAg                                                                                    >  1:2332379/1‑97 (MQ=255)
                   gTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTaa                                                                              >  1:2815491/1‑97 (MQ=255)
                    ttCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAg                                                                             <  1:1376166/97‑1 (MQ=255)
                    ttCATGTGACACGTCGATACCTACGGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAg                                                                             <  2:3316843/97‑1 (MQ=255)
                     tCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGa                                                                            >  1:2498013/1‑97 (MQ=255)
                      cATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGaa                                                                           >  2:739984/1‑97 (MQ=255)
                      cATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAACATCCAAATGCTTTAAGTGTTAAGaa                                                                           >  2:254956/1‑97 (MQ=255)
                        tgtgACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAaca                                                                         <  1:3308580/97‑1 (MQ=255)
                         gtgACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAacac                                                                        >  2:36406/1‑97 (MQ=255)
                         gtgACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAacac                                                                        >  2:1716773/1‑97 (MQ=255)
                         gtgACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAacac                                                                        >  1:204388/1‑97 (MQ=255)
                          tgACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACt                                                                       >  1:948385/1‑97 (MQ=255)
                          tgACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACt                                                                       <  1:1131005/97‑1 (MQ=255)
                            acacTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTaa                                                                     >  2:1275445/1‑97 (MQ=255)
                             cacTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATGACGTTACAATGACGTCTTTAGAATATACAAATGCTTTACGTGGTAAGAACACTAAc                                                                    >  2:447469/1‑97 (MQ=255)
                             cacTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAAc                                                                    <  1:2386025/97‑1 (MQ=255)
                             cacTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAAc                                                                    <  2:3052546/97‑1 (MQ=255)
                               cTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACAt                                                                  >  1:2286728/1‑97 (MQ=255)
                                 tCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGa                                                                >  2:644116/1‑97 (MQ=255)
                                 tCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGa                                                                <  1:3246651/97‑1 (MQ=255)
                                  cGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGaa                                                               >  2:2330811/1‑97 (MQ=255)
                                       ccTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCa                                                          >  1:2158805/1‑97 (MQ=255)
                                        cTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAg                                                         <  2:1756413/97‑1 (MQ=255)
                                        cTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAg                                                         <  1:2340780/97‑1 (MQ=255)
                                         tACTGCGGCTGTTATAGATCGGGGCACTAACGCTAAAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAAAAATAACATGAAATTCAGc                                                        >  2:2266828/1‑97 (MQ=255)
                                         tACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGc                                                        <  1:1703386/97‑1 (MQ=255)
                                         tACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGc                                                        >  1:287903/1‑97 (MQ=255)
                                                gCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGc                                                 <  1:1739248/97‑1 (MQ=255)
                                                 cTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGcc                                                <  2:3323699/97‑1 (MQ=255)
                                                 cTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGcc                                                <  2:1867443/97‑1 (MQ=255)
                                                  tGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCa                                               <  2:2754254/97‑1 (MQ=255)
                                                  tGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCa                                               >  1:1400683/1‑97 (MQ=255)
                                                  tGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCa                                               >  1:2994974/1‑97 (MQ=255)
                                                    ttATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGt                                             >  2:898640/1‑97 (MQ=255)
                                                     tataGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGtt                                            >  1:1527819/1‑97 (MQ=255)
                                                     tataGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGtt                                            >  2:1403748/1‑97 (MQ=255)
                                                      ataGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTa                                           >  2:2438506/1‑97 (MQ=255)
                                                      ataGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTa                                           <  1:1782713/97‑1 (MQ=255)
                                                      ataGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTa                                           <  1:408374/97‑1 (MQ=255)
                                                      ataGATCGCGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTa                                           <  2:220065/97‑1 (MQ=255)
                                                        aGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAAt                                         >  2:328745/1‑97 (MQ=255)
                                                         gATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGataata                                        >  1:692381/1‑97 (MQ=255)
                                                         gATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATa                                        >  2:1564384/1‑97 (MQ=255)
                                                         gATCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATa                                        >  1:2774073/1‑97 (MQ=255)
                                                          aTCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATac                                       <  2:2060351/97‑1 (MQ=255)
                                                           tCGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAGATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATaca                                      <  2:148538/97‑1 (MQ=255)
                                                            cGGGGCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATacac                                     <  1:1111453/97‑1 (MQ=255)
                                                               ggCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTg                                  >  2:1181093/1‑97 (MQ=255)
                                                               ggCAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTg                                  <  1:1990128/97‑1 (MQ=255)
                                                                 cAATAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCa                                >  2:2603444/1‑97 (MQ=255)
                                                                   aTAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCaaa                              <  2:2452586/97‑1 (MQ=255)
                                                                   aTAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCaaa                              <  1:483559/97‑1 (MQ=255)
                                                                   aTAACGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCaaa                              <  1:3052713/97‑1 (MQ=255)
                                                                   aTAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCaaa                              <  2:1591319/97‑1 (MQ=255)
                                                                      aCGTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAAt                           <  2:944447/97‑1 (MQ=255)
                                                                       cGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATg                          >  1:1724611/1‑97 (MQ=255)
                                                                        gTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGa                         <  2:695018/97‑1 (MQ=255)
                                                                        gTTACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGa                         >  1:241001/1‑97 (MQ=255)
                                                                         ttACAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGAt                        <  1:371380/97‑1 (MQ=255)
                                                                          taacaTGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGCTa                       <  2:1503777/93‑1 (MQ=255)
                                                                          tACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATa                       <  1:919802/97‑1 (MQ=255)
                                                                           aCAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATag                      <  1:636496/97‑1 (MQ=255)
                                                                            cAATGACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATaga                     <  1:2146376/97‑1 (MQ=255)
                                                                               ggACGTCGTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGAc                  <  2:600602/96‑1 (MQ=255)
                                                                                gACGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACa                 <  2:2942509/97‑1 (MQ=255)
                                                                                 aCGTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACAt                >  2:1576158/1‑97 (MQ=255)
                                                                                  cGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACAtt               <  2:1543328/97‑1 (MQ=255)
                                                                                   gTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTg              <  1:2691733/97‑1 (MQ=255)
                                                                                   gTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTg              >  2:1776531/1‑97 (MQ=255)
                                                                                   gTCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTg              >  1:224212/1‑97 (MQ=255)
                                                                                   gTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTg              <  2:2816332/97‑1 (MQ=255)
                                                                                   gTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTg              <  1:533397/97‑1 (MQ=255)
                                                                                   gTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTg              <  1:2129897/97‑1 (MQ=255)
                                                                                    tCTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGt             <  2:113788/97‑1 (MQ=255)
                                                                                     cTTTAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTc            >  1:444688/1‑97 (MQ=255)
                                                                                       ttAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCtt          >  2:2934149/1‑97 (MQ=255)
                                                                                        tAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCCGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCttt         >  2:508119/1‑97 (MQ=255)
                                                                                        tAGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCttt         >  1:1800088/1‑97 (MQ=255)
                                                                                        tAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCttt         >  2:1373564/1‑97 (MQ=255)
                                                                                         aGAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCTTTa        >  1:1473857/1‑97 (MQ=255)
                                                                                          gAATATCCAAATGCTTTAAGTGTTAAGAACACTAACGTGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCTTTaa       <  1:1034958/97‑1 (MQ=255)
                                                                                          gAATATCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCTTTaa       <  2:2385471/97‑1 (MQ=255)
                                                                                            atatCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCTTTAAtt     <  2:800690/97‑1 (MQ=255)
                                                                                             tAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCTTTAAttt    >  1:2622647/1‑97 (MQ=255)
                                                                                               |                                                                                                
ATCGTATCTGTACCACCTAGTTCATGTGACACTTCGATACCTACTGCGGCTGTTATAGATCGGGGCAATAACGTTACAATGACGTCTTTAGAATAGCCAAATGCTTTAAGTGTTAAGAACACTAACATGAAATTCAGCATGACGCCAGTTAATACACTTGCAAAAATGATAGAGACATTGTCTTTAATTTTC  >  BX571857/2602851‑2603042

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