Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,434,357 C→G 12.0% A141G (GCT→GGT)  SAS1338 → putative tetrahydrodipicolinate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,434,3570CG12.0% 764.7 / 23.2 220A141G (GCT→GGT) SAS1338putative tetrahydrodipicolinate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (93/98);  new base G (26/2);  total (119/101)
Fisher's exact test for biased strand distribution p-value = 3.83e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

TGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGT  >  BX571857/1434262‑1434453
                                                                                               |                                                                                                
tGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCt                                                                                                 >  1:1619232/1‑97 (MQ=255)
 ggcgCTGTCGTATTGATGGGCGCAACGATGAGTATGGGGGCAGGCGTTTGCGAAGGTACGATGTTTGATATGAAGGCTACTCTCGGTGGTCGTGCTa                                                                                                <  2:723061/94‑1 (MQ=255)
  gTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTAc                                                                                               >  1:1456929/1‑97 (MQ=255)
    gCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACaa                                                                                             <  1:399325/97‑1 (MQ=255)
     cTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAAc                                                                                            >  1:2951759/1‑97 (MQ=255)
     cTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAAc                                                                                            >  1:2952070/1‑97 (MQ=255)
         cGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                        >  1:2910533/1‑97 (MQ=255)
         cGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                        <  1:2445100/97‑1 (MQ=255)
         cGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                        >  1:1841909/1‑97 (MQ=255)
          gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                       >  1:443615/1‑97 (MQ=255)
          gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACAACTAGTa                                                                                       >  2:577661/1‑97 (MQ=255)
            tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                                     <  1:2775753/97‑1 (MQ=255)
            tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                                     <  1:816535/97‑1 (MQ=255)
            tATGATGGGCGCAACAATTAATATTGGCGCAGCCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                                     >  1:3001942/1‑97 (MQ=255)
             atgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACCACTGGAaaaa                                                                                    >  2:1774634/1‑97 (MQ=255)
             atgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                                    <  1:660969/97‑1 (MQ=255)
             atgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                                    >  2:2908514/1‑97 (MQ=255)
              tgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATACGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                                   <  1:2789057/97‑1 (MQ=255)
                atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                 >  1:48250/1‑97 (MQ=255)
                atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                 >  1:959321/1‑97 (MQ=255)
                atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCCTACCACTGGaaaaaacg                                                                                 >  2:467960/1‑95 (MQ=255)
                atgGGCGCAACAATGAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGCTCGTGCTACAACTGGTAAAAATg                                                                                 <  2:3006055/97‑1 (MQ=255)
                  gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTa                                                                               <  1:2615323/97‑1 (MQ=255)
                  gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTAGAAAAAAAGTa                                                                               >  2:1393165/1‑97 (MQ=255)
                           aaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                      >  1:1698219/1‑97 (MQ=255)
                           aaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                      >  1:3335779/1‑97 (MQ=255)
                              tAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGcg                                                                   >  1:2392752/1‑96 (MQ=255)
                               aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTg                                                                  >  1:1792229/1‑97 (MQ=255)
                               aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTg                                                                  >  1:1175089/1‑97 (MQ=255)
                               aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCACTGAAAAAAAAGTACATGTAGGGGCTg                                                                  >  2:1410810/1‑97 (MQ=255)
                                  atTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCCTACCACTAGAAAAAAAGTACATGTAGGGGCTGgcg                                                               >  2:1797428/1‑97 (MQ=255)
                                   ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCACTAAAAAAAAAATACATGTAGGGGCTGgcga                                                              >  2:123072/1‑96 (MQ=255)
                                    tGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCa                                                             >  1:3260645/1‑97 (MQ=255)
                                     ggCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGCCTGGCGCAg                                                            >  1:902612/1‑97 (MQ=255)
                                      gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGt                                                           <  1:252304/97‑1 (MQ=255)
                                       cgcAGGCGTTGGCGAAGGTACAATGATTGGTATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                          <  2:791864/97‑1 (MQ=255)
                                           gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                      >  1:589886/1‑97 (MQ=255)
                                           gTCGTGGGCGAAGGGACAATGATTGATATGAATGCAACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                      <  2:64382/97‑1 (MQ=255)
                                             cGTTGGCGAAGGTACAATGATTGATATGAATGCTGCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                    <  2:1253283/97‑1 (MQ=255)
                                             cGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                    <  2:3219506/97‑1 (MQ=255)
                                             cGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                    >  1:1971550/1‑97 (MQ=255)
                                              gTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAg                                                   <  1:2426784/97‑1 (MQ=255)
                                               ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAgg                                                  <  1:2886862/97‑1 (MQ=255)
                                                tGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                 <  2:329787/97‑1 (MQ=255)
                                                 ggCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                                >  1:1001174/1‑97 (MQ=255)
                                                 ggCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCATTGGAAAAAATATACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                                >  2:1689411/1‑97 (MQ=255)
                                                    gAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                             >  2:2453234/1‑97 (MQ=255)
                                                      aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAtt                                           <  1:2029247/97‑1 (MQ=255)
                                                      aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAtt                                           >  1:2839615/1‑97 (MQ=255)
                                                      aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAtt                                           >  1:2679563/1‑97 (MQ=255)
                                                       ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                          <  1:2176900/97‑1 (MQ=255)
                                                       ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGGGaggg                                          >  2:511875/1‑94 (MQ=255)
                                                        gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                         >  1:914960/1‑97 (MQ=255)
                                                        gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                         >  1:562452/1‑97 (MQ=255)
                                                        gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                         >  1:1931130/1‑97 (MQ=255)
                                                        gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                         <  1:2240424/97‑1 (MQ=255)
                                                         tACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                        >  1:1249460/1‑97 (MQ=255)
                                                         tACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                        <  1:2316751/97‑1 (MQ=255)
                                                           cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACCACTAATAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                      >  2:521645/1‑97 (MQ=255)
                                                            aaTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccc                                     <  1:1512373/97‑1 (MQ=255)
                                                            aaTGAGTGATATGAATGCTGCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccc                                     <  2:2263217/97‑1 (MQ=255)
                                                             aTGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTGAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                    >  2:281923/1‑97 (MQ=255)
                                                              tgattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccccc                                   >  1:1278118/1‑97 (MQ=255)
                                                               gattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACACCCt                                  >  1:775705/1‑97 (MQ=255)
                                                                attgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTa                                 <  1:1563298/97‑1 (MQ=255)
                                                                 ttgatATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAg                                >  2:2571667/1‑97 (MQ=255)
                                                                  tgatATGAATGCTACTCTCGGTGGTCGTGGTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                               <  2:930113/97‑1 (MQ=255)
                                                                  tgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                               >  1:1207144/1‑97 (MQ=255)
                                                                  tgatATGAATGCTACTCTCGGTGGTCGCCGTACCACTGGAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                               >  2:2330376/1‑97 (MQ=255)
                                                                   gatATGAATGCTACTCTCGGTGGTCGCCGTACCACTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                              >  2:2285960/1‑97 (MQ=255)
                                                                    atatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGc                             <  2:1902263/97‑1 (MQ=255)
                                                                    atatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGc                             <  1:1543288/97‑1 (MQ=255)
                                                                       tGAATGCTACTCTCGGTGGTCGCCGTACCACTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATc                          >  2:2726195/1‑97 (MQ=255)
                                                                        gAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCa                         <  1:1105731/97‑1 (MQ=255)
                                                                         aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCAc                        >  1:1779856/1‑97 (MQ=255)
                                                                           tGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGACCCCCCTAGTGCATCACCg                      <  2:2682254/97‑1 (MQ=255)
                                                                           tGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCg                      >  2:1339235/1‑97 (MQ=255)
                                                                            gCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCgg                     >  1:2007506/1‑97 (MQ=255)
                                                                            gCTACTCTCGGTGGTCGCCGTACCACTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCgg                     >  2:2181917/1‑97 (MQ=255)
                                                                               aCTCTCGGTGGTCGTGCTACAGCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTta                  <  2:1577511/97‑1 (MQ=255)
                                                                               aCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTta                  >  1:3298794/1‑97 (MQ=255)
                                                                               aCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTta                  <  1:304660/97‑1 (MQ=255)
                                                                                ctctCGGTGGTCGTGCTACAGCTGGTAAAAATGTACATGGAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTtat                 <  2:767914/97‑1 (MQ=255)
                                                                                ctctCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTtat                 >  1:4540/1‑97 (MQ=255)
                                                                                 tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGGAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTtata                <  2:2412995/97‑1 (MQ=255)
                                                                                     ggtggtCGCCCTACCACTAAAAAAAAAGCACATGTAGGGGCTGGCGCAGTATTAGCACGTGTGACTGAACCCCCTAGTGCATCACCGGTTATAATCg            >  2:1898072/1‑97 (MQ=255)
                                                                                     ggtggtCGCCCTACAACTAGAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCg            >  2:1969691/1‑97 (MQ=255)
                                                                                     ggtgggCGTGTTACAACTGGGAAAAATGTACATGTAGGGGCGGGCGCAGTATTAGCAGGGGTGATTGAGGCCCCGAGTGCATCACCGGTTATAATCg            <  2:2953277/97‑1 (MQ=255)
                                                                                       tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGATCCCCCTAGTGCATCACCGGTTATAATCGAg          <  2:2499305/97‑1 (MQ=255)
                                                                                       tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAg          <  1:2084476/97‑1 (MQ=255)
                                                                                        ggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCgggg         <  1:3159560/97‑1 (MQ=255)
                                                                                          tCGTGGTAAAACTGGTAAAAATGTACGTGTAGGGGCTGGGGCAGTATTAGCAGGTGGGATTGAGCCCCCTAGTGCATAACCGGTTATAATCGAGgat       <  2:843612/97‑1 (MQ=255)
                                                                                          tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGgat       <  1:2855331/97‑1 (MQ=255)
                                                                                           cGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGgctg      >  2:1129952/1‑97 (MQ=255)
                                                                                           cGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGgatg      >  2:876275/1‑97 (MQ=255)
                                                                                              gCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGatgatg   <  1:271504/97‑1 (MQ=255)
                                                                                               cTACAACTGGTAAAAATGTACGTGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGACCCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGt  <  2:1893588/97‑1 (MQ=255)
                                                                                               cTACAACTGGTAAAAATGTACGTGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGt  <  2:2541840/97‑1 (MQ=255)
                                                                                               cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGt  >  1:3087228/1‑97 (MQ=255)
                                                                                               cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGt  >  2:3076062/1‑97 (MQ=255)
                                                                                               |                                                                                                
TGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCATCACCGGTTATAATCGAGGATGATGT  >  BX571857/1434262‑1434453

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