Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,358 A→T 10.3% 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,3580AT10.3% 744.1 / 58.6 232intergenic (+325/+23)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (98/109);  new base T (7/17);  total (105/127)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 251 total aligned reads displayed.

AGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTA  >  USA300TCH1516_ALE/168268‑168448
                                                                                          |                                                                                          
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                         >  2:3011718/1‑91 (MQ=255)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                         <  2:1095413/91‑1 (MQ=38)
    ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTc                                                                                        >  2:1480903/1‑91 (MQ=255)
    ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                        <  2:1656975/91‑1 (MQ=37)
    ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                        <  1:1873767/91‑1 (MQ=37)
    ttATGTTAAGGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTc                                                                                        <  1:3721691/91‑1 (MQ=255)
     tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCa                                                                                       <  2:284944/91‑1 (MQ=255)
     tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCa                                                                                       <  2:508079/91‑1 (MQ=255)
     tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCa                                                                                       <  2:2927205/91‑1 (MQ=37)
      aTGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAg                                                                                      >  2:268696/1‑91 (MQ=255)
      aTGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAg                                                                                      <  1:1044630/91‑1 (MQ=37)
       tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGt                                                                                     <  1:1686070/91‑1 (MQ=255)
       tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGt                                                                                     <  2:3350247/91‑1 (MQ=255)
        gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGtt                                                                                    <  1:1490186/91‑1 (MQ=255)
        gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGtt                                                                                    <  1:2850737/91‑1 (MQ=255)
         ttAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTc                                                                                   <  2:1730720/91‑1 (MQ=255)
         ttAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAgctc                                                                                   <  1:1522258/91‑4 (MQ=25)
          tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                                  <  1:56403/91‑1 (MQ=255)
            aTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                                                <  1:2026942/91‑1 (MQ=255)
            aTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                                                <  1:2384570/91‑1 (MQ=255)
             tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                               <  1:3276600/91‑1 (MQ=255)
              gAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGt                                                                              >  2:2886482/1‑91 (MQ=255)
                aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                            >  2:3334393/1‑91 (MQ=255)
                 aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                           <  2:3534314/91‑1 (MQ=255)
                 aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGttta                                                                           <  2:3628975/91‑2 (MQ=21)
                  aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                          <  1:702022/91‑1 (MQ=255)
                  aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                          >  2:845214/1‑91 (MQ=255)
                   tAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttttt                                                                         >  1:1287639/1‑91 (MQ=255)
                   tAGAGAATATAATATAACGAATCTATATAAAGAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttttt                                                                         <  2:1005416/91‑1 (MQ=255)
                     gagaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttat                                                                       >  1:272710/1‑91 (MQ=255)
                      agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                      >  1:3525812/1‑91 (MQ=255)
                      agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                      >  1:2081499/1‑91 (MQ=255)
                      agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                      <  2:3757066/91‑1 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                     >  2:751043/1‑91 (MQ=255)
                        aatataatatAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                                    >  1:941707/1‑91 (MQ=255)
                         atataatataACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                                   <  1:1634384/91‑1 (MQ=255)
                            tgtcataaGGAATCTTTATAAAGAAGGCACCAATATTTAAACCTCTAGCTGAACAGATGATCATTTCAGTTGGAGGTTTTTTAGTTTATcc                                                                <  2:1351096/87‑1 (MQ=255)
                             aatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCct                                                               >  1:1744867/1‑91 (MQ=255)
                              atataaCGAATCGATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                              <  2:3223603/91‑1 (MQ=255)
                               tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctct                                                             >  2:686463/1‑91 (MQ=255)
                               tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctct                                                             <  1:1996409/91‑1 (MQ=255)
                                ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                                            <  1:1881156/91‑1 (MQ=255)
                                    cGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTa                                                        <  1:2581850/91‑1 (MQ=255)
                                    cGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTa                                                        <  2:2060376/91‑1 (MQ=255)
                                    cGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTa                                                        >  1:2340691/1‑91 (MQ=255)
                                     gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTAt                                                       <  2:2036627/91‑1 (MQ=255)
                                      aaTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATg                                                      >  1:3648555/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)
                                           atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                                                 >  1:1970733/1‑91 (MQ=255)
                                           atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                                                 <  1:226000/91‑1 (MQ=255)
                                            tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                                >  2:1390431/1‑91 (MQ=255)
                                              taaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGtcttct                                              >  1:812577/1‑91 (MQ=255)
                                                 ataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCtttt                                           <  2:2471270/91‑1 (MQ=255)
                                                   aaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                                         <  2:2340691/91‑1 (MQ=255)
                                                    aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                        <  2:2240816/91‑1 (MQ=255)
                                                     ggCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACtt                                       <  1:1395319/91‑1 (MQ=255)
                                                      gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                                      >  1:1619125/1‑91 (MQ=255)
                                                      gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                                      <  1:3790298/91‑1 (MQ=255)
                                                      gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                                      >  2:2021308/1‑91 (MQ=255)
                                                      gCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTg                                      >  2:1136470/1‑91 (MQ=255)
                                                       cACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                                     <  2:3087359/91‑1 (MQ=255)
                                                       cACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                                     <  2:1088794/91‑1 (MQ=255)
                                                         ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                                   <  1:1480903/91‑1 (MQ=255)
                                                         ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGCTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaca                                   >  2:399283/1‑89 (MQ=255)
                                                          cAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAt                                  <  2:2776028/91‑1 (MQ=255)
                                                           aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                 >  2:3730550/1‑91 (MQ=255)
                                                           aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                 <  2:1902397/91‑1 (MQ=255)
                                                           aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                 <  1:3874723/91‑1 (MQ=255)
                                                           aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                 <  1:1239339/91‑1 (MQ=255)
                                                              atTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGt                              <  1:523069/91‑1 (MQ=255)
                                                                ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                            >  2:1042926/1‑91 (MQ=255)
                                                                ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                            <  2:941707/91‑1 (MQ=255)
                                                                ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                            <  2:2277329/91‑1 (MQ=255)
                                                                 tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                           <  1:1827434/91‑1 (MQ=255)
                                                                 aaaaCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                           >  1:294306/2‑91 (MQ=21)
                                                                  aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc                          <  2:1452253/91‑1 (MQ=255)
                                                                   aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                         >  1:562912/1‑91 (MQ=255)
                                                                   aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCCTTTACTTGAAATTTGTTTTcc                         >  1:324210/1‑91 (MQ=255)
                                                                    aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                            <  1:73463/39‑1 (MQ=255)
                                                                    aCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                        >  1:1955646/1‑91 (MQ=21)
                                                                        cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGttt                    <  1:1567777/91‑1 (MQ=255)
                                                                         gaacTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttt                   >  1:3570715/4‑91 (MQ=25)
                                                                              gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt              <  1:2512659/91‑1 (MQ=255)
                                                                               aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt             >  1:989654/1‑91 (MQ=255)
                                                                               aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt             >  1:914130/1‑91 (MQ=255)
                                                                               aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCGTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt             >  1:3441789/1‑91 (MQ=255)
                                                                                aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg            >  2:1787055/1‑91 (MQ=255)
                                                                                aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg            >  2:2311010/1‑91 (MQ=255)
                                                                                 aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt           <  1:2187716/91‑1 (MQ=255)
                                                                                 aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt           <  1:1605290/91‑1 (MQ=255)
                                                                                  aGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTCTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTa          >  2:3190592/1‑91 (MQ=255)
                                                                                    atcatcATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaaa        <  2:2484059/91‑1 (MQ=255)
                                                                                     tcatcaTTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGGTTTTTTATTGTAAAt       >  2:1980309/1‑91 (MQ=255)
                                                                                      catcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtt      <  2:644858/91‑1 (MQ=255)
                                                                                      catcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGATATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtt      <  1:3087933/91‑1 (MQ=255)
                                                                                       atcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt     >  1:430580/1‑91 (MQ=255)
                                                                                       atcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt     <  1:1042926/91‑1 (MQ=255)
                                                                                          atTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTtta  <  1:2445417/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTA  >  USA300TCH1516_ALE/168268‑168448

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: