Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,725,352 A→T 11.3% 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,3520AT11.3% 305.0 / 26.8 97intergenic (‑58/+276)USA300HOU_RS14015/USA300HOU_RS14020Baeyer‑Villiger flavin‑containing monooxygenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (57/29);  new base T (3/8);  total (60/37)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 108 total aligned reads displayed.

GGCACCGATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTTC  >  USA300TCH1516_ALE/2725262‑2725442
                                                                                          |                                                                                          
ggcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaata                                                                                            >  1:768852/1‑91 (MQ=255)
 gcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaataa                                                                                           >  1:1818718/1‑91 (MQ=255)
 gcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaataa                                                                                           >  1:3908537/1‑91 (MQ=255)
 gcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaataa                                                                                           <  1:3147057/91‑1 (MQ=255)
  caccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTaataat                                                                                          >  1:1830458/1‑91 (MQ=255)
  caccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTaataat                                                                                          <  1:3098845/91‑1 (MQ=255)
    ccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTa                                                                                        >  1:524399/1‑91 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                       <  1:884926/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                       <  1:4080354/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                       <  1:3158692/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTAc                                                                                       >  1:2312826/1‑91 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTAc                                                                                       >  1:2372409/1‑91 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTAc                                                                                       >  1:3069559/1‑91 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTAc                                                                                       >  1:3065793/1‑91 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTAc                                                                                       >  1:3065805/1‑91 (MQ=255)
       ataataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGa                                                                                     <  1:2156025/91‑1 (MQ=255)
         aataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGAtt                                                                                   <  1:1129591/91‑1 (MQ=255)
          ataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGAttt                                                                                  <  1:4090536/91‑1 (MQ=255)
          ataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGAttt                                                                                  <  1:4090484/91‑1 (MQ=255)
          ataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAttt                                                                                  >  1:3777459/1‑91 (MQ=255)
            aGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGATTtat                                                                                <  1:2202510/91‑1 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                              >  1:3364300/1‑91 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                              >  1:3368299/1‑91 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                              <  1:2268356/91‑1 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                              >  1:3554170/1‑91 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                              >  1:130197/1‑91 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                              <  1:792285/91‑1 (MQ=255)
              ccACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaa                                                                              <  1:279725/91‑1 (MQ=255)
                  tttATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGta                                                                          >  1:2117105/1‑91 (MQ=255)
                     atgatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAg                                                                       >  1:1534811/1‑91 (MQ=255)
                      tgatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGa                                                                      <  1:1676468/91‑1 (MQ=255)
                      tgatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGa                                                                      <  1:2852768/91‑1 (MQ=255)
                          gTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTa                                                                  <  1:1739092/91‑1 (MQ=255)
                           ttGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTAt                                                                 >  1:2773117/1‑91 (MQ=255)
                           ttGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTAt                                                                 >  1:661695/1‑91 (MQ=255)
                               atatTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCtt                                                             <  1:173097/91‑1 (MQ=255)
                                tatTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTg                                                            <  1:1050438/91‑1 (MQ=255)
                                 atTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCATTga                                                           >  1:3458885/1‑91 (MQ=255)
                                  ttAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTgag                                                          <  1:118625/91‑1 (MQ=255)
                                  ttAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTgag                                                          <  1:2493828/91‑1 (MQ=255)
                                    aaCACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGaa                                                        >  1:3141720/1‑91 (MQ=255)
                                    aaCACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGaa                                                        >  1:294357/1‑91 (MQ=255)
                                     acacTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGaaa                                                       >  1:1565239/1‑91 (MQ=255)
                                      cacTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAt                                                      >  1:2812646/1‑91 (MQ=255)
                                      cacTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAt                                                      >  1:878311/1‑91 (MQ=255)
                                      cacTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAt                                                      >  1:2813488/1‑91 (MQ=255)
                                       acTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAtt                                                     >  1:549871/1‑91 (MQ=255)
                                       acTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAtt                                                     <  1:1808447/91‑1 (MQ=255)
                                        cTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCCTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAttt                                                    >  1:619608/1‑91 (MQ=255)
                                         tCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTc                                                   <  1:1074077/91‑1 (MQ=255)
                                         tCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTc                                                   <  1:212444/91‑1 (MQ=255)
                                         tCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTc                                                   <  1:2841993/91‑1 (MQ=255)
                                            cTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAat                                                >  1:2555585/1‑91 (MQ=255)
                                             tCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAata                                               <  1:812864/91‑1 (MQ=255)
                                              cAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAatat                                              >  1:2391334/1‑91 (MQ=255)
                                               aaTGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATAtt                                             >  1:4115326/1‑91 (MQ=255)
                                                  gATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTaaa                                          >  1:4035971/1‑91 (MQ=255)
                                                  gATCTATTTCAACTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATGaaa                                          >  1:3209593/1‑91 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                        >  1:2286464/1‑91 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                        >  1:2286259/1‑91 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                        >  1:2263498/1‑91 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                        >  1:2263445/1‑91 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                        >  1:2262785/1‑91 (MQ=255)
                                                        tttCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaaaa                                    >  1:1599716/1‑91 (MQ=255)
                                                        tttCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaaaa                                    >  1:3537276/1‑91 (MQ=255)
                                                        tttCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaaaa                                    >  1:3251196/1‑91 (MQ=255)
                                                           cAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACt                                 <  1:1333642/91‑1 (MQ=255)
                                                             aGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTc                               >  1:1880949/1‑91 (MQ=255)
                                                             aGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTc                               >  1:110942/1‑91 (MQ=255)
                                                               ttATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCg                             <  1:1151807/91‑1 (MQ=255)
                                                               ttATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCg                             >  1:123219/1‑91 (MQ=255)
                                                                tATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGt                            <  1:78995/91‑1 (MQ=255)
                                                                   tGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTAt                         >  1:1992971/1‑91 (MQ=255)
                                                                    gTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTAtt                        >  1:2749329/1‑91 (MQ=255)
                                                                    gTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTAtt                        >  1:1644814/1‑91 (MQ=255)
                                                                    gTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTAtt                        >  1:286116/1‑91 (MQ=255)
                                                                     tttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTg                       >  1:3900170/1‑91 (MQ=255)
                                                                      ttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGt                      >  1:3950008/1‑91 (MQ=255)
                                                                      ttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGt                      <  1:2635107/91‑1 (MQ=255)
                                                                       tttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTa                     >  1:1657053/1‑91 (MQ=255)
                                                                       tttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTa                     >  1:3974464/1‑91 (MQ=255)
                                                                       tttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTa                     >  1:3644264/1‑91 (MQ=255)
                                                                        ttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTag                    <  1:2455257/91‑1 (MQ=255)
                                                                        ttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTag                    >  1:312135/1‑91 (MQ=255)
                                                                        ttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTag                    <  1:698533/91‑1 (MQ=255)
                                                                         tGATTATAAATCGTAATAATTACGATTTATAAAGTATAGATCTATCTTTGAGAAATTTAAATATTAAAAGAAAAATTTGCGTTATTGTaga                   >  1:3810162/1‑91 (MQ=38)
                                                                         tGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTaga                   <  1:2542027/91‑1 (MQ=255)
                                                                          gATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTagag                  <  1:454388/91‑1 (MQ=255)
                                                                           tttATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAgaga                 <  1:2860196/90‑1 (MQ=255)
                                                                           aTTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAgaga                 <  1:2860181/91‑1 (MQ=255)
                                                                           aTTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAgaga                 <  1:1523697/91‑1 (MQ=255)
                                                                             tataAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAAt               >  1:2579439/1‑91 (MQ=255)
                                                                                 aaTCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCAt           <  1:2863597/91‑1 (MQ=255)
                                                                                    cGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTa        >  1:2281927/1‑91 (MQ=255)
                                                                                     gTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTAc       <  1:1286544/91‑1 (MQ=255)
                                                                                     gTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTAc       >  1:3514850/1‑91 (MQ=255)
                                                                                       aataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAt     >  1:4147679/1‑91 (MQ=255)
                                                                                       aataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAt     >  1:1356509/1‑91 (MQ=255)
                                                                                          aatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTTc  <  1:2180803/91‑1 (MQ=255)
                                                                                          aatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTTc  >  1:782702/1‑91 (MQ=255)
                                                                                          |                                                                                          
GGCACCGATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTTC  >  USA300TCH1516_ALE/2725262‑2725442

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