Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,350 A→T 12.6% intergenic (+317/+31) 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,3500AT12.6% 716.5 / 76.5 237intergenic (+317/+31)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (100/107);  new base T (7/23);  total (107/130)
Fisher's exact test for biased strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 255 total aligned reads displayed.

ACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTA  >  USA300TCH1516_ALE/168260‑168440
                                                                                          |                                                                                          
aCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGaaaa                                                                                            <  1:270333/91‑1 (MQ=255)
  tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                          <  1:2606811/91‑1 (MQ=255)
  tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                          >  2:146603/1‑91 (MQ=255)
   ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGat                                                                                         >  1:644858/1‑91 (MQ=255)
    tttAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatc                                                                                        >  2:768083/1‑91 (MQ=255)
      tAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATAGTTAAACCTCGAGCTGAAAAGatcat                                                                                      >  1:3619000/1‑91 (MQ=255)
       aaGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatcatc                                                                                     >  2:2581850/1‑91 (MQ=255)
        aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                    >  2:1397559/1‑91 (MQ=255)
          gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                  >  2:1705298/1‑91 (MQ=255)
          gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                  <  1:3696987/91‑1 (MQ=255)
          gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGatgatcttt                                                                                  <  1:1257351/91‑4 (MQ=38)
           aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                 >  2:3011718/1‑91 (MQ=255)
           aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                 <  1:1509025/91‑1 (MQ=255)
            ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTc                                                                                >  2:1480903/1‑91 (MQ=255)
            ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                <  1:1668340/91‑1 (MQ=38)
            ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                <  1:1873767/91‑1 (MQ=37)
            ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                <  2:1019492/91‑1 (MQ=37)
            ttATGTTAAGGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTc                                                                                <  1:3721691/91‑1 (MQ=255)
              aTGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACTAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAg                                                                              >  1:1561057/1‑91 (MQ=255)
              aTGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAg                                                                              >  2:268696/1‑91 (MQ=255)
                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)
                   aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                         >  2:2810523/1‑91 (MQ=255)
                     tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                       <  2:2021986/91‑1 (MQ=255)
                      gAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGt                                                                      >  2:2886482/1‑91 (MQ=255)
                       aaaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtt                                                                     >  2:3132258/1‑91 (MQ=255)
                        aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                    >  2:3334393/1‑91 (MQ=255)
                         aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                   >  2:3845083/1‑91 (MQ=255)
                         aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                   >  1:1043633/1‑91 (MQ=255)
                         aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                   <  2:3534314/91‑1 (MQ=255)
                         aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                   >  1:1991919/1‑91 (MQ=255)
                          agagagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                  <  2:3270034/89‑1 (MQ=255)
                          aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                  <  2:664026/91‑1 (MQ=255)
                          aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                  <  2:364183/91‑1 (MQ=255)
                          aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                  >  2:845214/1‑91 (MQ=255)
                           tAGAGAATATAATATAACGAATCTATATAAAGAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttttt                                                                 <  2:1005416/91‑1 (MQ=255)
                            agagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTtttta                                                                >  1:3264109/1‑91 (MQ=255)
                               gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                             >  1:1424088/1‑91 (MQ=255)
                               gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                             <  1:3476990/91‑1 (MQ=255)
                               gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                             >  2:2860066/1‑91 (MQ=255)
                                aatataatatAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                            >  1:941707/1‑91 (MQ=255)
                                 atataatataACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                           <  1:1634384/91‑1 (MQ=255)
                                  tataatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTtttattttat                                                          <  2:2264384/91‑1 (MQ=255)
                                   ataatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATc                                                         >  1:1489390/1‑91 (MQ=255)
                                     aatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCct                                                       >  1:1744867/1‑91 (MQ=255)
                                       tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctct                                                     <  1:1616304/91‑1 (MQ=255)
                                        ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                                    >  1:3535213/1‑91 (MQ=255)
                                          aaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTg                                                  <  1:945610/91‑1 (MQ=255)
                                            cGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTa                                                <  2:272710/91‑1 (MQ=255)
                                             gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTAt                                               <  2:2036821/91‑1 (MQ=255)
                                               aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                             <  1:1879915/91‑1 (MQ=255)
                                                tCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATggg                                            >  1:293887/1‑91 (MQ=255)
                                                  tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                          <  1:838027/91‑1 (MQ=255)
                                                  tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                          <  1:1958354/91‑1 (MQ=255)
                                                    tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                        >  1:1351091/1‑91 (MQ=255)
                                                         ataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCtttt                                   <  2:2471270/91‑1 (MQ=255)
                                                          taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                                  >  1:3596815/1‑91 (MQ=255)
                                                           aaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                                 <  2:2340691/91‑1 (MQ=255)
                                                            aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                <  2:2240696/91‑1 (MQ=255)
                                                              gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                              >  2:477290/1‑91 (MQ=255)
                                                              gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                              >  2:1136470/1‑91 (MQ=255)
                                                              gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                              >  1:2073560/1‑91 (MQ=255)
                                                              gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                              <  1:3790298/91‑1 (MQ=255)
                                                              gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                              >  2:3362425/1‑91 (MQ=255)
                                                               cACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                             <  2:1088794/91‑1 (MQ=255)
                                                               cACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                             <  2:1089005/91‑1 (MQ=255)
                                                                aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                            >  2:1371604/1‑91 (MQ=255)
                                                                aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                            <  1:2145229/91‑1 (MQ=255)
                                                                  cAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAt                          >  1:1758269/1‑91 (MQ=255)
                                                                   aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                         <  2:1902397/91‑1 (MQ=255)
                                                                   aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                         <  1:3874688/91‑1 (MQ=255)
                                                                   aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                         >  2:3732170/1‑91 (MQ=255)
                                                                   aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                         <  1:1239339/91‑1 (MQ=255)
                                                                   aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                         >  1:2346987/1‑91 (MQ=255)
                                                                      atTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                  <  1:1082109/47‑1 (MQ=255)
                                                                      atTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGt                      <  1:523069/91‑1 (MQ=255)
                                                                      atTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGt                      <  1:41219/91‑1 (MQ=255)
                                                                        ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                    <  2:941707/91‑1 (MQ=255)
                                                                        ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                    <  1:1589464/91‑1 (MQ=255)
                                                                        aaaaaCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                    >  2:1536106/3‑91 (MQ=21)
                                                                         aaaaCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                   >  1:294553/2‑91 (MQ=21)
                                                                          aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc                  <  1:81015/91‑1 (MQ=255)
                                                                           aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                 >  2:733234/1‑91 (MQ=255)
                                                                           aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                 >  1:562911/1‑91 (MQ=255)
                                                                           aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCCTTTACTTGAAATTTGTTTTcc                 >  1:324210/1‑91 (MQ=255)
                                                                            aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                    <  1:73463/39‑1 (MQ=255)
                                                                            aCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                >  1:1955187/1‑91 (MQ=21)
                                                                              cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt              <  1:2312174/91‑1 (MQ=255)
                                                                               tCGTGCTGATATGTTCATTTTTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGTGTTTTCTTTTACTTGAAATTTGTTTTCCTGtt             <  2:2015140/91‑1 (MQ=255)
                                                                                cgaacTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt            >  1:381483/5‑91 (MQ=25)
                                                                                cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt            <  2:515402/91‑1 (MQ=255)
                                                                                      gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt      <  1:2512659/91‑1 (MQ=255)
                                                                                      gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt      >  1:2197221/1‑91 (MQ=255)
                                                                                       aaaaGATCATCATTTCAGTTCGAGGTTTTTTCTTTTATCCTCTTTGTATGGGTATTCGTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt     >  1:3441771/1‑91 (MQ=255)
                                                                                        aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg    >  2:2311010/1‑91 (MQ=255)
                                                                                        aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg    >  2:1787055/1‑91 (MQ=255)
                                                                                          aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTa  <  1:2431547/91‑1 (MQ=255)
                                                                                          aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTCTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTa  >  2:3190592/1‑91 (MQ=255)
                                                                                          |                                                                                          
ACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTA  >  USA300TCH1516_ALE/168260‑168440

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