Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,343 G→A 11.8% intergenic (+310/+38) 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,3430GA11.8% 710.2 / 61.6 245intergenic (+310/+38)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (103/113);  new base A (7/22);  total (110/135)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 247 total aligned reads displayed.

ATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTT  >  USA300TCH1516_ALE/168253‑168430
                                                                                          |                                                                                       
  ttAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCt                                                                                       <  1:1008810/91‑1 (MQ=255)
   tAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTg                                                                                      >  2:1522321/1‑91 (MQ=255)
    aaaaCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGa                                                                                     >  1:1928778/1‑91 (MQ=255)
     aaaCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGaa                                                                                    <  2:3821934/91‑1 (MQ=255)
      aaCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGaaa                                                                                   <  1:2686687/91‑1 (MQ=255)
         tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATGAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                <  1:3259981/91‑1 (MQ=255)
         tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                <  1:2606811/91‑1 (MQ=255)
         tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                <  1:112193/91‑1 (MQ=255)
          ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATCTAACGAATCTATATAAATAAGGCACCACTATTTAAACCTCGAGCTGAGAAGat                                                                               >  1:3604754/1‑91 (MQ=255)
          ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGat                                                                               >  2:2627246/1‑91 (MQ=255)
           tttAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatc                                                                              >  2:2819072/1‑91 (MQ=255)
           tttAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatc                                                                              <  1:3733720/91‑1 (MQ=255)
              aaGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatcatc                                                                           >  2:2581850/1‑91 (MQ=255)
                 gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                        >  2:1705298/1‑91 (MQ=255)
                 gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGatgatcttt                                                                        <  1:1257351/91‑4 (MQ=38)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                       <  1:1509025/91‑1 (MQ=255)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                       <  1:3502724/91‑1 (MQ=255)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                       >  2:1734529/1‑91 (MQ=255)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                       >  2:3011718/1‑91 (MQ=255)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                       <  2:1095413/91‑1 (MQ=38)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                       <  2:845288/91‑1 (MQ=38)
                   ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                      <  1:1873767/91‑1 (MQ=37)
                   ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                      <  2:1926896/91‑1 (MQ=38)
                   ttATGTTAAGGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTc                                                                      <  1:3721691/91‑1 (MQ=255)
                    tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCa                                                                     <  2:508079/91‑1 (MQ=255)
                    tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCa                                                                     <  1:696509/91‑1 (MQ=255)
                    tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCa                                                                     <  2:3857547/91‑1 (MQ=37)
                     aTGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACTAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAg                                                                    >  1:1561057/1‑91 (MQ=255)
                     aTGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAg                                                                    <  1:1044630/91‑1 (MQ=37)
                       gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGtt                                                                  <  1:2850737/91‑1 (MQ=255)
                        ttAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTc                                                                 <  1:202627/91‑1 (MQ=255)
                         tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                <  1:56403/91‑1 (MQ=255)
                         tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCAGAAAAGATCATGATTTCAGTTCg                                                                >  2:2958569/1‑91 (MQ=255)
                         tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCACGAGCTGAAAAGATCATCATTTCAGTTCg                                                                >  1:3110739/1‑91 (MQ=255)
                           aTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                              <  1:2384597/91‑1 (MQ=255)
                           aTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCCCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                              >  1:2509394/1‑91 (MQ=255)
                            tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTGAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                             <  1:3276551/91‑1 (MQ=255)
                               aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                          >  2:3334393/1‑91 (MQ=255)
                               aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGttt                                                          <  1:927436/91‑1 (MQ=21)
                                aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                         <  2:3534314/91‑1 (MQ=255)
                                 agagagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                        <  2:3270034/89‑1 (MQ=255)
                                 aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAg                                                                    <  2:1242611/79‑1 (MQ=255)
                                 aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                        <  2:364183/91‑1 (MQ=255)
                                    gagaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttat                                                     >  1:272710/1‑91 (MQ=255)
                                     agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                    >  1:2081499/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                   >  2:2860066/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                   >  1:1424088/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                   <  2:605988/91‑1 (MQ=255)
                                          ataatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATc                                               >  1:1489390/1‑91 (MQ=255)
                                           taatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATcc                                              <  2:2728671/91‑1 (MQ=255)
                                            aatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCct                                             <  1:415861/91‑1 (MQ=255)
                                              tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctct                                           >  2:686463/1‑91 (MQ=255)
                                              tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctct                                           <  1:1996409/91‑1 (MQ=255)
                                               ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                          >  1:3535213/1‑91 (MQ=255)
                                                 aaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTg                                        <  1:945379/91‑1 (MQ=255)
                                                   cGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTa                                      >  1:2340691/1‑91 (MQ=255)
                                                   cGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTa                                      <  2:2060376/91‑1 (MQ=255)
                                                    gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTAt                                     <  2:2036627/91‑1 (MQ=255)
                                                    gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTAt                                     >  2:2899368/1‑91 (MQ=255)
                                                      aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                   <  1:1879915/91‑1 (MQ=255)
                                                       tCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATggg                                  >  1:293887/1‑91 (MQ=255)
                                                       tCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATggg                                  <  1:3821328/91‑1 (MQ=255)
                                                         tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                <  2:2718622/91‑1 (MQ=255)
                                                          atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                               >  1:3635882/1‑91 (MQ=255)
                                                           tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                              >  1:1351091/1‑91 (MQ=255)
                                                            ataaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTc                             >  1:1645289/1‑91 (MQ=255)
                                                             taaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCt                            >  2:3156291/1‑91 (MQ=255)
                                                                  aaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                       <  2:2340691/91‑1 (MQ=255)
                                                                   aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                      <  2:2240816/91‑1 (MQ=255)
                                                                    ggCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACtt                     >  2:1477288/1‑91 (MQ=255)
                                                                     gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                    >  2:3362425/1‑91 (MQ=255)
                                                                     gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                    >  1:1619125/1‑91 (MQ=255)
                                                                     gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                    >  2:2021308/1‑91 (MQ=255)
                                                                      cACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                   <  2:3087359/91‑1 (MQ=255)
                                                                      cACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                   <  2:1088960/91‑1 (MQ=255)
                                                                       aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                  <  1:2145229/91‑1 (MQ=255)
                                                                       aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                  >  2:1371604/1‑91 (MQ=255)
                                                                        ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                 <  1:1705298/91‑1 (MQ=255)
                                                                         cAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAt                >  1:1758269/1‑91 (MQ=255)
                                                                         cAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTCTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAt                >  2:151692/1‑91 (MQ=255)
                                                                          aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt               >  2:3730550/1‑91 (MQ=255)
                                                                          aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt               <  1:1239339/91‑1 (MQ=255)
                                                                          aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt               <  2:1902397/91‑1 (MQ=255)
                                                                          aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt               <  1:3874688/91‑1 (MQ=255)
                                                                             atTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGt            <  1:523069/91‑1 (MQ=255)
                                                                              tttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtt           >  2:101023/1‑91 (MQ=255)
                                                                               ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt          <  1:1589464/91‑1 (MQ=255)
                                                                               ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt          <  2:2277329/91‑1 (MQ=255)
                                                                               ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt          >  2:1042926/1‑91 (MQ=255)
                                                                                aaaaCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt         >  1:294553/2‑91 (MQ=21)
                                                                                 aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc        <  2:1452253/91‑1 (MQ=255)
                                                                                 aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc        <  1:81015/91‑1 (MQ=255)
                                                                                  aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc       >  1:562912/1‑91 (MQ=255)
                                                                                  aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc       <  1:1970487/91‑1 (MQ=255)
                                                                                  aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc       >  2:237856/1‑91 (MQ=255)
                                                                                   aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                          <  1:73463/39‑1 (MQ=255)
                                                                                   aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                          >  2:73463/1‑39 (MQ=255)
                                                                                   aCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt      >  1:1955646/1‑91 (MQ=21)
                                                                                      tCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtt   >  1:2306385/1‑91 (MQ=255)
                                                                                       cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt  >  2:80959/1‑91 (MQ=255)
                                                                                          |                                                                                       
ATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTT  >  USA300TCH1516_ALE/168253‑168430

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