Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,169,390 T→C 10.6% G206G (GGA→GGG SAS1097 ← alpha‑hemolysin precursor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,169,3900TC10.6% 864.1 / 70.4 232G206G (GGA→GGGSAS1097alpha‑hemolysin precursor
Reads supporting (aligned to +/- strand):  ref base T (119/88);  new base C (12/13);  total (131/101)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 234 total aligned reads displayed.

ATCAAGGAAGTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTGAAC  >  BX571857/1169294‑1169485
                                                                                                |                                                                                               
 tCAAGGAAGTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTc                                                                                                <  2:2866059/97‑1 (MQ=255)
    aGGAAGTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTcccc                                                                                             <  2:63330/97‑1 (MQ=255)
     ggAAGTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCa                                                                                            <  1:1322000/97‑1 (MQ=255)
     ggAAGTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCa                                                                                            >  2:3542228/1‑97 (MQ=255)
      gAAGTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCaa                                                                                           <  1:2377590/97‑1 (MQ=255)
       aaGTTATCTGCTGCTTTCATAGAACCATTTCTCGTTTTCATGAAACGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAAt                                                                                          >  2:214638/1‑97 (MQ=255)
         gTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAAttt                                                                                        >  2:2937895/1‑97 (MQ=255)
         gTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAAttt                                                                                        >  1:1460470/1‑97 (MQ=255)
           tATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTg                                                                                      <  1:3317878/97‑1 (MQ=255)
           tATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTg                                                                                      >  1:282529/1‑97 (MQ=255)
           tATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTg                                                                                      >  2:1552310/1‑97 (MQ=255)
           tATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTg                                                                                      <  2:2133724/97‑1 (MQ=255)
           tATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTg                                                                                      >  1:1436511/1‑97 (MQ=255)
            aTCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGa                                                                                     <  1:334982/97‑1 (MQ=255)
               tgctgcTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTc                                                                                  >  1:1112141/1‑97 (MQ=255)
               tgctgcTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTc                                                                                  >  1:2016329/1‑97 (MQ=255)
               tgatgctTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTc                                                                                  >  2:2381482/4‑97 (MQ=255)
                gctgctTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCa                                                                                 <  1:1658403/97‑1 (MQ=255)
                  tgctTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCAcc                                                                               >  2:1068046/1‑97 (MQ=255)
                  tgctTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCAcc                                                                               >  2:2805246/1‑97 (MQ=255)
                   gctTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCa                                                                              <  1:3148514/97‑1 (MQ=255)
                     tttCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCata                                                                            <  1:2973739/97‑1 (MQ=255)
                      ttCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCatat                                                                           >  1:2594903/1‑97 (MQ=255)
                       tCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATAtt                                                                          >  1:1597839/1‑97 (MQ=255)
                        cATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTg                                                                         >  2:3557008/1‑97 (MQ=255)
                          tAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGtt                                                                       <  2:2702399/97‑1 (MQ=255)
                           aGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTa                                                                      >  2:2654834/1‑97 (MQ=255)
                            gAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTaa                                                                     >  2:2888320/1‑97 (MQ=255)
                            gAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTaa                                                                     >  2:1199257/1‑97 (MQ=255)
                            gAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTaa                                                                     <  2:2283773/97‑1 (MQ=255)
                            gAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTaa                                                                     >  1:2255888/1‑97 (MQ=255)
                             aaCCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTaaa                                                                    <  1:3478659/97‑1 (MQ=255)
                             aaCCATTCCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTaaa                                                                    >  1:2413984/1‑97 (MQ=255)
                              aCCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATCCCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAat                                                                   <  1:1881232/97‑1 (MQ=255)
                              aCCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAat                                                                   <  1:2181965/97‑1 (MQ=255)
                              aCCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAat                                                                   <  1:3632347/97‑1 (MQ=255)
                              aCCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAat                                                                   <  2:200809/97‑1 (MQ=255)
                               ccATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAata                                                                  >  2:1800173/1‑97 (MQ=255)
                               ccATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAata                                                                  <  1:3012227/97‑1 (MQ=255)
                               ccATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAata                                                                  <  2:3185043/97‑1 (MQ=255)
                                cATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAatat                                                                 >  2:8670/1‑97 (MQ=255)
                                cATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAatat                                                                 >  1:32011/1‑97 (MQ=255)
                                    tCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACt                                                             >  1:3445105/1‑97 (MQ=255)
                                    tCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACt                                                             >  1:3245609/1‑97 (MQ=255)
                                    tCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAATATCACt                                                             >  2:3289623/1‑97 (MQ=255)
                                       aGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTc                                                          <  1:1871397/97‑1 (MQ=255)
                                               tGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTAc                                                  <  1:1207721/97‑1 (MQ=255)
                                                gAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACt                                                 >  2:3376858/1‑97 (MQ=255)
                                                 aaaaGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACtt                                                >  2:784800/1‑97 (MQ=255)
                                                 aaaaGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACtt                                                >  1:1709383/1‑97 (MQ=255)
                                                 aaaaGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACtt                                                >  2:2432298/1‑97 (MQ=255)
                                                   aaGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACtttt                                              <  2:3394022/97‑1 (MQ=255)
                                                    aGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACttttt                                             <  1:1187007/97‑1 (MQ=255)
                                                    aGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACttttt                                             <  1:1935811/97‑1 (MQ=255)
                                                     gTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACtttttt                                            <  1:936903/97‑1 (MQ=255)
                                                     gTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACtttttt                                            <  1:1822807/97‑1 (MQ=255)
                                                       tGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTAt                                          >  2:744825/1‑97 (MQ=255)
                                                       tGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTAt                                          >  2:2613807/1‑97 (MQ=255)
                                                        gATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATc                                         <  1:1383587/97‑1 (MQ=255)
                                                           tGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGt                                      >  1:3657546/1‑97 (MQ=255)
                                                            gCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGtt                                     >  2:3481699/1‑97 (MQ=255)
                                                            gCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGtt                                     >  2:1389045/1‑97 (MQ=255)
                                                            gCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGtt                                     <  2:445823/97‑1 (MQ=255)
                                                            gCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAATATCCCTTTCCAGCCTACTTTTTTATCAGtt                                     >  2:1721986/1‑97 (MQ=255)
                                                            gCCATATACCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATCGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGtt                                     <  2:1341373/97‑1 (MQ=255)
                                                             ccATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTg                                    <  2:42055/97‑1 (MQ=255)
                                                                tataCCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGc                                 <  2:1896781/97‑1 (MQ=255)
                                                                tataCCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGc                                 <  2:1143334/97‑1 (MQ=255)
                                                                 ataCCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGct                                <  1:2843574/97‑1 (MQ=255)
                                                                  taCCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTAAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGctc                               >  1:2926273/1‑97 (MQ=255)
                                                                  taCCGGGTTCCAAGAATCTCTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGctc                               <  1:1721335/97‑1 (MQ=255)
                                                                   aCCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGctct                              >  1:54280/1‑97 (MQ=255)
                                                                       ggTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTaa                          <  2:3035923/97‑1 (MQ=255)
                                                                       ggTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTaa                          >  1:1562890/1‑97 (MQ=255)
                                                                         ttCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTaaaa                        <  2:2517573/97‑1 (MQ=255)
                                                                         gtCCAAGAATCACTATCATATGGCCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTaaaa                        >  2:3144507/2‑97 (MQ=255)
                                                                          tCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAAt                       <  2:1082176/97‑1 (MQ=255)
                                                                          tCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAAt                       >  1:268301/1‑97 (MQ=255)
                                                                           ccAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAAtt                      <  1:375869/97‑1 (MQ=255)
                                                                           ccAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAAtt                      >  1:2109126/1‑97 (MQ=255)
                                                                             aaGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGt                    <  1:1015291/97‑1 (MQ=255)
                                                                                aaTCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGtttt                 >  2:1915379/1‑97 (MQ=255)
                                                                                 aTCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTg                >  1:2000838/1‑97 (MQ=255)
                                                                                  tctcTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGa               <  2:877558/97‑1 (MQ=255)
                                                                                  tctcTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGa               <  1:353828/97‑1 (MQ=255)
                                                                                  tctcTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGa               >  1:1125525/1‑97 (MQ=255)
                                                                                      tATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATc           <  1:783846/97‑1 (MQ=255)
                                                                                       aTCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCa          <  2:3633958/97‑1 (MQ=255)
                                                                                       aTCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCa          <  1:2598199/97‑1 (MQ=255)
                                                                                        tCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAg         >  1:896689/1‑97 (MQ=255)
                                                                                         cATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAgg        <  2:1687035/97‑1 (MQ=255)
                                                                                         cATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAgg        <  1:3368893/97‑1 (MQ=255)
                                                                                          atatGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGt       >  2:47791/1‑97 (MQ=255)
                                                                                          atatGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGt       >  1:2196928/1‑97 (MQ=255)
                                                                                            atGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTg     <  1:428629/97‑1 (MQ=255)
                                                                                             tGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTGa    <  2:2963196/97‑1 (MQ=255)
                                                                                             tGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTGa    >  2:1578115/1‑97 (MQ=255)
                                                                                              ggTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTGaa   <  1:2753617/97‑1 (MQ=255)
                                                                                               gTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTGAAc  >  2:1002605/1‑97 (MQ=255)
                                                                                               gTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTGAAc  >  1:2412422/1‑97 (MQ=255)
                                                                                                |                                                                                               
ATCAAGGAAGTTATCTGCTGCTTTCATAGAACCATTTCTAGTTTTCATGAAAAGTTGATTGCCATATACCGGGTTCCAAGAATCTCTATCATATGGTCCCCAATTTTGATTCACCATATTGTTAAATATCACTTTCCAGCCTACTTTTTTATCAGTTGGGCTCTCTAAAATTGTTTTGAAATCAGGTTGAAC  >  BX571857/1169294‑1169485

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