Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,358 A→T 7.7% intergenic (+325/+23) yfkN_1 → / ← USA300TCH1516_00148 Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE168,3580AT 7.7% 740.4 / 39.2 229intergenic (+325/+23)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (91/120);  new base T (5/13);  total (96/133)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 265 total aligned reads displayed.

AGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTA  >  USA300TCH1516_ALE/168268‑168448
                                                                                          |                                                                                          
 ggATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATcatcat                                                                                           >  2:1000384/1‑91 (MQ=255)
  gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                          >  2:3250611/1‑91 (MQ=255)
  gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                          <  2:523659/91‑1 (MQ=255)
   gttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                         <  2:832650/90‑1 (MQ=38)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                         <  2:1622923/91‑1 (MQ=255)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                         >  1:3225347/1‑91 (MQ=255)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                         <  1:1104578/91‑1 (MQ=38)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                         <  2:2198305/91‑1 (MQ=38)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                         <  2:2950725/91‑1 (MQ=38)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                         <  1:3924399/91‑1 (MQ=38)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                         <  1:2760895/91‑1 (MQ=38)
    ttATGTTAATGCAAATAGAGAATGTAATATAACGAATCTATATAAATAAGGCACCAATAGTTAAGCCTCGAACTGAAATGATGATCTTTTc                                                                                        <  2:3019663/91‑1 (MQ=37)
    ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTc                                                                                        <  1:1998545/91‑1 (MQ=255)
    ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                        <  1:1669497/91‑1 (MQ=37)
        gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAgct                                                                                    <  2:7859/91‑3 (MQ=25)
         ttAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTc                                                                                   <  1:61843/91‑1 (MQ=255)
          tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                                  >  1:3401134/1‑91 (MQ=255)
          tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAGCTCg                                                                                  <  1:2909092/91‑1 (MQ=25)
           aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                                 >  1:1453829/1‑91 (MQ=255)
           aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                                 <  2:2844815/91‑1 (MQ=255)
            aTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                                                <  1:564731/91‑1 (MQ=255)
             tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                               <  1:2079866/91‑1 (MQ=255)
             tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                               <  1:2627034/91‑1 (MQ=255)
              gAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGt                                                                              <  1:3807141/91‑1 (MQ=255)
                 aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                           >  2:1117518/1‑91 (MQ=255)
                 aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGttta                                                                           <  1:2772660/91‑2 (MQ=21)
                  aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCTAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                          <  2:1356598/91‑1 (MQ=255)
                  aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                          >  2:2882725/1‑91 (MQ=255)
                  aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                          >  2:228145/1‑91 (MQ=255)
                    agagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTtttta                                                                        >  2:3834606/1‑91 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                     >  1:1249166/1‑91 (MQ=255)
                       gaATATAATATAACGAAGCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttagtg                                                                     <  2:2620885/91‑4 (MQ=255)
                         atataatataACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAGTTCAGTTCGAGGTTttttatttta                                                                   >  2:3677292/1‑91 (MQ=255)
                           ataatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATc                                                                 <  1:2396489/91‑1 (MQ=255)
                           ataatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATc                                                                 <  1:1017083/91‑1 (MQ=255)
                            taatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATcc                                                                <  2:3110401/91‑1 (MQ=255)
                              atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                              <  1:3731678/91‑1 (MQ=255)
                              atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                              >  2:2831470/1‑91 (MQ=255)
                              atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                              <  1:282/91‑1 (MQ=255)
                              atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                              <  2:755528/91‑1 (MQ=255)
                                ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                                            <  1:2940196/91‑1 (MQ=255)
                                 taaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCttt                                                           <  2:2588696/91‑1 (MQ=255)
                                      aaTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATg                                                      >  1:679454/1‑91 (MQ=255)
                                       aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                                     >  2:3529229/1‑91 (MQ=255)
                                       aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                                     >  2:2773925/1‑91 (MQ=255)
                                          tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                                  >  1:1321788/1‑91 (MQ=255)
                                          tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                                  <  1:334524/91‑1 (MQ=255)
                                           atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGCCTTTTATTTTATCCTCTTTGTATGGGTAt                                                 >  2:2981524/1‑91 (MQ=255)
                                              taaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTActct                                              >  1:3058627/1‑91 (MQ=255)
                                                aataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCttt                                            >  2:663964/1‑91 (MQ=255)
                                                aataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCttt                                            >  2:2229976/1‑91 (MQ=255)
                                                aataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCttt                                            <  1:2356/91‑1 (MQ=255)
                                                  taaGGCATCAATTTTTAAACCTCGAGCTGAAATGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                                          <  2:612883/91‑1 (MQ=255)
                                                  taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                                          <  2:2439299/91‑1 (MQ=255)
                                                    aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTTCt                                        >  2:1749989/1‑91 (MQ=255)
                                                    aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                        >  1:1116520/1‑91 (MQ=255)
                                                    aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                        <  2:3719180/91‑1 (MQ=255)
                                                    aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAATTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                        <  2:3550940/91‑1 (MQ=255)
                                                     ggCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACtt                                       <  1:2146810/91‑1 (MQ=255)
                                                      gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                                      >  2:2698133/1‑91 (MQ=255)
                                                        aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                                    >  1:416795/1‑91 (MQ=255)
                                                         ccACTATTTAAACCTCGAGCTGCAAAGATCATCATTTCAGTTCGAGGTTATTTATATTATCCTATTTGAATGGGTATTCTTTTACTTGaaa                                   >  2:3873411/1‑91 (MQ=255)
                                                         ccAATATTGAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                                   <  2:891857/91‑1 (MQ=255)
                                                           aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                 <  2:2057585/91‑1 (MQ=255)
                                                             tatTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTg                               >  1:1196288/1‑91 (MQ=255)
                                                                 tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                           >  1:164985/1‑91 (MQ=255)
                                                                   aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                         <  2:166496/91‑1 (MQ=255)
                                                                   aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                         >  1:1910271/1‑91 (MQ=255)
                                                                      cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt                      <  1:1467430/91‑1 (MQ=255)
                                                                      cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt                      <  1:3474751/91‑1 (MQ=255)
                                                                       tCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtt                     >  2:1316662/1‑91 (MQ=255)
                                                                        cGATCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTCTGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt                    >  2:2521343/1‑91 (MQ=255)
                                                                        cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt                    <  1:2182157/91‑1 (MQ=255)
                                                                         gaacTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATCTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtctt                   >  2:2349909/4‑91 (MQ=25)
                                                                         gaacTGAAATGATGATCTTTTCAGCTCGAGGTATTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttt                   >  2:3890723/4‑91 (MQ=25)
                                                                            cTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttttttt                >  2:326845/1‑91 (MQ=255)
                                                                             tGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTa               >  2:256654/1‑91 (MQ=37)
                                                                             tGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTa               <  2:3411542/91‑1 (MQ=255)
                                                                              gAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt              >  1:1484291/1‑91 (MQ=38)
                                                                              gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGTGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt              <  2:405324/91‑1 (MQ=255)
                                                                               aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTCTTCCTGTTTTTTTAtt             >  2:3777803/1‑91 (MQ=255)
                                                                                aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAGTTGTTTTCCTGTTTTTTTATTg            <  2:288680/91‑1 (MQ=255)
                                                                                 aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt           <  1:846162/91‑1 (MQ=255)
                                                                                 aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt           <  2:3225298/91‑1 (MQ=255)
                                                                                 aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt           >  1:1857679/1‑91 (MQ=255)
                                                                                 aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt           <  2:1481714/91‑1 (MQ=255)
                                                                                 aaGATCATCATTTCAGTTCGAGGTTGTTTGTTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt           <  1:3689538/91‑1 (MQ=255)
                                                                                  aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTCTa          >  2:1842023/1‑91 (MQ=255)
                                                                                   gATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaa         <  2:1844975/91‑1 (MQ=255)
                                                                                   gATCATCATTTCAGTTCGAGGTTTTGTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaa         <  2:1167965/91‑1 (MQ=255)
                                                                                    atcatcATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaaa        >  1:3513571/1‑91 (MQ=255)
                                                                                     tcatcaTTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAt       >  2:3334517/1‑91 (MQ=255)
                                                                                       atcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt     >  1:249073/1‑91 (MQ=255)
                                                                                       atcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt     >  2:1351297/1‑91 (MQ=255)
                                                                                        tcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtttt    <  1:2560863/91‑1 (MQ=255)
                                                                                         catTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttttt   <  1:3250611/91‑1 (MQ=255)
                                                                                         catTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttttt   <  1:1115441/91‑1 (MQ=255)
                                                                                          atTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTtta  <  1:2190004/91‑1 (MQ=255)
                                                                                          atTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTtta  <  1:952727/91‑1 (MQ=255)
                                                                                          atTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTtta  >  1:1433328/1‑91 (MQ=255)
                                                                                          |                                                                                          
AGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTA  >  USA300TCH1516_ALE/168268‑168448

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: