Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,725,352 A→T 19.9% 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,3520AT19.9% 495.8 / 119.3 213intergenic (‑58/+276)USA300HOU_RS14015/USA300HOU_RS14020Baeyer‑Villiger flavin‑containing monooxygenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (91/78);  new base T (24/18);  total (116/97)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 238 total aligned reads displayed.

GGCACCGATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTTC  >  USA300TCH1516_ALE/2725262‑2725442
                                                                                          |                                                                                          
ggcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaata                                                                                            >  2:939980/1‑91 (MQ=255)
ggcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaata                                                                                            >  1:181741/1‑91 (MQ=255)
 gcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaataa                                                                                           >  1:1083771/1‑91 (MQ=255)
 gcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaataa                                                                                           >  2:1738010/1‑91 (MQ=255)
 gcaccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGtaataa                                                                                           >  1:3254327/1‑91 (MQ=255)
  caccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTaataat                                                                                          <  1:1001336/91‑1 (MQ=255)
  caccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTaataat                                                                                          <  1:447708/91‑1 (MQ=255)
    ccgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAAttatta                                                                                        <  2:3677313/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                       <  1:3745149/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                       <  2:2571478/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                       <  1:364535/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTAc                                                                                       <  2:3822621/91‑1 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTAc                                                                                       >  2:2881074/1‑91 (MQ=255)
     cgATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTATTTGATTATAAATCGTAATAATTAc                                                                                       >  1:3597568/1‑91 (MQ=255)
      gATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACg                                                                                      <  1:3626334/91‑1 (MQ=255)
      gATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACg                                                                                      <  1:3784186/91‑1 (MQ=255)
       ataataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGa                                                                                     <  1:525316/91‑1 (MQ=255)
       ataataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGa                                                                                     >  2:2514107/1‑91 (MQ=255)
        taataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAt                                                                                    >  2:1337098/1‑91 (MQ=255)
        taataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAt                                                                                    <  2:2229865/91‑1 (MQ=255)
        taataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAt                                                                                    >  1:581824/1‑91 (MQ=255)
        taataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAt                                                                                    >  2:304453/1‑91 (MQ=255)
         aataGCCACTTTATGATGTTGCATATTAACACCCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGAtt                                                                                   <  2:2250508/91‑1 (MQ=255)
          ataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGAttt                                                                                  <  1:2014669/91‑1 (MQ=255)
          ataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAttt                                                                                  >  2:17446/1‑91 (MQ=255)
          ataGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGAttt                                                                                  >  1:3856267/1‑91 (MQ=255)
           taGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTATTACGATTta                                                                                 <  1:1675133/91‑1 (MQ=255)
              ccACTTTATGATGTTGCATCCAAACACGCCCTCAAGGATCTATTGCAGGTGATTGCTTTCGATTATAAATCGTAGTGATCACGAGTTATaa                                                                              >  2:672769/1‑91 (MQ=255)
               cACTTTATGATGTTGCATATTAATTCTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATTAGTACGATGTATaat                                                                             <  2:3527147/91‑2 (MQ=38)
               cACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATaaa                                                                             >  1:1326114/1‑91 (MQ=255)
                aCTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAg                                                                            <  1:2926197/91‑1 (MQ=255)
                 cTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGt                                                                           <  1:3501595/91‑1 (MQ=255)
                    tATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGtata                                                                        <  2:112710/91‑1 (MQ=255)
                    tATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGtata                                                                        <  2:2172003/91‑1 (MQ=255)
                     atgatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGa                                                                                     >  1:734434/1‑77 (MQ=255)
                       gatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGAc                                                                     >  2:3286523/1‑91 (MQ=255)
                       gatgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGAc                                                                     >  1:1983101/1‑91 (MQ=255)
                        atgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGAcc                                                                    >  1:3143666/1‑91 (MQ=255)
                         tgTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCt                                                                   >  1:2264466/1‑91 (MQ=255)
                              cATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCt                                                              <  2:1302569/91‑1 (MQ=255)
                               atatTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCtt                                                             <  1:520629/91‑1 (MQ=255)
                               atatTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCtt                                                             >  1:1384293/1‑91 (MQ=255)
                               atatTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCtt                                                             <  2:1319421/91‑1 (MQ=255)
                               atatTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCtt                                                             <  1:252629/91‑1 (MQ=255)
                                    aaCACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGaa                                                        >  2:3168114/1‑91 (MQ=255)
                                      cacTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAt                                                      >  2:3732117/1‑91 (MQ=255)
                                       acTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAtt                                                     >  2:691893/1‑91 (MQ=255)
                                       acTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAAtt                                                     >  1:2922537/1‑91 (MQ=255)
                                           ccTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCaa                                                 >  2:492720/1‑91 (MQ=255)
                                             tCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAata                                               >  2:3020702/1‑91 (MQ=255)
                                              cAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAatat                                              >  2:2298295/1‑91 (MQ=255)
                                              cAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAatat                                              >  2:1581440/1‑91 (MQ=255)
                                                aTGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTa                                            <  2:594813/91‑1 (MQ=255)
                                                 tGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTaa                                           <  1:935269/91‑1 (MQ=255)
                                                 tGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTaa                                           >  1:2545334/1‑91 (MQ=255)
                                                   aTCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTaaaa                                         >  1:491409/1‑91 (MQ=255)
                                                   aTCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTaaaa                                         <  2:283149/91‑1 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAg                                        <  2:2202897/91‑1 (MQ=255)
                                                    tCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATCTATAAAGTATAGACCTATCCTTGAGACATTTCAATATTAAAAg                                        >  2:2304521/1‑91 (MQ=255)
                                                     cTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGa                                       <  1:2169667/91‑1 (MQ=255)
                                                      tATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaa                                      >  2:2773051/1‑91 (MQ=255)
                                                      tATTGCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaa                                      <  1:1641025/91‑1 (MQ=255)
                                                       aTTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGaaa                                     <  1:2691093/91‑1 (MQ=255)
                                                          tCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAAc                                  <  1:1741582/91‑1 (MQ=255)
                                                           cAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACt                                 <  2:2705519/91‑1 (MQ=255)
                                                             aGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTc                               >  1:3451118/1‑91 (MQ=255)
                                                             aGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTc                               <  2:145915/91‑1 (MQ=255)
                                                             aGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTc                               <  1:798451/91‑1 (MQ=255)
                                                                 aTTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGtt                           >  2:1637642/1‑91 (MQ=255)
                                                                 aTTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGtt                           >  2:194815/1‑91 (MQ=255)
                                                                  ttGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTa                          <  2:3185305/91‑1 (MQ=255)
                                                                    gTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTAtt                        >  1:471781/1‑91 (MQ=255)
                                                                     tttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTg                       >  2:18977/1‑91 (MQ=255)
                                                                      ttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGt                      <  1:1901492/91‑1 (MQ=255)
                                                                      ttttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGt                      <  1:2298295/91‑1 (MQ=255)
                                                                        ttGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTag                    <  2:1613468/91‑1 (MQ=255)
                                                                             tataAATCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAAt               >  2:2575410/1‑91 (MQ=255)
                                                                              ataAATCGTAATTATTACGATTTATAAANTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAAtt              >  1:1456168/1‑91 (MQ=255)
                                                                               taAATCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTc             >  2:2126112/1‑91 (MQ=255)
                                                                               taAATCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATCACAAGAAAAACTTCCGTTATTGTAGAGAATTc             >  2:1372941/1‑91 (MQ=255)
                                                                                aaaTCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCa            >  2:2191263/1‑91 (MQ=255)
                                                                                aaaTCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCa            <  2:18112/91‑1 (MQ=255)
                                                                                aaaTCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCa            >  1:2912880/1‑91 (MQ=255)
                                                                                aaaTCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCa            <  1:633067/91‑1 (MQ=255)
                                                                                 aaTCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTACGTTATTGTAGAGAATTCAt           >  1:870963/1‑91 (MQ=255)
                                                                                 aaTCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCAt           >  2:868475/1‑91 (MQ=255)
                                                                                   tCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCt         >  1:556733/1‑91 (MQ=255)
                                                                                   tCGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCt         >  1:2503088/1‑91 (MQ=255)
                                                                                   tCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCt         <  2:2922656/91‑1 (MQ=255)
                                                                                   tCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCt         <  1:1337098/91‑1 (MQ=255)
                                                                                    cGTAATTATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTa        >  1:3207704/1‑91 (MQ=255)
                                                                                     gTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTAc       <  1:3564271/91‑1 (MQ=255)
                                                                                     gTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTAc       <  1:1199560/91‑1 (MQ=255)
                                                                                      taataaTTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACa      >  2:1017584/1‑91 (MQ=255)
                                                                                       aataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAt     <  1:2854911/91‑1 (MQ=255)
                                                                                       aataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAt     >  1:2779862/1‑91 (MQ=255)
                                                                                       aataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAt     <  1:151265/91‑1 (MQ=255)
                                                                                        ataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTACAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAtt    >  2:1082662/1‑91 (MQ=255)
                                                                                        ataatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAtt    >  1:2281392/1‑91 (MQ=255)
                                                                                         taatTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACAttt   <  2:115085/91‑1 (MQ=255)
                                                                                          |                                                                                          
GGCACCGATAATAGCCACTTTATGATGTTGCATATTAACACTCCCTCAATGATCTATTTCAAGTTATTGTTTTTGATTATAAATCGTAATAATTACGATTTATAAAGTATAGACCTATCCTTGAGAAATTTCAATATTAAAAGAAAAACTTCCGTTATTGTAGAGAATTCATCTACATTTC  >  USA300TCH1516_ALE/2725262‑2725442

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