Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,725,352 A→T 12.1% intergenic (‑58/+276) USA300HOU_RS14015 ← / ← USA300HOU_RS14020 Baeyer‑Villiger flavin‑containing monooxygenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,725,3520AT12.1% 730.2 / 69.5 240intergenic (‑58/+276)USA300HOU_RS14015/USA300HOU_RS14020Baeyer‑Villiger flavin‑containing monooxygenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (104/106);  new base T (11/18);  total (115/125)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

GGCACCGATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTT  >  USA300TCH1516_ALE/2725262‑2725441
                                                                                          |                                                                                         
 gcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaataa                                                                                          >  1:442268/1‑91 (MQ=255)
  caccgATAATAGCCACTTTATGATGTTGCATATTAACGCTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTaattat                                                                                         <  2:2247635/91‑1 (MQ=255)
  caccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTaataat                                                                                         <  2:1859118/91‑1 (MQ=255)
  caccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTaataat                                                                                         >  1:2088414/1‑91 (MQ=255)
   accgATAGTAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAAtt                                                                                        <  2:382208/91‑1 (MQ=255)
   accgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAattatt                                                                                        <  2:3495138/91‑1 (MQ=255)
    ccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTa                                                                                       >  1:1309152/1‑91 (MQ=255)
    ccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTa                                                                                       <  1:3390113/91‑1 (MQ=255)
    ccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCCATTTCAAGTTATTGTTCTTGATTATAAATCGTAATAATTa                                                                                       >  1:596848/1‑91 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                      <  2:1352250/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTAc                                                                                      >  2:1210756/1‑91 (MQ=255)
      gATAATAGGCACTTTATGATGTTGCATATTAAGAGGTCCTCAATGATTTAGTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACg                                                                                     <  2:2075921/91‑1 (MQ=255)
      gATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACg                                                                                     <  2:1915422/91‑1 (MQ=255)
      gATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACg                                                                                     <  2:3039791/91‑1 (MQ=255)
       ataataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGa                                                                                    <  2:3603468/91‑1 (MQ=255)
        taataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAt                                                                                   >  1:2097930/1‑91 (MQ=255)
        taataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAt                                                                                   >  1:1711418/1‑91 (MQ=255)
        taataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAACTACGAt                                                                                   <  1:3471801/91‑1 (MQ=255)
         aataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGAtt                                                                                  <  1:646200/91‑1 (MQ=255)
         aataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAtt                                                                                  <  2:3333654/91‑1 (MQ=255)
           tggCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTTTTTTTGATTATAAATCGTAATAATTACGATTta                                                                                <  2:178707/89‑1 (MQ=255)
            aGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGATTtat                                                                               <  1:531250/91‑1 (MQ=255)
            aGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGATTtat                                                                               <  1:3478740/91‑1 (MQ=255)
             gCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTtata                                                                              <  2:1164562/91‑1 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                             <  1:3917318/91‑1 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                             <  1:3308165/91‑1 (MQ=255)
               cACTTTATGATGTTGCATATTACCACTCCCTCAATGCTCTCTTTCACGTCATTGTTTTTGATTATTAATCGTACTAATTACGATTTATaaa                                                                            >  2:1556309/1‑91 (MQ=255)
               cACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCCAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaaa                                                                            >  2:1679528/1‑91 (MQ=255)
               cACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaaa                                                                            >  1:671761/1‑91 (MQ=255)
                aCTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAg                                                                           <  1:1070361/91‑1 (MQ=255)
                aCTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAg                                                                           >  1:1909992/1‑91 (MQ=255)
                   ttATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGtat                                                                        >  1:1935112/1‑91 (MQ=255)
                    tATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGtata                                                                       <  2:410340/91‑1 (MQ=255)
                    tATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGtata                                                                       <  2:2768155/91‑1 (MQ=255)
                     atgatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGCTTATAACTCGTAATAATTACGATTTATAAAGTATAg                                                                      >  2:1869512/1‑91 (MQ=255)
                     atgatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAg                                                                      >  1:893736/1‑91 (MQ=255)
                     atgatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAg                                                                      >  2:2926017/1‑91 (MQ=255)
                       gatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGAc                                                                    >  1:931278/1‑91 (MQ=255)
                                 atTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTga                                                          <  2:610037/91‑1 (MQ=255)
                                 atTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTga                                                          <  1:2904340/91‑1 (MQ=255)
                                 atTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTga                                                          >  2:2423105/1‑91 (MQ=255)
                                 atTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTga                                                          <  2:2479909/91‑1 (MQ=255)
                                  ttAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTgag                                                         <  1:919336/91‑1 (MQ=255)
                                  ttAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTgag                                                         <  1:1529983/91‑1 (MQ=255)
                                   tAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATACATCGTAATAAGTACGATTTATAAAGTATAGACCTATCCTTgaga                                                        >  1:2832645/1‑91 (MQ=255)
                                   tAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTgaga                                                        >  1:2724069/1‑91 (MQ=255)
                                    aaCACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGaa                                                       >  1:3814792/1‑91 (MQ=255)
                                    aaCACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGaa                                                       >  2:640811/1‑91 (MQ=255)
                                     acacTCCCTCAATGATCTATTTCAAGTTATTGCTTTTGTTTATAAATCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGaaa                                                      >  2:3877274/1‑91 (MQ=38)
                                       acTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAtt                                                    <  1:1851883/91‑1 (MQ=255)
                                         tCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTc                                                  >  1:3090506/1‑91 (MQ=255)
                                          cccTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCa                                                 >  1:484304/1‑91 (MQ=255)
                                          cccTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCa                                                 >  1:2536769/1‑91 (MQ=255)
                                          cccTCAATGATCTATTTCAAGATATTGTATTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCa                                                 >  1:1693784/1‑91 (MQ=255)
                                           ccTCAATGATCTATTTCAAGTTATTGTTTTTGTTTATTTTTCGTAATAATTACGATTTATAAAGTATAGACCTACCCTACAGCAATTTCaa                                                >  2:589597/1‑91 (MQ=255)
                                              cAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAatat                                             <  1:2401993/91‑1 (MQ=255)
                                              cAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAatat                                             >  1:3000558/1‑91 (MQ=255)
                                              cAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAatat                                             >  1:3228876/1‑91 (MQ=255)
                                                aTGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTa                                           <  1:160172/91‑1 (MQ=255)
                                                  gATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTaaa                                         >  1:2132666/1‑91 (MQ=255)
                                                   aTCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTaaaa                                        >  1:2309573/1‑91 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                       >  2:1116358/1‑91 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                       <  2:1815176/91‑1 (MQ=255)
                                                     cTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGa                                      <  1:364293/91‑1 (MQ=255)
                                                     cTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGa                                      >  1:3017536/1‑91 (MQ=255)
                                                       aTTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaaa                                    <  1:3006057/91‑1 (MQ=255)
                                                        tttCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaaaa                                   >  2:2370651/1‑91 (MQ=255)
                                                         ttCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaaaaa                                  <  1:2313964/91‑1 (MQ=255)
                                                          tCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAAc                                 >  1:984858/1‑91 (MQ=255)
                                                          tCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAAc                                 >  1:3749675/1‑91 (MQ=255)
                                                           cAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACt                                >  1:2363135/1‑91 (MQ=255)
                                                            aaGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACtt                               <  1:1238630/91‑1 (MQ=255)
                                                               ttATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCg                            <  2:2985973/91‑1 (MQ=255)
                                                                      ttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGt                     <  1:762772/91‑1 (MQ=255)
                                                                      ttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGt                     <  1:2351855/91‑1 (MQ=255)
                                                                      ttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGt                     <  1:569708/91‑1 (MQ=255)
                                                                       tttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTa                    >  1:101131/1‑91 (MQ=255)
                                                                           aTTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAgaga                >  1:1639918/1‑91 (MQ=255)
                                                                              ataAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATCTTAAAAGAAAAACTTCCGTTATTGTAGAGAAtt             >  2:2706322/1‑91 (MQ=255)
                                                                              ataAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAAtt             >  2:2909493/1‑91 (MQ=255)
                                                                              ataAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAAtt             >  1:3766878/1‑91 (MQ=255)
                                                                               taAATCGTAATTATTACGATTTATAAAGTATAGACCTATCCTGGATAAATTTTAATATTACAAGACACACTTCCGTTATTGTCGCGGATTc            >  2:3715993/1‑91 (MQ=25)
                                                                                 aaTCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCAt          >  1:3194319/1‑91 (MQ=255)
                                                                                   tCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCt        <  2:928848/91‑1 (MQ=255)
                                                                                   tCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCt        <  2:1140121/91‑1 (MQ=255)
                                                                                   tCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCt        <  1:3305438/91‑1 (MQ=255)
                                                                                   tCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATGCATCt        >  2:3643116/1‑91 (MQ=255)
                                                                                    cGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTa       >  2:766382/1‑91 (MQ=255)
                                                                                    cGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATACTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTa       >  2:220166/1‑91 (MQ=255)
                                                                                    cGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTa       >  1:3728207/1‑91 (MQ=255)
                                                                                    cGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTa       <  2:3263803/91‑1 (MQ=255)
                                                                                    cGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTa       >  1:2634850/1‑91 (MQ=255)
                                                                                     gTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATCAAAAGAAAAACTTCCGTTATTGTAGAGCATTCATCTaa      >  2:1486941/1‑90 (MQ=255)
                                                                                     gTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTAc      <  1:1472711/91‑1 (MQ=255)
                                                                                     gTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTAc      <  2:775578/91‑1 (MQ=255)
                                                                                      taattaTTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAACAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACa     >  2:3535779/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACa     >  2:1065080/1‑91 (MQ=255)
                                                                                        ataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAtt   <  1:2016325/91‑1 (MQ=255)
                                                                                        ataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAtt   <  1:219532/91‑1 (MQ=255)
                                                                                         taatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAttt  <  1:3565909/91‑1 (MQ=255)
                                                                                          |                                                                                         
GGCACCGATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTT  >  USA300TCH1516_ALE/2725262‑2725441

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: