Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,432,243 A→T 8.9% intergenic (+590/+57) birA_2 → / ← USA300HOU_RS12525 Bifunctional ligase/repressor BirA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,432,2430AT 8.9% 657.0 / 38.6 199intergenic (+590/+57)birA_2/USA300HOU_RS12525Bifunctional ligase/repressor BirA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (88/93);  new base T (5/13);  total (93/106)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

GTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAAA  >  USA300TCH1516_ALE/2432153‑2432333
                                                                                          |                                                                                          
gTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAaataa                                                                                            <  1:1683953/91‑1 (MQ=255)
   aaTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                         <  1:2607100/91‑1 (MQ=255)
   aaTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                         <  1:475902/91‑1 (MQ=255)
   aaTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                         >  1:1969634/1‑91 (MQ=255)
    aTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                        <  1:2659424/91‑1 (MQ=255)
     tCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                       <  1:435218/91‑1 (MQ=255)
     tCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                       >  1:165032/1‑91 (MQ=255)
      ccTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCACAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATa                                                                                      <  1:3088009/91‑1 (MQ=255)
        tAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGa                                                                                    >  1:107643/1‑91 (MQ=255)
        tAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGa                                                                                    >  2:3518365/1‑91 (MQ=255)
        tAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGa                                                                                    >  2:3038485/1‑91 (MQ=255)
         aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGaa                                                                                   <  1:2857653/91‑1 (MQ=37)
         aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                   >  1:2180554/1‑91 (MQ=255)
         aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                   >  2:770791/1‑91 (MQ=255)
          aTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGaaa                                                                                  <  2:1363696/91‑1 (MQ=25)
           tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                 <  1:1809088/91‑1 (MQ=37)
           tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                 <  1:513666/91‑1 (MQ=25)
           tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                 <  1:561678/91‑1 (MQ=25)
           tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                                 <  1:278695/91‑1 (MQ=255)
           tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                                 <  1:597543/91‑1 (MQ=255)
             cTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCa                                                                               <  1:1883575/91‑1 (MQ=25)
             cTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCa                                                                               <  1:3213355/91‑1 (MQ=25)
             cTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                               >  2:1884613/1‑91 (MQ=255)
               gATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                             >  1:1867126/1‑91 (MQ=255)
                aTGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCt                                                                            >  1:2374691/1‑91 (MQ=255)
                aTGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCt                                                                            <  1:1959510/91‑1 (MQ=255)
                aTGCGCTAATTGAATCTTGATTTTTTCAGTTTCATATTATAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCAGCt                                                                            <  2:1423709/91‑1 (MQ=18)
                 tGCGCTAATTGAATCTTTATTTATTCAGTTCCATAATAAAATGCACATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                           >  2:339835/1‑91 (MQ=255)
                  gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                          <  1:3038596/91‑1 (MQ=255)
                  gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                          <  2:3260602/91‑1 (MQ=255)
                  gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                          <  2:2904302/91‑1 (MQ=255)
                  gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                          <  2:2354753/91‑1 (MQ=255)
                     cTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                       >  1:1222033/1‑91 (MQ=255)
                       aaTTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCa                                                                     <  1:458152/91‑1 (MQ=255)
                       aaTTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCa                                                                     >  1:2251267/1‑91 (MQ=255)
                         ttGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCAtt                                                                   >  1:276066/1‑91 (MQ=255)
                           gAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGc                                                                 <  2:753210/91‑1 (MQ=255)
                           gAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGc                                                                 >  2:2594561/1‑91 (MQ=255)
                            aaTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGcc                                                                >  2:896527/1‑91 (MQ=255)
                             aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                               >  1:2109141/1‑91 (MQ=255)
                             aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                               >  1:3537279/1‑91 (MQ=255)
                             aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                               <  1:789552/91‑1 (MQ=255)
                              tCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTa                                                              <  1:1395120/91‑1 (MQ=255)
                                 ttatttatTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTGTCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaa                                                           <  1:611043/91‑1 (MQ=255)
                                 ttatttatTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaa                                                           >  2:1949581/1‑91 (MQ=255)
                                 ttatttatTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaa                                                           >  1:1324582/1‑91 (MQ=255)
                                 ttatttatTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaa                                                           <  1:2516682/91‑1 (MQ=255)
                                   atttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAc                                                         <  2:2659348/91‑1 (MQ=255)
                                   atttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAc                                                         <  1:978730/91‑1 (MQ=255)
                                    tttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACt                                                        <  2:215197/91‑1 (MQ=255)
                                       attCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTaa                                                     >  1:2673554/1‑91 (MQ=255)
                                        ttCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAAt                                                    >  1:1460255/1‑91 (MQ=255)
                                         tCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATg                                                   >  1:1601539/1‑91 (MQ=255)
                                         tCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATg                                                   >  2:2396301/1‑91 (MQ=255)
                                         tCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATg                                                   >  1:3584610/1‑91 (MQ=255)
                                          cAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAGTTAATGt                                                  >  1:1951409/1‑91 (MQ=255)
                                          cAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGt                                                  <  2:1600076/91‑1 (MQ=255)
                                          cAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGt                                                  >  1:971565/1‑91 (MQ=255)
                                          cAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGt                                                  >  1:1244398/1‑91 (MQ=255)
                                           aGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTc                                                 <  2:1325972/91‑1 (MQ=255)
                                                cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                            <  1:3744345/91‑1 (MQ=255)
                                                cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                            >  1:2196528/1‑91 (MQ=255)
                                                cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                            >  1:97291/1‑91 (MQ=255)
                                                  taataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACct                                          <  1:2948759/91‑1 (MQ=255)
                                                   aataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctc                                         <  2:2828761/91‑1 (MQ=255)
                                                   aataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctc                                         <  2:629162/91‑1 (MQ=255)
                                                   aataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctc                                         <  1:2491720/91‑1 (MQ=255)
                                                     taaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCtt                                       <  1:3123800/91‑1 (MQ=255)
                                                          tGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCt                                  <  1:2606807/91‑1 (MQ=255)
                                                            cATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTAc                                >  2:2445467/1‑91 (MQ=255)
                                                             atatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACg                               >  2:2476963/1‑91 (MQ=255)
                                                                 cAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATa                           <  2:2926904/91‑1 (MQ=255)
                                                                  aaaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATaa                          <  2:1068980/91‑1 (MQ=255)
                                                                   aaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATaaa                         <  2:1801879/91‑1 (MQ=255)
                                                                   aaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATaaa                         <  1:3035090/91‑1 (MQ=255)
                                                                   aaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTACCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATaaa                         <  1:3356272/91‑1 (MQ=255)
                                                                    aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAAt                        <  1:292287/91‑1 (MQ=255)
                                                                     aaaGCTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATa                       >  1:1808353/1‑91 (MQ=21)
                                                                     aaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATa                       <  2:2075492/91‑1 (MQ=255)
                                                                      aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCTTTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc                      >  2:2991773/1‑91 (MQ=255)
                                                                      aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc                      <  2:2767543/91‑1 (MQ=255)
                                                                      aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc                      <  2:1022093/91‑1 (MQ=255)
                                                                       aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACt                     >  1:2511678/1‑91 (MQ=255)
                                                                           gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatt                 >  2:2858716/1‑91 (MQ=255)
                                                                            aTTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatta                >  2:3100402/1‑91 (MQ=25)
                                                                            aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatta                <  1:2439322/91‑1 (MQ=255)
                                                                            aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatta                >  2:3035573/1‑91 (MQ=255)
                                                                             tttCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTAc               <  1:1765216/91‑1 (MQ=255)
                                                                              ttCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACt              <  2:3719665/91‑1 (MQ=255)
                                                                               tCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTg             >  2:1623205/1‑91 (MQ=255)
                                                                                cTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGa            >  2:2302829/1‑91 (MQ=255)
                                                                                 tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat           <  2:1189521/91‑1 (MQ=255)
                                                                                 tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat           <  2:3321288/91‑1 (MQ=255)
                                                                                 tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat           <  1:2604366/91‑1 (MQ=255)
                                                                                 tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTACCTACGCATAAATACTTATTACTGat           >  1:3620873/1‑91 (MQ=255)
                                                                                  aTCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGata          >  1:1110728/1‑91 (MQ=255)
                                                                                   tCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGatat         <  1:3816708/91‑1 (MQ=255)
                                                                                     aaaTAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATaa       >  1:3795709/1‑91 (MQ=255)
                                                                                        taattaatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACga    <  2:2627429/91‑1 (MQ=255)
                                                                                         aattaatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaa   >  1:2184408/1‑91 (MQ=255)
                                                                                          |                                                                                          
GTGAATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAAA  >  USA300TCH1516_ALE/2432153‑2432333

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: