Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 228,312 T→C 8.0% intergenic (+158/+36) ybbH_1 → / ← USA300TCH1516_00200 putative HTH‑type transcriptional regulator YbbH/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE228,3120TC 8.0% 822.2 / 45.4 276intergenic (+158/+36)ybbH_1/USA300TCH1516_00200putative HTH‑type transcriptional regulator YbbH/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (133/121);  new base C (9/13);  total (142/134)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 288 total aligned reads displayed.

TTGCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAA  >  USA300TCH1516_ALE/228222‑228402
                                                                                          |                                                                                          
ttGCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGt                                                                                            >  1:504408/1‑91 (MQ=255)
 tGCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGtt                                                                                           >  1:694781/1‑91 (MQ=255)
 tGCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGtt                                                                                           <  1:3409608/91‑1 (MQ=255)
 tGCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGtt                                                                                           >  2:2159775/1‑91 (MQ=255)
  gCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTc                                                                                          <  1:2809136/91‑1 (MQ=255)
  gCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTc                                                                                          >  1:1710580/1‑91 (MQ=255)
    ttGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAg                                                                                        >  1:3695537/1‑91 (MQ=255)
      gatgaAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGAt                                                                                      <  1:29839/91‑1 (MQ=255)
      gatgaAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGAt                                                                                      <  1:30088/91‑1 (MQ=255)
       atgaAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCaaatc                                                                                     <  2:3295349/91‑5 (MQ=255)
           aaaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAg                                                                                 >  2:2201846/1‑91 (MQ=255)
            aaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgc                                                                                >  1:1826088/1‑91 (MQ=255)
            aaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgc                                                                                >  1:1352173/1‑91 (MQ=255)
             aaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgct                                                                               <  2:698469/91‑1 (MQ=255)
             aaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgct                                                                               >  2:3557347/1‑91 (MQ=255)
               ttgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAacttt                                                                             <  2:2889488/91‑1 (MQ=38)
                tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttg                                                                            <  2:1863146/91‑1 (MQ=255)
                tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttg                                                                            >  2:1390769/1‑91 (MQ=255)
                 gttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgc                                                                           >  2:3357007/1‑91 (MQ=255)
                  ttgCTACGAGTTTGTAAAGCTGGCATAGAGGCTATTTTGAAGATGATAGAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgct                                                                          <  2:2932580/91‑1 (MQ=25)
                  ttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgct                                                                          <  2:3884953/91‑1 (MQ=38)
                   tgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgctt                                                                         >  1:3548650/1‑91 (MQ=255)
                   tgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgctt                                                                         <  1:886213/91‑1 (MQ=255)
                     cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgcttta                                                                       <  2:2750062/91‑2 (MQ=37)
                      tACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCttttt                                                                      <  2:1662090/91‑1 (MQ=255)
                       aCAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgctttaat                                                                     <  1:1178159/91‑4 (MQ=25)
                        cAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgctttaatt                                                                    <  2:1011984/91‑5 (MQ=25)
                          aGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCAt                                                                  <  2:1594218/91‑1 (MQ=255)
                          aGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCAt                                                                  <  1:630507/91‑1 (MQ=255)
                           gTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATa                                                                 >  1:3432651/1‑91 (MQ=255)
                            tttGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAg                                                                <  1:793546/91‑1 (MQ=255)
                             ttGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGa                                                               <  2:1497422/91‑1 (MQ=255)
                              ttacaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGAt                                                              >  2:1819550/5‑91 (MQ=255)
                              tGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGAt                                                              <  1:2983634/91‑1 (MQ=255)
                              tGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGAt                                                              >  1:892055/1‑91 (MQ=255)
                               gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                             >  2:3157173/1‑91 (MQ=255)
                               gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                             >  1:1522353/1‑91 (MQ=255)
                               gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                             >  1:3727520/1‑91 (MQ=255)
                                tAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATgg                                                            >  1:396207/1‑91 (MQ=255)
                                tAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATgg                                                            >  1:951424/1‑91 (MQ=255)
                                 aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                           <  1:2451795/91‑1 (MQ=255)
                                 aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                           >  1:455924/1‑91 (MQ=255)
                                 aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                           <  1:2465009/91‑1 (MQ=255)
                                 aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTAGAGCTTTGCTTTTTTCATAGATGGa                                                           >  1:3621370/1‑91 (MQ=255)
                                     cTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAAtt                                                       <  1:327476/91‑1 (MQ=255)
                                      tGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTa                                                      <  1:3386274/91‑1 (MQ=255)
                                      tGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTa                                                      <  2:3496697/91‑1 (MQ=255)
                                      tGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTa                                                      >  2:573948/1‑91 (MQ=255)
                                       ggCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTaa                                                     <  1:3669299/91‑1 (MQ=255)
                                       ggCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTaa                                                     >  2:2070575/1‑91 (MQ=255)
                                        gCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTaaa                                                    >  2:1236936/1‑91 (MQ=255)
                                          aTAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATg                                                  >  1:146986/1‑91 (MQ=255)
                                              agGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTga                                              <  2:1753981/91‑1 (MQ=255)
                                              agGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTga                                              <  2:1363583/91‑1 (MQ=255)
                                               ggCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaa                                             <  1:2642545/91‑1 (MQ=255)
                                               ggCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaa                                             <  2:843139/91‑1 (MQ=255)
                                                 cTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaaga                                           >  1:2200879/1‑91 (MQ=255)
                                                 cTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaaga                                           >  1:1010672/1‑91 (MQ=255)
                                                  tATTTTGAAAATGATAAAGGCGTGACAAATTTAAGCAAAGTTCAGATTTGAGCTTTGGTTTTTTCATAGATGGAAATTAAATGCTgaagaa                                          <  1:3102920/91‑1 (MQ=255)
                                                  tATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaagaa                                          <  2:1947025/91‑1 (MQ=255)
                                                   aTTTTGAAAATGCTAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                         >  1:2969920/1‑91 (MQ=255)
                                                   aTTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                         >  2:301551/1‑91 (MQ=255)
                                                   aTTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                         >  1:2774412/1‑91 (MQ=255)
                                                    ttttGAAAATGATAAAGGCGTGACAAATTAAAGCACAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaata                                        >  2:641774/1‑91 (MQ=255)
                                                       tGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtt                                     >  2:843559/1‑91 (MQ=255)
                                                        gAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAttt                                    <  2:1295396/91‑1 (MQ=255)
                                                         aaaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtttt                                   >  2:3112109/1‑91 (MQ=255)
                                                         aaaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtttt                                   <  1:1977184/91‑1 (MQ=255)
                                                         aaaaTGATAAAGACGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGTAGAATAAtatt                                   >  2:2523688/1‑91 (MQ=255)
                                                          aaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTa                                  <  2:2619944/91‑1 (MQ=255)
                                                           aaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTaa                                 <  2:170237/91‑1 (MQ=255)
                                                              gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                              <  2:668014/91‑1 (MQ=255)
                                                               aTAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCt                             >  2:1137810/1‑91 (MQ=255)
                                                                 aaaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaa                           >  2:3221857/1‑91 (MQ=255)
                                                                 aaaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaa                           >  1:2347041/1‑91 (MQ=255)
                                                                   aGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAAt                         >  2:2823308/1‑91 (MQ=255)
                                                                    ggCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATg                        >  2:2537372/1‑91 (MQ=255)
                                                                      cGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGaa                      <  2:2324760/91‑1 (MQ=255)
                                                                      cGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTGTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGaa                      <  2:459943/91‑1 (MQ=255)
                                                                       gTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAAc                     <  1:3867034/91‑1 (MQ=255)
                                                                       gTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAAc                     >  1:3033434/1‑91 (MQ=255)
                                                                          aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                  >  2:2998680/1‑91 (MQ=255)
                                                                          aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                  >  1:2949915/1‑91 (MQ=255)
                                                                          aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                  >  2:1584387/1‑91 (MQ=255)
                                                                           cAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAAt                 >  2:1071399/1‑91 (MQ=255)
                                                                                 aaaGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAg           >  1:2578500/1‑91 (MQ=255)
                                                                                  aaGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg          >  2:191347/1‑91 (MQ=255)
                                                                                  aaGCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg          >  1:1599712/1‑91 (MQ=37)
                                                                                   aGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGt         <  1:2514560/91‑1 (MQ=255)
                                                                                    gCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTg        >  1:3339206/1‑91 (MQ=38)
                                                                                    gCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTg        >  2:368552/1‑91 (MQ=38)
                                                                                     cAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc       <  2:454193/91‑1 (MQ=255)
                                                                                     cAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc       >  1:1765671/1‑91 (MQ=255)
                                                                                      aaaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCt      <  1:1324262/91‑1 (MQ=255)
                                                                                      aaaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCt      <  1:1715374/91‑1 (MQ=255)
                                                                                       aaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                         >  2:3635995/1‑55 (MQ=255)
                                                                                       aaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                         <  1:3635995/55‑1 (MQ=255)
                                                                                       aaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTa     <  1:1883646/91‑1 (MQ=255)
                                                                                        aGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTaa    >  2:877846/1‑91 (MQ=255)
                                                                                         gTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTaaa   >  1:1273419/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTGCTTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAA  >  USA300TCH1516_ALE/228222‑228402

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: