Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,380,051 A→C 14.6% D73A (GAT→GCT)  SAS1289 → aconitate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,380,0510AC14.6% 615.8 / 111.2 203D73A (GAT→GCT) SAS1289aconitate hydratase
Reads supporting (aligned to +/- strand):  ref base A (107/66);  new base C (16/14);  total (123/80)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 206 total aligned reads displayed.

ATTCAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTTTACGTA  >  BX571857/1379955‑1380146
                                                                                                |                                                                                               
aTTCAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGa                                                                                                 >  2:2957599/1‑97 (MQ=255)
 ttCAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGAt                                                                                                >  1:3211356/1‑97 (MQ=255)
 ttCAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGAt                                                                                                <  1:3053663/97‑1 (MQ=255)
  tCAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGCTg                                                                                               >  2:1084379/1‑97 (MQ=255)
   cAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATgg                                                                                              >  1:1166629/1‑97 (MQ=255)
   cAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATgg                                                                                              >  1:1767551/1‑97 (MQ=255)
    aaTTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGa                                                                                             <  2:2636473/97‑1 (MQ=255)
    aaTTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGa                                                                                             <  2:2570640/97‑1 (MQ=255)
     aTTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGCTGGaa                                                                                            <  1:1987588/97‑1 (MQ=255)
       tCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAAt                                                                                          >  1:3030133/1‑97 (MQ=255)
       tCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAAt                                                                                          >  2:1517736/1‑97 (MQ=255)
        cGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATg                                                                                         <  2:488164/97‑1 (MQ=255)
        cGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATg                                                                                         <  2:434429/97‑1 (MQ=255)
         gtgtTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGa                                                                                        <  2:2272569/97‑1 (MQ=255)
          tgtTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGaa                                                                                       <  1:2283317/97‑1 (MQ=255)
            ttttGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAgg                                                                                     >  2:2678011/1‑97 (MQ=255)
              ttGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGCTGGAAATGAAGGCg                                                                                   >  1:2562341/1‑97 (MQ=255)
              ttGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCg                                                                                   >  2:3355639/1‑97 (MQ=255)
              ttGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCg                                                                                   >  1:2392272/1‑97 (MQ=255)
              ttGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCg                                                                                   >  1:2860557/1‑97 (MQ=255)
                  tAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGCTTGGAAAAGCTGGAAATGAAGGCGAGGt                                                                               >  1:1011552/1‑97 (MQ=255)
                    gAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTAc                                                                             <  1:345610/97‑1 (MQ=255)
                    gAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTAc                                                                             <  2:1177553/97‑1 (MQ=255)
                     aaTCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTAcc                                                                            <  2:2477186/97‑1 (MQ=255)
                      aTCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGCTGGAAATGAAGGCGAGGTACCa                                                                           <  2:1966133/97‑1 (MQ=255)
                      aTCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCa                                                                           >  2:787043/1‑97 (MQ=255)
                       tCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCCTATTAAAGCTTTACGTCAGTTTGGATACGCTGGAAATGAAGGCGAGGTACCAt                                                                          >  2:1536523/1‑97 (MQ=255)
                        cTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCAtt                                                                         >  2:2388217/1‑97 (MQ=255)
                         tttACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCAttt                                                                        >  1:1526210/1‑97 (MQ=255)
                          ttACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGCTGGAAATGAAGGCGAGGTACCATTTa                                                                       <  2:1442914/97‑1 (MQ=255)
                          ttACTTCGTCAAGAAGATGATATTGAAAATAAAGACGATCATATTAAAGCTTTAAGTCAGATTGGAACAGATGGAAATGAAGGCGTGGTACCactca                                                                       >  2:1495137/1‑93 (MQ=255)
                           tACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTaa                                                                      >  2:1337605/1‑97 (MQ=255)
                           tACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTaa                                                                      >  2:3275070/1‑97 (MQ=255)
                              ttCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAAcc                                                                   >  2:1712276/1‑97 (MQ=255)
                               tCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGCTGGAAATGAAGGCGAGGTACCATTTAAACCt                                                                  <  1:1879920/97‑1 (MQ=255)
                               tCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCt                                                                  >  1:2149536/1‑97 (MQ=255)
                               tCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGCTTGGCAAAGATGGAAATGAAGGCGAGGTACCATTTAAAccc                                                                  >  2:254430/1‑96 (MQ=255)
                                 gTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTtc                                                                <  2:2580798/97‑1 (MQ=255)
                                  tCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTtct                                                               <  1:2495024/97‑1 (MQ=255)
                                    aagaagATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACTTTATCg                                                             >  2:1824568/1‑97 (MQ=255)
                                    aagaagATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCg                                                             >  2:1564732/1‑97 (MQ=255)
                                    aagaagATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCg                                                             >  2:512483/1‑97 (MQ=255)
                                     agaagaTGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCTAGGTACCATTTAAACCTTCTCgt                                                            <  2:3010531/97‑1 (MQ=255)
                                     agaagaTGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCgt                                                            <  1:2104241/97‑1 (MQ=255)
                                      gaagaTGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCgtg                                                           >  2:356806/1‑97 (MQ=255)
                                       aagaTGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCgtgt                                                          >  1:179055/1‑97 (MQ=255)
                                       aagaTGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCgtgt                                                          <  1:1357806/97‑1 (MQ=255)
                                        agaTGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGtt                                                         >  2:2521545/1‑97 (MQ=255)
                                          atgatTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTAt                                                       >  1:2650472/1‑97 (MQ=255)
                                            gatTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTAttt                                                     <  1:429685/97‑1 (MQ=255)
                                               tttGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTAc                                                  >  2:1275658/1‑97 (MQ=255)
                                                ttGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACa                                                 >  1:1672739/1‑97 (MQ=255)
                                                ttGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACa                                                 >  1:2324318/1‑97 (MQ=255)
                                                ttGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACa                                                 >  2:295405/1‑97 (MQ=255)
                                                 tGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACaa                                                <  2:1229075/97‑1 (MQ=255)
                                                  gTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAg                                               <  1:516651/97‑1 (MQ=255)
                                                     aTTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGAtt                                            >  1:68711/1‑97 (MQ=255)
                                                     aTTACAGACGATCATATTAAAGCTTTAAGTCAGATTGGAAAAGATGGAAATGAAGGCGAGGTACCATTCAAACCTTCTCGTGTTATTTTATAAGAtg                                            >  1:775609/1‑96 (MQ=255)
                                                       tACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGCTGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTc                                          >  2:2163650/1‑97 (MQ=255)
                                                       tACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTc                                          <  2:479363/97‑1 (MQ=255)
                                                       tACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTc                                          >  2:392411/1‑97 (MQ=255)
                                                             cGATCATATTAAAGCTTTACGTCAGTTTGGAAAAGATGGAACTGAAGGCGAGGTCCCATTTAAACCTTCTCGTGTTATTTTACAAGATGTCACAGgt                                    >  2:3153011/1‑97 (MQ=255)
                                                             cGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGgt                                    >  2:956786/1‑97 (MQ=255)
                                                             cGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTCAAACATTCTCGTGTTATTTCACAAGATTTCACAggg                                    >  1:1659373/1‑96 (MQ=255)
                                                             cGATCATATTAAAGCTTTAAGCTTGTGTGCACAAGATGGAAATGAAGGCGAGGGCCCATTTACCCCTTCTCGTGTTATTATACACGATTTCACAGgt                                    >  2:2735922/1‑97 (MQ=255)
                                                                tCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTa                                 <  2:3030133/97‑1 (MQ=255)
                                                                 cATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTAc                                >  2:2819240/1‑97 (MQ=255)
                                                                 cATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTAc                                >  2:839998/1‑97 (MQ=255)
                                                                  atatTAAAGCTTTAAGTCAGTTTGGAAAAGCTGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTAcc                               >  1:2033847/1‑97 (MQ=255)
                                                                  atatTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTAcc                               >  1:160598/1‑97 (MQ=255)
                                                                   tatTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCa                              <  1:743322/97‑1 (MQ=255)
                                                                     ttAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGc                            >  1:2904002/1‑97 (MQ=255)
                                                                        aaGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGt                         >  1:2846225/1‑97 (MQ=255)
                                                                        aaGCTTTAAGCCAGTATGGAAAACCAGGAAATGAAGGCGAGGTACCATTTAAACCCTCTCGTGCTATCTTACAAGATTTCCCCGGCGTACCAGCCgc                         >  2:815669/1‑96 (MQ=255)
                                                                         aGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTa                        >  2:119399/1‑97 (MQ=255)
                                                                         aGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTa                        <  2:1199444/97‑1 (MQ=255)
                                                                          gCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAg                       >  1:143084/1‑97 (MQ=255)
                                                                           cTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGt                      >  2:1051652/1‑97 (MQ=255)
                                                                            tttAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGtt                     <  1:1942423/97‑1 (MQ=255)
                                                                               aaGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGAt                  >  1:631261/1‑97 (MQ=255)
                                                                                aGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGAtt                 >  1:2297391/1‑97 (MQ=255)
                                                                                 gTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAATCCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGAttt                <  1:97018/97‑1 (MQ=255)
                                                                                 gTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGAttt                <  1:1136353/97‑1 (MQ=255)
                                                                                  tCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTa               >  2:2298240/1‑97 (MQ=255)
                                                                                  tCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTa               >  2:2978570/1‑97 (MQ=255)
                                                                                   cAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAg              >  2:1766872/1‑97 (MQ=255)
                                                                                   cAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAg              <  2:826048/97‑1 (MQ=255)
                                                                                    aGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGc             >  1:986928/1‑97 (MQ=255)
                                                                                    aGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGc             <  1:2752471/97‑1 (MQ=255)
                                                                                    aGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGc             >  2:541996/1‑97 (MQ=255)
                                                                                       ttGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGcttc          >  1:2961545/1‑97 (MQ=255)
                                                                                        tGGAAAAGCTGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGcttct         >  1:941265/1‑97 (MQ=255)
                                                                                        tGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGcttct         >  2:2463246/1‑97 (MQ=255)
                                                                                         ggAAAAGAAGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGcttctt        <  1:998809/97‑1 (MQ=255)
                                                                                              aGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTTTACGt   >  2:3041223/1‑97 (MQ=255)
                                                                                              aGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTTTACGt   >  1:2806036/1‑97 (MQ=255)
                                                                                              aGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTTTACGt   >  1:2750071/1‑97 (MQ=255)
                                                                                              aGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGATATTTTACAAGATATCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTGTACGt   >  1:2300965/1‑97 (MQ=255)
                                                                                               gATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTTTACGTa  >  2:1046194/1‑97 (MQ=255)
                                                                                               gATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTTTACGTa  <  2:1866598/97‑1 (MQ=255)
                                                                                                |                                                                                               
ATTCAATTCGTGTTTTGTTAGAATCTTTACTTCGTCAAGAAGATGATTTTGTAATTACAGACGATCATATTAAAGCTTTAAGTCAGTTTGGAAAAGATGGAAATGAAGGCGAGGTACCATTTAAACCTTCTCGTGTTATTTTACAAGATTTCACAGGTGTACCAGCCGTAGTTGATTTAGCTTCTTTACGTA  >  BX571857/1379955‑1380146

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: ‑