Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 495,714 T→C 6.1% intergenic (+552/‑352) unknown → / → SAS0438 –/Orn/Lys/Arg decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857495,7140TC 6.1% 762.9 / 19.9 205intergenic (+552/‑352)unknown/SAS0438–/Orn/Lys/Arg decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (73/118);  new base C (8/6);  total (81/124)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

CGAAAAAGAGGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACTAGAAGTATA  >  BX571857/495618‑495810
                                                                                                 |                                                                                                
cGAAAAAGAGGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCt                                                                                                  >  1:2884477/1‑97 (MQ=255)
cGAAAAAGAGGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGCTCCGCTAGAGTAGAACGTTGCt                                                                                                  >  2:1140562/1‑97 (MQ=255)
    ggagagGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGc                                                                                              <  2:1652523/95‑1 (MQ=255)
    aaaGAGGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGc                                                                                              >  1:882741/1‑97 (MQ=255)
    aaaGAGGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGc                                                                                              >  1:86499/1‑97 (MQ=255)
      agagGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCaa                                                                                            >  1:1589919/1‑97 (MQ=255)
         ggAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGAAGAACGTTGCTAGGCAAGaa                                                                                         >  2:1035935/1‑97 (MQ=255)
          gAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGaaa                                                                                        <  2:912896/97‑1 (MQ=255)
           aaaGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAAt                                                                                       >  1:2513326/1‑97 (MQ=255)
           aaaGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAAt                                                                                       >  1:2088946/1‑97 (MQ=255)
           aaaGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAAt                                                                                       <  1:1954202/97‑1 (MQ=255)
            aagaagCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATg                                                                                      <  1:684538/97‑1 (MQ=255)
             agaagCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATgg                                                                                     >  1:554022/1‑97 (MQ=255)
              gaagCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGa                                                                                    <  1:1349728/97‑1 (MQ=255)
               aagcaagcAATTGTCACAATTCAAGAAAGGCCTTTAGCGACGATGGTAGCCAACTTCCGTTCCGCTCGAGCAGACCGTTGCTCGGCTAGACATgcac                                                                                   >  2:1446464/1‑94 (MQ=14)
               aagcaagcAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGAt                                                                                   >  1:1025153/1‑97 (MQ=255)
                 gcaagcaaTTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGc                                                                                 <  1:675027/97‑1 (MQ=255)
                    agcaaTTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGAt                                                                              >  1:44998/1‑97 (MQ=255)
                    agcaaTTGTCACAAGTCAAGAAAGGTCCTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGAt                                                                              <  2:2801829/97‑1 (MQ=255)
                      cgaTTGTCGCAAGTCAAGAAAGGTCTGTAGCGACGATGGTAGCCTACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGTAATGGATGCGATGa                                                                            <  2:1491926/95‑1 (MQ=14)
                      caaTTGTCACAAGTCAAGAAAGGTCTTTAGCGCCGATGGTAGCCAACTTACGTTCCGCTAAAGTAGAACGTTGCTAGGCAAGAAAGGGATGCGATGa                                                                            >  1:1264919/1‑97 (MQ=255)
                      caaTTGTCACAAGTCAAGAAAGGTCTTTAGCGCCGACGGTAGGCAACTTACGTTTCGCTAGAGCAGCGCGTTGCTAGTCAAGAAATGGATGCGATGa                                                                            <  2:2340886/97‑1 (MQ=17)
                      caaTTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGa                                                                            >  1:1929436/1‑97 (MQ=255)
                       aaTTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAg                                                                           >  2:2483747/1‑97 (MQ=255)
                       aaTTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAg                                                                           >  1:506729/1‑97 (MQ=255)
                          tGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCg                                                                        <  1:1753711/97‑1 (MQ=255)
                          tGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCg                                                                        <  1:2029509/97‑1 (MQ=255)
                           gTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGc                                                                       >  2:1691333/1‑97 (MQ=255)
                            tCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCa                                                                      <  2:765329/97‑1 (MQ=255)
                              acaAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCAtt                                                                    >  2:2087682/1‑97 (MQ=255)
                               caAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTg                                                                   <  2:1370740/97‑1 (MQ=255)
                               caAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTg                                                                   <  1:254320/97‑1 (MQ=255)
                                aaGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGa                                                                  <  1:1795777/97‑1 (MQ=255)
                                 aGTCAAGAAAGATCTTTAGCGCCGTTGTTAGCCAACTTACGTTCCTCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAg                                                                 <  1:1651839/97‑1 (MQ=17)
                                  gTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGt                                                                <  2:1364920/97‑1 (MQ=255)
                                    cAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGttt                                                              <  2:2969434/97‑1 (MQ=255)
                                    cAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGttt                                                              >  2:725945/1‑97 (MQ=255)
                                       gAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGCTTggg                                                           >  2:1389507/1‑97 (MQ=255)
                                         aaGGACTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGCAGAACGTCGCTAGGCGAGAAATGGATGCGCTCAGCCGCATTGAGTATTGGaa                                                         >  2:1928644/1‑97 (MQ=17)
                                           ggTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAAt                                                       >  1:779887/1‑97 (MQ=255)
                                           ggTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAAt                                                       <  1:2737270/97‑1 (MQ=255)
                                           ggTCCTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAAt                                                       <  2:637448/97‑1 (MQ=255)
                                            gTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATg                                                      >  2:331849/1‑97 (MQ=255)
                                            gTCCTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATg                                                      <  2:1888692/97‑1 (MQ=255)
                                             tCTTTAGCGACGATGGGAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTCGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGa                                                     <  2:1335962/97‑1 (MQ=255)
                                              ttttAGCGACGATGGTAGCCAACTTACGGTCCGCTAGAGGAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGat                                                    <  2:1761404/96‑1 (MQ=255)
                                              cTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGat                                                    <  1:2972206/97‑1 (MQ=255)
                                              cTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGat                                                    <  1:2805572/97‑1 (MQ=255)
                                               tttAGCGTCGAGGGTAGCAAACTAACGTTCCGCTAGAGAAGGACGTTGCTAGGCAAGAAAGGGTTACGCAGAGCCGCATGGAGGTTTGGAAGGGata                                                   >  2:1582658/1‑97 (MQ=9)
                                               tttAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGata                                                   <  2:2742312/97‑1 (MQ=255)
                                                 tAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGatat                                                  <  2:598317/97‑1 (MQ=11)
                                                 tAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATAtt                                                 >  1:424767/1‑97 (MQ=255)
                                                   gCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATAtttt                                               <  2:987665/97‑1 (MQ=255)
                                                    cgacgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTta                                              <  1:117814/97‑1 (MQ=255)
                                                    cgacgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTta                                              <  1:2263584/97‑1 (MQ=255)
                                                     gacgaTGGTATCCAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAGAAAGGGATCCGAGGGCCCCCATGGGGCTTTGGAATTGATATTTtaa                                             >  2:657001/1‑97 (MQ=9)
                                                     gacgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaa                                             >  1:1099825/1‑97 (MQ=255)
                                                     gacgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaa                                             <  2:1281589/97‑1 (MQ=255)
                                                      gcgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaat                                            <  1:2117334/96‑1 (MQ=255)
                                                      acgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaat                                            <  2:782941/97‑1 (MQ=255)
                                                       cgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaata                                           <  2:1133386/97‑1 (MQ=255)
                                                       cgaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaata                                           <  1:2181897/97‑1 (MQ=255)
                                                       agaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaata                                           <  1:562865/96‑1 (MQ=255)
                                                        gaTGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTtaataa                                          >  1:701852/1‑97 (MQ=255)
                                                                 ccACCTTCCGTTCGGCTAAGGTAGAACGTTGTCAGGAAAAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACtt                                 <  2:1254619/97‑1 (MQ=18)
                                                                  cAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGCAATGATATTTTAATAATGTGCACtt                                 >  1:1303194/1‑97 (MQ=255)
                                                                   aaCTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGGTGAGGAAATGATATTTTAATAATGTGCACttt                                >  2:2685076/1‑97 (MQ=255)
                                                                    aCTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTg                                                                   >  2:1406884/1‑61 (MQ=255)
                                                                    aCTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTg                                                                   <  1:1406885/61‑1 (MQ=255)
                                                                    aCTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACtttt                               >  2:173030/1‑97 (MQ=255)
                                                                    aCTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACtttt                               <  2:28315/97‑1 (MQ=255)
                                                                     cTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCAtt                                                                    <  1:341508/59‑1 (MQ=255)
                                                                     cTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTg                              <  2:1952450/97‑1 (MQ=255)
                                                                     cTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTg                              <  2:2619502/97‑1 (MQ=255)
                                                                     cTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTg                              <  1:1431629/97‑1 (MQ=255)
                                                                      ttACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGTTCGGGAAGTGATATTTTTTTAATGTGCACATTTGa                             >  1:1135321/1‑97 (MQ=14)
                                                                       tACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGAt                            <  1:2303210/97‑1 (MQ=255)
                                                                           ttCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGtt                        >  2:895922/1‑97 (MQ=255)
                                                                            tCCGCTAGAGTAGAACGTTGCCAGGCGAGAAATGGATGCGATGGGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGttt                       <  2:470841/97‑1 (MQ=255)
                                                                             ccGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCAtt                                                                    <  2:2718644/51‑1 (MQ=255)
                                                                             ccGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTa                      <  1:1825535/97‑1 (MQ=255)
                                                                             ccGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTa                      >  1:496372/1‑97 (MQ=255)
                                                                             ccGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTa                      >  2:434963/1‑97 (MQ=255)
                                                                                 tAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATc                  >  2:1227585/1‑97 (MQ=255)
                                                                                 tAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATc                  <  2:366251/97‑1 (MQ=255)
                                                                                  agagTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCt                 <  1:2788249/97‑1 (MQ=255)
                                                                                   gagTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCtg                <  1:1211393/97‑1 (MQ=255)
                                                                                    agTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCtgt               <  1:1150682/97‑1 (MQ=255)
                                                                                       aGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGAc            >  2:1960162/1‑97 (MQ=255)
                                                                                        gAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACt           <  1:1238960/97‑1 (MQ=255)
                                                                                        gAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACt           <  2:3008840/97‑1 (MQ=255)
                                                                                         aaCGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACt                                  <  1:2437221/73‑1 (MQ=255)
                                                                                         aaCGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACt                                  >  2:2437220/1‑73 (MQ=21)
                                                                                          aCGTTGCCAGGCAAAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAAt                                                       <  1:1243894/50‑1 (MQ=11)
                                                                                           cGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTg                          >  2:1822643/1‑79 (MQ=25)
                                                                                            gTTGGTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACTAGaa       <  2:859268/97‑1 (MQ=255)
                                                                                             ttGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACTAGAAg      <  1:892582/97‑1 (MQ=255)
                                                                                             ttGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACTAGAAg      <  2:201362/97‑1 (MQ=255)
                                                                                              tGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACTAGAAGt     >  1:2295196/1‑97 (MQ=255)
                                                                                                 tAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACTAGAAGtata  <  1:610314/97‑1 (MQ=255)
                                                                                                 |                                                                                                
CGAAAAAGAGGAAAGAAGCAAGCAATTGTCACAAGTCAAGAAAGGTCTTTAGCGACGATGGTAGCCAACTTACGTTCCGCTAGAGTAGAACGTTGCTAGGCAAGAAATGGATGCGATGAGCCGCATTGAGTTTTGGAAATGATATTTTAATAATGTGCACTTTTGATTGTTTAAATCTGTGACTAGAAGTATA  >  BX571857/495618‑495810

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: