Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,354 C→G 12.5% intergenic (+321/+27) 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,3540CG12.5% 755.5 / 74.8 240intergenic (+321/+27)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (90/119);  new base G (5/25);  total (95/145)
Fisher's exact test for biased strand distribution p-value = 5.22e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

Only 100 of 261 total aligned reads displayed.

TTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATT  >  USA300TCH1516_ALE/168264‑168444
                                                                                          |                                                                                          
  tAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatcat                                                                                          <  1:1461929/91‑1 (MQ=255)
   aaGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatcatc                                                                                         <  2:48077/91‑1 (MQ=255)
   aaGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatcatc                                                                                         <  2:1896660/91‑1 (MQ=255)
    aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGatgatct                                                                                        <  1:3085270/91‑2 (MQ=39)
     ggATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATcatcat                                                                                       >  2:1764573/1‑91 (MQ=255)
     ggATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATcatcat                                                                                       >  2:1000384/1‑91 (MQ=255)
      gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                      <  2:523659/91‑1 (MQ=255)
       gttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                     <  2:832650/90‑1 (MQ=38)
       aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAATGATGATCtttt                                                                                     <  1:2767718/91‑1 (MQ=39)
       aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                     >  2:2548391/1‑91 (MQ=255)
       aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                     >  1:3225347/1‑91 (MQ=255)
       aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                     <  2:2950725/91‑1 (MQ=38)
       aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                     <  1:3924399/91‑1 (MQ=38)
       aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                     <  1:1104578/91‑1 (MQ=38)
        ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                    <  1:1669398/91‑1 (MQ=37)
        ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                    <  2:3389713/91‑1 (MQ=38)
         tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTcc                                                                                   >  1:64289/1‑90 (MQ=255)
         tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATAGTTAAACCTCGAACTGAAATGATGATCTTTTCa                                                                                   <  2:3885298/91‑1 (MQ=37)
            gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGtt                                                                                <  1:2866895/91‑1 (MQ=255)
            gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAgct                                                                                <  2:7859/91‑3 (MQ=25)
              tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                              <  1:2101279/91‑1 (MQ=255)
              tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                              <  1:3395390/91‑1 (MQ=255)
              tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                              <  2:2984955/91‑1 (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)
                aTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                                            <  2:3791566/91‑1 (MQ=255)
                 tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                           <  1:2297687/91‑1 (MQ=255)
                 tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                           <  1:2079866/91‑1 (MQ=255)
                  gAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGt                                                                          <  1:3807141/91‑1 (MQ=255)
                     aaTAGAGAATATAATTTAACGAATCTATATAAATATGGCACCACTATTTAATCCTCGAGCTGTAACGATAATCATTTCAGTTCGAGGtttt                                                                       <  2:2301674/91‑1 (MQ=255)
                     aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                       >  2:1117518/1‑91 (MQ=255)
                      aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCTAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                      <  2:1356598/91‑1 (MQ=255)
                      aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                      >  2:2882725/1‑91 (MQ=255)
                        agagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTtttta                                                                    >  2:3834606/1‑91 (MQ=255)
                          agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                  <  2:3870215/91‑1 (MQ=255)
                          agaATATAATATAACGAATCTAGATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                  <  2:26540/91‑1 (MQ=255)
                           gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                 >  1:1249138/1‑91 (MQ=255)
                           gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                 >  1:2471406/1‑91 (MQ=255)
                            aatataatatAACGAATCTATATAAATAAGGCACCAATAGTTAAACCGCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                                <  2:3555737/91‑1 (MQ=255)
                               ataatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATc                                                             <  1:1017083/91‑1 (MQ=255)
                                taatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATcc                                                            <  2:3110401/91‑1 (MQ=255)
                                  atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                          >  2:2831470/1‑91 (MQ=255)
                                  atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                          <  1:1424596/91‑1 (MQ=255)
                                  atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                          <  2:755528/91‑1 (MQ=255)
                                   tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctct                                                         <  2:2014042/91‑1 (MQ=255)
                                    ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                                        <  2:2141612/91‑1 (MQ=255)
                                    ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                                        <  2:2619372/91‑1 (MQ=255)
                                           aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                                 >  2:3529377/1‑91 (MQ=255)
                                              tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                              <  1:334524/91‑1 (MQ=255)
                                              tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                              >  1:1321788/1‑91 (MQ=255)
                                              tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                              <  2:1675834/91‑1 (MQ=255)
                                               atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                                             >  2:3241283/1‑91 (MQ=255)
                                               atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                                             >  2:3318838/1‑91 (MQ=255)
                                               atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGCCTTTTATTTTATCCTCTTTGTATGGGTAt                                             >  2:2981524/1‑91 (MQ=255)
                                                tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                            <  2:2886497/91‑1 (MQ=255)
                                                   aaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCtt                                         <  1:277423/91‑1 (MQ=255)
                                                    aataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCttt                                        <  2:3195949/91‑1 (MQ=255)
                                                    aataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCttt                                        >  2:663964/1‑91 (MQ=255)
                                                      taaGGCATCAATTTTTAAACCTCGAGCTGAAATGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                                      <  2:612883/91‑1 (MQ=255)
                                                      taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                                      <  2:2439299/91‑1 (MQ=255)
                                                       aaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                                     <  1:1149385/91‑1 (MQ=255)
                                                       aaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                                     <  2:3261240/91‑1 (MQ=255)
                                                        aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                    >  1:1116520/1‑91 (MQ=255)
                                                         ggCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACtt                                   <  1:1638190/91‑1 (MQ=255)
                                                          gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                                  <  1:192999/91‑1 (MQ=255)
                                                             ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                               >  2:3693261/1‑91 (MQ=255)
                                                             ccAATATTGAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                               <  2:891857/91‑1 (MQ=255)
                                                              cAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGGAAt                              >  2:3729888/1‑91 (MQ=255)
                                                               aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGGGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                             <  2:1317985/91‑1 (MQ=255)
                                                               aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                             >  1:148141/1‑91 (MQ=255)
                                                                   tttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtt                         >  2:117401/1‑91 (MQ=255)
                                                                    ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                        <  1:3733890/91‑1 (MQ=255)
                                                                    ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCCTTGTATGGGTATTCTTTTACTTGAAATTTGttt                        >  2:3914010/1‑91 (MQ=255)
                                                                     tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                       <  1:1862918/91‑1 (MQ=255)
                                                                     tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                       >  1:1855712/1‑91 (MQ=255)
                                                                       aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                     <  1:23613/91‑1 (MQ=255)
                                                                         ccTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTg                   >  2:1632286/1‑91 (MQ=255)
                                                                          cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt                  <  2:2211869/91‑1 (MQ=255)
                                                                          cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt                  <  1:3474751/91‑1 (MQ=255)
                                                                            cGATCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTCTGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt                >  2:2521343/1‑91 (MQ=255)
                                                                            cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt                <  1:2182157/91‑1 (MQ=255)
                                                                              aGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttttt              >  1:2741293/1‑91 (MQ=255)
                                                                                cTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttttttt            >  2:2278883/1‑91 (MQ=255)
                                                                                  gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt          <  2:2760910/91‑1 (MQ=255)
                                                                                  gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt          <  1:455705/91‑1 (MQ=255)
                                                                                   aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt         >  2:1525669/1‑91 (MQ=255)
                                                                                   aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt         >  1:2960371/1‑91 (MQ=255)
                                                                                   aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt         <  2:2358048/91‑1 (MQ=255)
                                                                                   aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt         <  2:3493451/91‑1 (MQ=255)
                                                                                     aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt       <  1:988181/91‑1 (MQ=255)
                                                                                     aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt       <  1:2567578/91‑1 (MQ=255)
                                                                                     aaGATCATCATTTCAGTTCGAGGTTGTTTGTTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt       <  1:3689538/91‑1 (MQ=255)
                                                                                      aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTCTTGtt      >  2:1450697/1‑90 (MQ=255)
                                                                                      aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTa      <  1:715543/91‑1 (MQ=255)
                                                                                      aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTa      >  1:770934/1‑91 (MQ=255)
                                                                                      aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTCTa      >  2:1842023/1‑91 (MQ=255)
                                                                                       gATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaa     <  2:2570222/91‑1 (MQ=255)
                                                                                        atcatcATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaaa    >  1:3513571/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATT  >  USA300TCH1516_ALE/168264‑168444

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: