Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 4,422 C→A 5.1% Q171K (CAA→AAA)  SAS0004 → DNA replication and repair protein RecF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718574,4220CA 5.1% 1001.7 / 3.3 250Q171K (CAA→AAA) SAS0004DNA replication and repair protein RecF
Reads supporting (aligned to +/- strand):  ref base C (137/95);  new base A (2/16);  total (139/111)
Fisher's exact test for biased strand distribution p-value = 8.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

GGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGC  >  BX571857/4328‑4518
                                                                                              |                                                                                                
    cAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCaaaaaa                                                                                            >  1:2305012/1‑97 (MQ=255)
     aaaTTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCaaaaaaa                                                                                           >  2:91275/1‑97 (MQ=255)
      aaTTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAg                                                                                          <  2:1193927/97‑1 (MQ=255)
         ttCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGAc                                                                                       >  1:1702924/1‑97 (MQ=255)
         ttCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGACAAAAAAAGGAc                                                                                       <  2:2365918/97‑1 (MQ=255)
          tCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACt                                                                                      >  2:2118931/1‑97 (MQ=255)
           ctgctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACtt                                                                                     >  1:2124594/1‑97 (MQ=255)
           ctgctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCAAAAAAAAGGACtt                                                                                     <  2:319053/97‑1 (MQ=255)
            tgctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTa                                                                                    <  2:2245264/97‑1 (MQ=255)
              ctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAAc                                                                                  >  2:1391228/1‑97 (MQ=255)
              ctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAAc                                                                                  <  1:891370/97‑1 (MQ=255)
                   tACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGt                                                                             >  2:649412/1‑97 (MQ=255)
                   tACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGt                                                                             >  2:1200962/1‑97 (MQ=255)
                    aCTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGtt                                                                            >  1:750234/1‑97 (MQ=255)
                      ttAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTgg                                                                          <  1:518517/97‑1 (MQ=255)
                      ttAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTgg                                                                          >  1:1742976/1‑97 (MQ=255)
                       tAAATGATATAGCTAAATACCAACGTATTATAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGAAAAAAAAAGGACTTAACAATGTTGGa                                                                         <  2:1161998/97‑1 (MQ=255)
                        aaaTGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGaa                                                                        >  2:442799/1‑97 (MQ=255)
                         aaTGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAg                                                                       <  1:546254/97‑1 (MQ=255)
                         aaTGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAg                                                                       <  1:1947852/97‑1 (MQ=255)
                         aaTGACTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAAAAAAAAAAAAAAAAGGACTTAACAATGTTGGAAg                                                                       <  2:864852/97‑1 (MQ=255)
                          aTGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGt                                                                      >  1:1228688/1‑97 (MQ=255)
                          aTGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAAAAAAAAAAAAGGACTTAACAATGTTGGAAGt                                                                      <  2:1164325/97‑1 (MQ=255)
                           tGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTa                                                                     >  2:1545629/1‑97 (MQ=255)
                             aTTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTAtt                                                                   >  2:2430627/1‑97 (MQ=255)
                                tAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTaaa                                                                >  2:2108208/1‑97 (MQ=255)
                                    tCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcag                                                            >  1:112538/1‑97 (MQ=255)
                                     cAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagc                                                           <  1:455772/97‑1 (MQ=255)
                                      aaTACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagca                                                          <  1:782572/97‑1 (MQ=255)
                                      aaTACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagca                                                          >  1:1575837/1‑97 (MQ=255)
                                       aTACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagcag                                                         <  2:2170940/97‑1 (MQ=255)
                                        tACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGt                                                        <  2:2331625/97‑1 (MQ=255)
                                         aCCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                       <  1:867214/97‑1 (MQ=255)
                                         aCCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                       >  1:29934/1‑97 (MQ=255)
                                         aCCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                       >  1:283201/1‑97 (MQ=255)
                                         aCCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                       <  1:1870508/97‑1 (MQ=255)
                                          ccAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGttt                                                      >  2:2381284/1‑97 (MQ=255)
                                          ccAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATGAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGttt                                                      >  1:1082648/1‑97 (MQ=255)
                                              cGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTg                                                  <  2:1267401/97‑1 (MQ=255)
                                              cGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTg                                                  <  1:2249737/97‑1 (MQ=255)
                                              cGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTg                                                  <  2:2056308/97‑1 (MQ=255)
                                              cGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTg                                                  >  1:1350045/1‑97 (MQ=255)
                                               gTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGa                                                 <  1:2418068/97‑1 (MQ=255)
                                                tATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGaa                                                >  1:1660054/1‑97 (MQ=255)
                                                tATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGaa                                                <  1:633914/97‑1 (MQ=255)
                                                tATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGaa                                                >  1:705548/1‑97 (MQ=255)
                                                tATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGaa                                                >  2:1845235/1‑97 (MQ=255)
                                                 aTTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAAATAGGCAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAat                                               <  2:754172/97‑1 (MQ=255)
                                                  ttttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAata                                              >  2:2410236/1‑97 (MQ=255)
                                                  ttttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAata                                              >  2:1988469/1‑97 (MQ=255)
                                                  ttttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAata                                              >  2:37446/1‑97 (MQ=255)
                                                  ttttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAAGACTTAACAATGTTGGAAGAATTAAATCAGCAGTTAGCTGAata                                              >  1:143775/1‑97 (MQ=255)
                                                    ttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATg                                            >  1:69177/1‑97 (MQ=255)
                                                     tAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGc                                           >  2:1058610/1‑97 (MQ=255)
                                                     tAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGc                                           <  1:1275302/97‑1 (MQ=255)
                                                        aGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAAt                                        >  2:1419830/1‑97 (MQ=255)
                                                         gCAAAAGAATAATTACTTAAAGCATTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATg                                       <  2:1468774/97‑1 (MQ=255)
                                                         gCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATg                                       >  2:155887/1‑97 (MQ=255)
                                                         gCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATg                                       >  1:2359204/1‑97 (MQ=255)
                                                            aaaGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGaaa                                    <  2:1668428/97‑1 (MQ=255)
                                                            aaaGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGaaa                                    >  2:96543/1‑97 (MQ=255)
                                                             aaGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATGAGCAATGAAAg                                   >  2:1326606/1‑97 (MQ=255)
                                                              aGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTGGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGt                                  <  1:1050213/97‑1 (MQ=255)
                                                               gAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTa                                 <  2:271749/97‑1 (MQ=255)
                                                                aataatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTaa                                >  1:778290/1‑97 (MQ=255)
                                                                 ataatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAAc                               >  1:1455738/1‑97 (MQ=255)
                                                                   aatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTg                             <  2:1713714/97‑1 (MQ=255)
                                                                   aatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTg                             >  1:1558428/1‑97 (MQ=255)
                                                                   aatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTg                             >  1:2167506/1‑97 (MQ=255)
                                                                       aCTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATaa                         >  2:2046201/1‑97 (MQ=255)
                                                                        cTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATaaa                        <  1:1804867/97‑1 (MQ=255)
                                                                         ttAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAAc                       >  1:1978575/1‑97 (MQ=255)
                                                                         ttAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAAc                       <  1:404557/97‑1 (MQ=255)
                                                                           aaaGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGt                     >  2:593160/1‑97 (MQ=255)
                                                                            aaGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTg                    <  1:2196918/97‑1 (MQ=255)
                                                                              gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGca                  >  1:1609392/1‑97 (MQ=255)
                                                                              gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGca                  <  2:2078156/97‑1 (MQ=255)
                                                                              gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGca                  <  1:981588/97‑1 (MQ=255)
                                                                              gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGca                  <  2:1619373/97‑1 (MQ=255)
                                                                              gCAGTGACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGca                  <  1:1623200/97‑1 (MQ=255)
                                                                               cAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGcac                 <  1:710590/97‑1 (MQ=255)
                                                                               cAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGcac                 <  1:275610/97‑1 (MQ=255)
                                                                                   tACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACAttt             <  1:725886/97‑1 (MQ=255)
                                                                                   tACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACAttt             <  2:567148/97‑1 (MQ=255)
                                                                                    aCAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACAtttt            >  2:491355/1‑97 (MQ=255)
                                                                                    aCAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACAtttt            >  2:668520/1‑97 (MQ=255)
                                                                                    aCAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACAtttt            <  1:200107/97‑1 (MQ=255)
                                                                                     cAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTa           <  1:74068/97‑1 (MQ=255)
                                                                                     cAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTa           >  1:1367473/1‑97 (MQ=255)
                                                                                      aaTTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTAt          <  2:52264/97‑1 (MQ=255)
                                                                                      aaTTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTAt          >  2:1778493/1‑97 (MQ=255)
                                                                                      aaTTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTAt          <  1:1028470/97‑1 (MQ=255)
                                                                                       aTTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGCAACTGATAAACGTGCACATTTTAtt         >  2:718605/1‑97 (MQ=255)
                                                                                        ttAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATACACGTGCACATTTTATTc        >  1:341308/1‑97 (MQ=255)
                                                                                        ttAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTc        >  2:807804/1‑97 (MQ=255)
                                                                                           ggCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAag     <  2:176382/97‑1 (MQ=255)
                                                                                            gCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAaga    <  1:2148418/97‑1 (MQ=255)
                                                                                            gCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAaga    <  1:1491481/97‑1 (MQ=255)
                                                                                              cAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGc  <  2:1101138/97‑1 (MQ=255)
                                                                                              cAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGc  >  2:841716/1‑97 (MQ=255)
                                                                                              |                                                                                                
GGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGC  >  BX571857/4328‑4518

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: