Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,243,499 T→A 17.4% intergenic (+77/+32) USA300HOU_RS11475 → / ← pyrG hypothetical protein/CTP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,243,4990TA17.4% 697.7 / 138.1 270intergenic (+77/+32)USA300HOU_RS11475/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base T (119/104);  new base A (20/27);  total (139/131)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 288 total aligned reads displayed.

ACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCG  >  USA300TCH1516_ALE/2243409‑2243589
                                                                                           |                                                                                          
gcTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGttt                                                                                            <  2:1612849/90‑1 (MQ=255)
aCTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGttt                                                                                            <  1:2109430/91‑1 (MQ=255)
 cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTa                                                                                           >  2:1092882/1‑91 (MQ=255)
 cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTa                                                                                           >  2:2752765/1‑91 (MQ=255)
 cTTAGAATTATAAAAGTCACTATTGTTATTTAATACCTTGTTACTATCATATTATTAAAAAAAAACAAAAAACCGGTACTTTCAATGTTTa                                                                                           <  2:773471/91‑1 (MQ=255)
 attAGAATTATAAAAGTCACTATTGATATTTCATACCTTGTCACTAACATAATATTAAAAAAAAAAAAAAAACCGGTACTTTCAATGTTTa                                                                                           <  2:3726680/90‑1 (MQ=255)
   aaGAAATATAAAAGTCACTATTGTTATTTCATAACTTGACAATAACATAATATTAAAAAAAAAAAAAAAACCGGTACTTTCATTGTTAAAc                                                                                         <  2:1160255/90‑1 (MQ=19)
    aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATGATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACa                                                                                        <  2:3055690/91‑1 (MQ=39)
    aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACa                                                                                        >  2:2130422/1‑91 (MQ=255)
     gAATTATACAAGTCCCTATTGTTATTTCATACCTTGTCACTAACATAATATTAGCATTATACAAAAAACCGGTACATTCAATGTTTAACaa                                                                                       >  2:3881277/1‑91 (MQ=255)
     gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                       <  2:461400/91‑2 (MQ=255)
     gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                       <  1:2080799/91‑2 (MQ=255)
      aaTTATAAAAGGCACTATTGTTATTTCATACCTTGTCACTATAATAATATTAAAAAAAAAAAAAAAACCGGTACTTTCAATGTTTAACAAt                                                                                      <  2:665970/91‑1 (MQ=255)
       aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATCTTAAAATTATACGAAAAACCGGTACTTTCGATGTTTAACAATg                                                                                     >  1:3961207/1‑91 (MQ=255)
       aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACattg                                                                                     <  1:1430296/91‑4 (MQ=38)
       aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACattg                                                                                     <  1:1112626/91‑4 (MQ=38)
       aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                     >  1:1419806/1‑91 (MQ=255)
       aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTGCATTGTTAAACattg                                                                                     <  1:3369305/91‑4 (MQ=38)
        ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                    <  2:3494552/91‑1 (MQ=38)
        ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGa                                                                                    >  1:1050373/1‑91 (MQ=255)
        ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGa                                                                                    >  2:1559704/1‑91 (MQ=255)
         tataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaa                                                                                   >  1:608139/1‑91 (MQ=255)
          ataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                  <  1:3536254/91‑1 (MQ=255)
           taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 >  1:2526441/1‑91 (MQ=255)
           taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 <  2:1998016/91‑1 (MQ=255)
             aaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                               >  1:1304883/1‑91 (MQ=255)
              aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                              <  1:1813072/91‑1 (MQ=255)
              aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                              <  1:3241410/91‑1 (MQ=255)
              aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                              >  1:3154342/1‑91 (MQ=255)
               aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             >  2:605897/1‑91 (MQ=255)
               aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             >  1:3845981/1‑91 (MQ=255)
                gTCACTATTGTTATTTCATACCCTTGTCACTAACATAATATCAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             <  2:222771/91‑1 (MQ=39)
                  cACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGt                                                                          <  1:2016468/91‑1 (MQ=255)
                       ttGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttat                                                                     >  1:425283/1‑91 (MQ=255)
                           tATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGc                                                                 <  2:3346001/91‑1 (MQ=255)
                              ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCACTGTTTAACAATGAATGTACCGGTTTTATTAGCaaa                                                              >  1:440578/1‑91 (MQ=255)
                              ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                              >  2:254640/1‑91 (MQ=255)
                              ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                              >  1:1365104/1‑91 (MQ=255)
                               tCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAat                                                             >  1:2180662/1‑91 (MQ=255)
                               tCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAat                                                             >  1:2405178/1‑91 (MQ=255)
                                cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                            >  2:167703/1‑91 (MQ=255)
                                cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                            >  1:203275/1‑91 (MQ=255)
                                  tACCTTGTCACTAACATAATATTAATATTATACGAAAACCCGGTACTTTCAATTTTTAACAATAAAAGTACGGGTTTTATTAGCAAattta                                                          <  2:1451879/91‑1 (MQ=255)
                                  tACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAattta                                                          >  1:1280919/1‑91 (MQ=255)
                                   aCCTTGTCACTAACATCAAATTAAAATTATACGAAAAACCGGTACTATAAATGTTTAACAATGAAAGAACTGGTTATATTAGCAAatgtat                                                         >  2:3914020/1‑91 (MQ=255)
                                   aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                         >  2:1742303/1‑91 (MQ=255)
                                       ttGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattta                                                      >  1:3058733/1‑91 (MQ=255)
                                         gTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAACGGTTTTATTAGCAAatttatttata                                                    >  2:26380/1‑90 (MQ=255)
                                             cTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTGACAATGAAAGTACCGGTTCAATTAGCAAATTTATTTATTttgt                                                >  2:3671059/1‑91 (MQ=39)
                                             cTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgt                                                >  1:2794821/1‑91 (MQ=255)
                                              tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgtt                                               >  2:2428160/1‑91 (MQ=255)
                                              tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgtt                                               <  1:2549236/91‑1 (MQ=255)
                                              tAACATAATATTAAAATTATACGAAAAACCGGTACTTCCAATGTTTAACAATGAAACTACCGGTTTTATTAGCCATTTTCCTTATTttgct                                               >  2:3437929/1‑89 (MQ=255)
                                                gtataattTTAATATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             >  1:3869440/3‑91 (MQ=38)
                                                aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             <  2:2362178/91‑1 (MQ=255)
                                                 cATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGat                                            <  1:1717498/91‑1 (MQ=255)
                                                   taataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGatat                                          <  2:3744855/91‑1 (MQ=255)
                                                    aataTTAAAATTATACGTAAAACCGGTACTTTGAATGTTTAACAATGAGAGTACCGGTTTTATTAGCATATTTATTTATTTTGTTGATAtt                                         <  2:1355595/91‑1 (MQ=38)
                                                     atatTAAAATTATATGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                        >  2:2805771/1‑91 (MQ=255)
                                                       atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                      <  1:2339792/91‑1 (MQ=255)
                                                        ttAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTaa                                     <  1:1022880/91‑1 (MQ=255)
                                                         tAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAAt                                    >  1:3773708/1‑91 (MQ=255)
                                                           aaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTaacga                                  <  1:2498525/91‑4 (MQ=255)
                                                           aaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGa                                  >  1:2076055/1‑91 (MQ=255)
                                                            aaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGaa                                 >  2:1859852/1‑91 (MQ=255)
                                                             aTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAg                                <  1:3734815/91‑1 (MQ=255)
                                                              ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                               <  1:2206878/91‑1 (MQ=255)
                                                               tataCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCt                              <  2:1222659/91‑1 (MQ=255)
                                                                ataCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCtt                             <  1:3216274/91‑1 (MQ=255)
                                                                ataCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCtt                             >  2:828765/1‑91 (MQ=255)
                                                                  aCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                           >  2:3609022/1‑91 (MQ=255)
                                                                  aCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATCAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                           >  2:3469186/1‑91 (MQ=255)
                                                                   cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                          <  2:2533266/91‑1 (MQ=255)
                                                                    gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                         <  1:632653/91‑1 (MQ=255)
                                                                       aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTGTTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                      >  2:3572113/1‑91 (MQ=255)
                                                                       aaaCCGGTACTTTCAATGTTTAACAAGGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                      <  1:3707378/91‑1 (MQ=255)
                                                                         aCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATg                    >  1:2607571/1‑91 (MQ=255)
                                                                         aCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATg                    >  1:2957660/1‑91 (MQ=255)
                                                                           cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                  <  1:3495221/91‑1 (MQ=255)
                                                                            ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAtt                 <  2:2031104/91‑1 (MQ=255)
                                                                             gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTCTTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  2:3089262/1‑91 (MQ=255)
                                                                             gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  1:1205521/1‑91 (MQ=255)
                                                                             gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                <  1:3324166/91‑1 (MQ=255)
                                                                             gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  2:538999/1‑91 (MQ=255)
                                                                               aCTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaa              >  2:3260782/1‑91 (MQ=255)
                                                                                cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa             >  1:2562312/1‑91 (MQ=255)
                                                                                    cattGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc         >  2:1006877/4‑91 (MQ=38)
                                                                                    cattGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc         >  2:2974432/4‑91 (MQ=38)
                                                                                     attGTTAAACATTGAAAGTACCGGTTTTATTAGCACATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg        >  2:1068138/3‑91 (MQ=38)
                                                                                     attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg        >  2:1291091/3‑91 (MQ=38)
                                                                                     attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg        >  1:1875481/3‑91 (MQ=255)
                                                                                     attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAACCg        >  2:893273/3‑91 (MQ=38)
                                                                                      ttGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg       >  1:2280614/2‑91 (MQ=255)
                                                                                      ttGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg       >  1:1606800/2‑91 (MQ=255)
                                                                                      aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAGATGATTTAAAAATCgg       >  2:3036939/1‑91 (MQ=255)
                                                                                      aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg       >  1:3556605/1‑91 (MQ=255)
                                                                                         tttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTg    >  1:36381/1‑91 (MQ=255)
                                                                                          ttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGc   >  1:861150/1‑91 (MQ=255)
                                                                                           tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg  <  1:2669666/91‑1 (MQ=255)
                                                                                           tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg  <  2:1923011/91‑1 (MQ=255)
                                                                                           |                                                                                          
ACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCG  >  USA300TCH1516_ALE/2243409‑2243589

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