Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,243,504 A→T 11.3% intergenic (+82/+27) 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,5040AT11.3% 791.7 / 71.5 259intergenic (+82/+27)USA300HOU_RS11475/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base A (122/107);  new base T (21/8);  total (144/115)
Fisher's exact test for biased strand distribution p-value = 7.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 267 total aligned reads displayed.

GAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGA  >  USA300TCH1516_ALE/2243414‑2243594
                                                                                           |                                                                                          
gAATTATACAAGTCCCTATTGTTATTTCATACCTTGTCACTAACATAATATTAGCATTATACAAAAAACCGGTACATTCAATGTTTAACaa                                                                                            >  2:3881277/1‑91 (MQ=255)
gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTCACATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                            >  2:159767/1‑91 (MQ=255)
 aaTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTACAAATATACGAAAAACCGGTACTTTCAATGTTTAACAAt                                                                                           >  2:3668464/1‑91 (MQ=255)
 aaTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAAt                                                                                           >  2:1723443/1‑91 (MQ=255)
  aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATCTTAAAATTATACGAAAAACCGGTACTTTCGATGTTTAACAATg                                                                                          >  1:3961207/1‑91 (MQ=255)
  aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                          >  1:1419806/1‑91 (MQ=255)
  aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                          >  1:3348063/1‑91 (MQ=255)
   ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                         <  1:1162901/91‑1 (MQ=38)
   ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                         <  1:138152/91‑1 (MQ=38)
   ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                         <  2:3494552/91‑1 (MQ=38)
     ataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                       >  1:1269539/1‑91 (MQ=255)
     ataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                       <  1:3536243/91‑1 (MQ=255)
      taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGACAg                                                                                      >  2:2542205/1‑91 (MQ=255)
      taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                      >  1:2526441/1‑91 (MQ=255)
      taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                      <  2:1998016/91‑1 (MQ=255)
        aaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                                    <  1:1938651/91‑1 (MQ=255)
         aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                                   <  1:3241410/91‑1 (MQ=255)
         aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                                   >  1:3154342/1‑91 (MQ=255)
         aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                                   <  1:1813072/91‑1 (MQ=255)
         aaGTCACTAGTATTATTTCATACGTTGTCACTAACATAATATTAAAATTATACGAAAAAGCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                                   <  2:3018242/91‑1 (MQ=255)
          aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAGCCGGTACTTTCATTGTTAAACATTGAAAGTAcc                                                                                  <  1:108736/91‑1 (MQ=37)
          aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                                  >  2:1185808/1‑91 (MQ=255)
           gTCACTATTGTTATTTCATACCCTTGTCACTAACATAATATCAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                                  <  2:222771/91‑1 (MQ=39)
            tCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAGAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCgg                                                                                <  2:1310491/91‑1 (MQ=255)
             cACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGt                                                                               <  1:2016635/91‑1 (MQ=255)
             cACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGt                                                                               <  1:3309216/91‑1 (MQ=255)
             cACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGt                                                                               <  1:3024082/91‑1 (MQ=255)
                  ttGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttat                                                                          >  1:425283/1‑91 (MQ=255)
                   tGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttatt                                                                         >  1:3835431/1‑91 (MQ=255)
                      tATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGc                                                                      <  1:618279/91‑1 (MQ=255)
                       aTTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCa                                                                     <  2:2469975/91‑1 (MQ=255)
                         ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                                   >  2:1033698/1‑91 (MQ=255)
                         ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                                   <  1:3692497/91‑1 (MQ=255)
                          tCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAat                                                                  >  1:2180662/1‑91 (MQ=255)
                           cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                                 <  2:3403532/91‑1 (MQ=255)
                           cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                                 >  2:680248/1‑91 (MQ=255)
                              aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                              >  2:7945/1‑91 (MQ=255)
                                 cTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattt                                                            >  1:2571985/1‑91 (MQ=255)
                                 cTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattt                                                            <  2:1476537/91‑1 (MQ=255)
                                    gTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttatt                                                         <  1:3940765/91‑1 (MQ=255)
                                     tCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttattt                                                        >  2:1825136/1‑91 (MQ=255)
                                      cACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttattct                                                       >  1:2423739/1‑89 (MQ=255)
                                      cACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTAtttt                                                       >  2:1000281/1‑91 (MQ=255)
                                       aCTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttg                                                      <  1:2076993/91‑1 (MQ=255)
                                       aCTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttg                                                      <  1:2500497/91‑1 (MQ=255)
                                        cTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTGACAATGAAAGTACCGGTTCAATTAGCAAATTTATTTATTttgt                                                     >  2:3671059/1‑91 (MQ=39)
                                        cTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgt                                                     >  1:2794821/1‑91 (MQ=255)
                                         tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgtt                                                    <  1:2549236/91‑1 (MQ=255)
                                            cATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGat                                                 <  1:1717498/91‑1 (MQ=255)
                                             ataataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTTCCGGTTTTATTAGCAAATTTATTTATTTTGTTGata                                                <  2:2256267/91‑1 (MQ=255)
                                             ataataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGata                                                <  2:3602375/91‑1 (MQ=255)
                                               aataTTAAAATTATACGTAAAACCGGTACTTTGAATGTTTAACAATGAGAGTACCGGTTTTATTAGCATATTTATTTATTTTGTTGATAtt                                              <  2:1355595/91‑1 (MQ=38)
                                               aataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtt                                              <  2:1228251/91‑1 (MQ=255)
                                                atatTAAAATTATATGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                             >  2:2805771/1‑91 (MQ=255)
                                                atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACACTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                             >  2:3048133/1‑91 (MQ=255)
                                                  atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                           <  1:772452/91‑1 (MQ=255)
                                                  atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                           <  1:2339792/91‑1 (MQ=255)
                                                   ttAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTaa                                          <  1:1022880/91‑1 (MQ=255)
                                                   ttAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTaa                                          <  2:163939/91‑1 (MQ=255)
                                                    tAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAAt                                         >  1:3773708/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)
                                                         ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                                    <  1:463719/91‑1 (MQ=255)
                                                           ataCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCtt                                  <  1:3216237/91‑1 (MQ=255)
                                                             aCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                                >  1:3617257/1‑91 (MQ=255)
                                                             aCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATCAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                                >  2:3469186/1‑91 (MQ=255)
                                                              cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                               >  1:3070715/1‑91 (MQ=255)
                                                                  aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                           >  1:2908102/1‑91 (MQ=255)
                                                                    aCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATg                         >  1:2607571/1‑91 (MQ=255)
                                                                     ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa                        >  2:2098731/1‑91 (MQ=255)
                                                                      cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                       >  1:131337/1‑91 (MQ=255)
                                                                      cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                       <  1:3495221/91‑1 (MQ=255)
                                                                       ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAAtaaat                      <  2:3735992/91‑5 (MQ=255)
                                                                       ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAtt                      <  2:2031104/91‑1 (MQ=255)
                                                                        gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTCTTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     >  2:3089262/1‑91 (MQ=255)
                                                                        gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     >  1:2118370/1‑91 (MQ=255)
                                                                        gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     >  2:2549915/1‑91 (MQ=255)
                                                                        gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     >  2:538999/1‑91 (MQ=255)
                                                                        gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                     <  2:1111795/91‑1 (MQ=255)
                                                                          aCTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaa                   <  2:3745020/91‑1 (MQ=255)
                                                                          aCTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaa                   >  2:3260782/1‑91 (MQ=255)
                                                                          aCTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaa                   >  2:618488/1‑91 (MQ=255)
                                                                           cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa                  >  1:2324801/1‑91 (MQ=255)
                                                                              tcattGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt               >  2:4105/5‑91 (MQ=38)
                                                                               cattGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc              >  2:3518260/4‑91 (MQ=38)
                                                                               cAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc              >  1:1282635/1‑91 (MQ=255)
                                                                                attGTTAAACATTGAAAGTACCGGTTTTATTAGCACATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg             >  2:1068138/3‑91 (MQ=38)
                                                                                attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg             >  1:2058285/3‑91 (MQ=38)
                                                                                attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg             >  1:1875232/3‑91 (MQ=38)
                                                                                attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg             >  2:575465/3‑91 (MQ=38)
                                                                                attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAACCg             >  2:893273/3‑91 (MQ=38)
                                                                                 ttGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTAAAAAATCgg            >  1:3552039/2‑91 (MQ=38)
                                                                                   gTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGt          <  2:2698082/91‑1 (MQ=255)
                                                                                     ttAACAATGAAAGTACCTGTTTTATTAGCAAATTTATGTATTTTGTTGATAGTTGAATGAAGCTTCAATAAGTGATTTAAGAATCGGGTGc        <  2:2396255/91‑1 (MQ=255)
                                                                                     ttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGc        >  1:3861744/1‑91 (MQ=255)
                                                                                      tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg       <  1:2669666/91‑1 (MQ=255)
                                                                                      tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg       <  2:1923011/91‑1 (MQ=255)
                                                                                       aaCAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCgg      <  2:3362461/91‑1 (MQ=255)
                                                                                           aTGAAAGTACCGGTTTTATTAGCAAATTGATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGa  <  2:836734/91‑1 (MQ=255)
                                                                                           |                                                                                          
GAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGA  >  USA300TCH1516_ALE/2243414‑2243594

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: