Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 549,806 A→G 70.5% intergenic (‑141/‑17) SAS0478 ← / → SAS0479 nucleoside permease/putative DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857549,8060AG70.5% 424.5 / 267.9 247intergenic (‑141/‑17)SAS0478/SAS0479nucleoside permease/putative DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (48/25);  new base G (103/71);  total (151/96)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 254 total aligned reads displayed.

AAAATGGAAGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTGAAATT  >  BX571857/549710‑549900
                                                                                                |                                                                                              
aaaaTGGAAGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCa                                                                                                >  2:2513669/1‑97 (MQ=255)
   aTGGAAGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGACAAAATAGGTCATATAATATAGTCAAAGAAGGTCgaaa                                                                                             >  1:1659934/1‑97 (MQ=255)
    tGGAAGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCgaaaa                                                                                            >  1:3097373/1‑97 (MQ=255)
    tGGAAGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCaaaaa                                                                                            >  1:2890915/1‑97 (MQ=255)
      gAAGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAgg                                                                                          <  1:1561076/97‑1 (MQ=255)
        aGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAgggg                                                                                        >  1:1309067/1‑97 (MQ=255)
        aGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAgggg                                                                                        >  1:193184/1‑97 (MQ=255)
          ggCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTg                                                                                      >  1:1182519/1‑97 (MQ=255)
          ggCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTg                                                                                      <  1:2247382/97‑1 (MQ=255)
           gCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGgaga                                                                                     >  1:27124/1‑97 (MQ=255)
            cTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGat                                                                                    >  2:393468/1‑97 (MQ=255)
               aCATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGAtata                                                                                 >  2:407399/1‑97 (MQ=255)
                 atcatAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACa                                                                               <  1:1687328/97‑1 (MQ=255)
                 atcatAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACa                                                                               >  2:3028311/1‑97 (MQ=255)
                 atcatAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACa                                                                               >  2:1351317/1‑97 (MQ=255)
                   catAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATg                                                                             >  1:2117167/1‑97 (MQ=255)
                   catAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATg                                                                             >  1:1656775/1‑97 (MQ=255)
                    atAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGc                                                                            >  2:2422288/1‑97 (MQ=255)
                     tAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGca                                                                           <  1:3193716/97‑1 (MQ=255)
                      aaaTTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGcac                                                                          <  2:2268008/97‑1 (MQ=255)
                      aaaTTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGcac                                                                          <  2:566209/97‑1 (MQ=255)
                      aaaTTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGcac                                                                          <  2:1516429/97‑1 (MQ=255)
                       aaTTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGcaca                                                                         >  1:915473/1‑97 (MQ=255)
                       aaTTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGcaca                                                                         >  2:3343341/1‑97 (MQ=255)
                        aTTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACaa                                                                        <  1:1252674/97‑1 (MQ=255)
                          ttttGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAata                                                                      >  1:2919910/1‑97 (MQ=255)
                           tttGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAatat                                                                     >  1:2775248/1‑97 (MQ=255)
                            ttGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATg                                                                    <  2:1194939/97‑1 (MQ=255)
                            ttGAAAATATAAAAACAAATTAATTGAAAAAGGGCAAAATAGGTCATATAATATAGTCAAAGCAGGTCAAAAAGGGGTGATATACAGGCACAATATg                                                                    >  2:2884372/1‑97 (MQ=255)
                             tGAAAATATAAAACCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATCGTCAAAGAAGGTCGAAAAGGGGAGATATACATGCACAATATGt                                                                   <  1:834295/97‑1 (MQ=255)
                             tGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGt                                                                   >  2:422388/1‑97 (MQ=255)
                               aaaaTATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCt                                                                 <  2:2632691/97‑1 (MQ=255)
                               aaaaTATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCt                                                                 >  1:1865674/1‑97 (MQ=255)
                                   tataAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGAca                                                             >  2:2865022/1‑97 (MQ=255)
                                      aaaaaCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGAcatca                                                          >  2:2494821/1‑97 (MQ=255)
                                       aaaaCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGAcatcat                                                         >  2:3185336/1‑97 (MQ=255)
                                       aaaaCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGAcatcat                                                         >  1:1772380/1‑97 (MQ=255)
                                         aaCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAg                                                       >  2:702107/1‑97 (MQ=255)
                                         aaCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAg                                                       <  1:2296940/97‑1 (MQ=255)
                                         aaCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAg                                                       >  2:1883194/1‑97 (MQ=255)
                                          aCAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGa                                                      >  2:2093538/1‑97 (MQ=255)
                                            aaaTTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAAc                                                    >  2:2329362/1‑97 (MQ=255)
                                            aaaTTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAAc                                                    <  2:1566890/97‑1 (MQ=255)
                                            aaaTTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAAc                                                    >  2:2694292/1‑97 (MQ=255)
                                             aattaattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCCCAATATGTCTGACATCATAGAACa                                                   >  2:1925387/1‑97 (MQ=255)
                                             aattaattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACa                                                   >  2:597160/1‑97 (MQ=255)
                                               ttaattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAAt                                                 >  1:779352/1‑97 (MQ=255)
                                               ttaattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAAt                                                 <  1:1721925/97‑1 (MQ=255)
                                               ttaattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAAt                                                 >  2:2259196/1‑97 (MQ=255)
                                               ttaattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAAt                                                 >  1:1955290/1‑97 (MQ=255)
                                                taattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATa                                                <  1:1171815/97‑1 (MQ=255)
                                                taattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATa                                                <  1:2648374/97‑1 (MQ=255)
                                                 aattGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATAc                                               >  2:1259730/1‑97 (MQ=255)
                                                  attGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACa                                              <  2:2386051/97‑1 (MQ=255)
                                                  attGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACa                                              >  2:3041717/1‑97 (MQ=255)
                                                     gAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCa                                           >  1:411021/1‑97 (MQ=255)
                                                     gAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCa                                           >  1:2063415/1‑97 (MQ=255)
                                                     gAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCa                                           >  1:2987565/1‑97 (MQ=255)
                                                     gAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCa                                           >  2:2998120/1‑97 (MQ=255)
                                                       aaaaGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCaaa                                         <  1:372189/97‑1 (MQ=255)
                                                       aaaaGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCaaa                                         <  2:785902/97‑1 (MQ=255)
                                                         aaGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACg                                       >  1:3323663/1‑97 (MQ=255)
                                                         aaGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACg                                       >  2:2372192/1‑97 (MQ=255)
                                                         aaGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACg                                       <  1:2166256/97‑1 (MQ=255)
                                                         aaGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACg                                       >  2:515341/1‑97 (MQ=255)
                                                           ggTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGtt                                     <  2:1782473/97‑1 (MQ=255)
                                                            gTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGttt                                    >  1:938577/1‑97 (MQ=255)
                                                               aaaaTAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTAtt                                 >  2:1093248/1‑97 (MQ=255)
                                                                 aaTAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTg                               >  2:1250715/1‑97 (MQ=255)
                                                                  aTAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGa                              >  1:1154471/1‑97 (MQ=255)
                                                                  aTAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGa                              >  1:1499221/1‑97 (MQ=255)
                                                                   tAGGTCATATAATATAGTCAAAGCAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCAGAGAACAATACATCAAACGTTTATTTgca                             >  2:2154335/1‑95 (MQ=255)
                                                                   tAGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGaa                             >  2:3112691/1‑97 (MQ=255)
                                                                    aGGTCATATAATATAGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGCACAATACATCAAACGTTTATTTGAag                            >  1:2941722/1‑97 (MQ=255)
                                                                         atataatataGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCg                       >  1:2900001/1‑97 (MQ=255)
                                                                          tataatataGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGACGa                      >  2:1144160/1‑97 (MQ=255)
                                                                           ataatataGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGaa                     >  2:2604109/1‑97 (MQ=255)
                                                                            taatataGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTAGTTGAAGAGTCGCAt                    >  2:1736752/1‑97 (MQ=255)
                                                                            taatataGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAAt                    >  2:1438805/1‑97 (MQ=255)
                                                                            taatataGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAAt                    >  1:2594707/1‑97 (MQ=255)
                                                                             aatataGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATg                   <  1:1131619/97‑1 (MQ=255)
                                                                               tataGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGaa                 >  1:591248/1‑97 (MQ=255)
                                                                               tataGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGaa                 <  2:2058358/97‑1 (MQ=255)
                                                                                 taGTCAAAGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGa               >  2:2811455/1‑97 (MQ=255)
                                                                                   gTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATg             >  1:2208660/1‑97 (MQ=255)
                                                                                   gTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATg             <  2:1850708/97‑1 (MQ=255)
                                                                                     cAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTc           >  2:65172/1‑97 (MQ=255)
                                                                                      aaaGTAGTTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACTATACATCATACGTTTATTCGAATAGTCGAATGATGATGTCg          >  2:163061/1‑97 (MQ=255)
                                                                                      aaaGAAGGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCg          >  1:2954233/1‑97 (MQ=255)
                                                                                      aaaGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCg          >  1:640468/1‑97 (MQ=255)
                                                                                      aaaGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCg          >  1:2685808/1‑97 (MQ=255)
                                                                                       aagaagGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGt         >  2:1008578/1‑97 (MQ=255)
                                                                                         gaagGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTg       >  1:2482425/1‑97 (MQ=255)
                                                                                         gaagGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTg       >  2:2016162/1‑97 (MQ=255)
                                                                                         gaagGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTg       >  2:2301616/1‑97 (MQ=255)
                                                                                          aagGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTGa      <  1:2457782/97‑1 (MQ=255)
                                                                                          aagGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTGa      >  2:1105482/1‑97 (MQ=255)
                                                                                          aagGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTGa      <  1:2012770/97‑1 (MQ=255)
                                                                                           agGTCGAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGCTTATTTGAAGAGTCGAATGAAGATGTCGTTGaa     >  2:1385457/1‑97 (MQ=255)
                                                                                              tCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTGAAAtt  <  2:2963162/97‑1 (MQ=255)
                                                                                                |                                                                                              
AAAATGGAAGGGCTAACATCATAAATTTTTGAAAATATAAAAACAAATTAATTGAAAAAGGTCAAAATAGGTCATATAATATAGTCAAAGAAGGTCAAAAAGGGGTGATATACATGCACAATATGTCTGACATCATAGAACAATACATCAAACGTTTATTTGAAGAGTCGAATGAAGATGTCGTTGAAATT  >  BX571857/549710‑549900

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