Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 228,316 G→A 11.3% intergenic (+162/+32) 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,3160GA11.3% 602.7 / 44.8 208intergenic (+162/+32)ybbH_1/USA300TCH1516_00200putative HTH‑type transcriptional regulator YbbH/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (95/89);  new base A (15/9);  total (110/98)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

TTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCAA  >  USA300TCH1516_ALE/228226‑228406
                                                                                          |                                                                                          
ttGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAg                                                                                            <  2:3163818/91‑1 (MQ=255)
     gaAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGCTTTg                                                                                       <  1:3128326/91‑1 (MQ=255)
       aaaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAg                                                                                     >  1:3213845/1‑91 (MQ=255)
       aaaaGTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTAGTTTGAAAATGATAAAGGCGTGACAAATGAAAGCAAAGCTCAAATCTGaa                                                                                     <  2:426571/91‑2 (MQ=255)
       aaaaGTGGTGCGACAGGTTTGTAAAGCTGGCATAGAGGCTAGTTTGAAAATGGTAACGGCGTGACAAGGTAGAGCAAAGCTCAAATCTGaa                                                                                     <  2:3466128/91‑2 (MQ=255)
         aaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgct                                                                                   >  1:871221/1‑91 (MQ=255)
          aTTGTTGGTACAAGTTTGGAGAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctt                                                                                  <  2:1299429/91‑1 (MQ=255)
          aTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctt                                                                                  <  2:2091657/91‑1 (MQ=255)
           ttgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgcttt                                                                                 <  1:1003991/91‑1 (MQ=255)
           ttgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgcttt                                                                                 >  1:2393361/1‑91 (MQ=255)
           ttgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgcttt                                                                                 >  1:2168985/1‑91 (MQ=255)
            tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttg                                                                                <  2:3659894/91‑1 (MQ=255)
            tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttg                                                                                >  1:1411439/1‑91 (MQ=255)
             gttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgc                                                                               >  2:3320724/1‑91 (MQ=255)
              ttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAATGCGTGACAAATTAAATCAAAGTTCAGATTTGAgcttttct                                                                              >  2:820603/1‑91 (MQ=255)
               tgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgctt                                                                             >  1:2429238/1‑91 (MQ=255)
               tgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgctt                                                                             >  1:771472/1‑91 (MQ=255)
                gCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgcttt                                                                            >  1:3810956/1‑91 (MQ=255)
                 cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCtttt                                                                           >  1:1796489/1‑91 (MQ=255)
                 cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCtttt                                                                           <  1:2435958/91‑1 (MQ=255)
                 cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCtttt                                                                           >  2:605790/1‑91 (MQ=255)
                  tACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgctttaa                                                                          <  2:3359784/91‑3 (MQ=37)
                    cAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTc                                                                        <  1:2614108/91‑1 (MQ=255)
                      aGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCAt                                                                      <  1:688080/91‑1 (MQ=255)
                        tttGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAg                                                                    >  2:2986924/1‑91 (MQ=255)
                        tttGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAg                                                                    >  2:2927408/1‑91 (MQ=255)
                          tGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGAt                                                                  >  2:874450/1‑91 (MQ=255)
                           gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                                 >  1:35227/1‑91 (MQ=255)
                           gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                                 >  1:2025697/1‑91 (MQ=255)
                            tAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATgg                                                                <  2:40608/91‑1 (MQ=255)
                            tAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATgg                                                                >  2:2728368/1‑91 (MQ=255)
                             aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                               >  2:2141611/1‑91 (MQ=255)
                             aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                               >  2:1520411/1‑91 (MQ=255)
                               aaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGaaa                                                             <  1:2728560/91‑1 (MQ=255)
                                aCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAAt                                                            <  1:3844544/91‑1 (MQ=255)
                                 cTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAAtt                                                           >  1:3928689/1‑91 (MQ=255)
                                 cTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATGAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAAtt                                                           <  1:3949151/91‑1 (MQ=255)
                                  tGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTa                                                          <  2:1849921/91‑1 (MQ=255)
                                   ggCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTaa                                                         <  1:3433214/91‑1 (MQ=255)
                                    gCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTaaa                                                        >  2:219599/1‑91 (MQ=255)
                                     gaTAGAGGGGATTGTGAGAATGATAAAGGCGTGGCAGATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTGTCATAGATGGAAATTAAAt                                                       <  2:2595971/90‑1 (MQ=255)
                                     cATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAAt                                                       >  1:3051443/1‑91 (MQ=255)
                                       tAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGc                                                     >  1:1880348/1‑91 (MQ=255)
                                        agagGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCt                                                    <  1:3598602/91‑1 (MQ=255)
                                          agGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgg                                                  >  2:1599575/1‑90 (MQ=255)
                                           ggCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaa                                                 <  1:1993107/91‑1 (MQ=255)
                                             cTATTTTGAAGATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTCTTTTCATAGATGGAAATTAAATGCTgaaga                                               >  2:3532347/1‑91 (MQ=255)
                                               aTTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                             >  2:1134185/1‑91 (MQ=255)
                                                ttttGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaata                                            >  1:2384681/1‑91 (MQ=255)
                                                 tttGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaataa                                           >  2:683949/1‑91 (MQ=255)
                                                    gAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAttt                                        >  1:876014/1‑91 (MQ=255)
                                                     aaaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtttt                                       <  2:2634258/91‑1 (MQ=255)
                                                     aaaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtttt                                       >  2:2397267/1‑91 (MQ=255)
                                                      aaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTa                                      <  2:1047170/91‑1 (MQ=255)
                                                       aaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTaa                                     <  2:2306201/91‑1 (MQ=255)
                                                          gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                                  <  2:1823747/91‑1 (MQ=255)
                                                          gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                                  <  2:2054120/91‑1 (MQ=255)
                                                          gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                                  <  2:3502254/91‑1 (MQ=255)
                                                            tAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTa                                >  1:959620/1‑91 (MQ=255)
                                                             aaaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaa                               <  1:1497503/91‑1 (MQ=255)
                                                             aaaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaa                               >  1:552656/1‑91 (MQ=255)
                                                              aaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGCATAATTTTAATGCTaaa                              >  1:3381231/1‑91 (MQ=255)
                                                              aaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaaa                              <  2:1775132/91‑1 (MQ=255)
                                                              aaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaaa                              >  2:2320283/1‑91 (MQ=255)
                                                               aGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTTTTATTTGGTTATTAAATGCTGAAGAATAATTTTAATGCTAAAt                             >  2:2704157/1‑91 (MQ=255)
                                                               aGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAAt                             >  1:1243698/1‑91 (MQ=255)
                                                                 gCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGa                           >  2:2282480/1‑91 (MQ=255)
                                                                 gCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGa                           <  2:3578514/91‑1 (MQ=255)
                                                                  cGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGAGGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGaa                          <  2:3126142/91‑1 (MQ=255)
                                                                   gTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAAc                         <  2:761239/91‑1 (MQ=255)
                                                                   gTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAAc                         >  2:49779/1‑91 (MQ=255)
                                                                      aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                      <  2:1639435/91‑1 (MQ=255)
                                                                       cAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAAt                     <  1:1051019/91‑1 (MQ=255)
                                                                         aaTTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATAc                   >  1:638854/1‑91 (MQ=255)
                                                                          aTTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTACTGCTAAATGAACTAATACg                  >  1:3070597/1‑91 (MQ=255)
                                                                          aTTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACg                  >  2:2136274/1‑91 (MQ=255)
                                                                           ttAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGa                 <  2:2562041/91‑1 (MQ=255)
                                                                           ttAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGa                 <  1:880400/91‑1 (MQ=255)
                                                                            aaaagcaaagcTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGaa                >  1:2495520/2‑91 (MQ=37)
                                                                              aaGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg              <  1:1120378/91‑1 (MQ=255)
                                                                              aaGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg              <  1:1299604/91‑1 (MQ=255)
                                                                              aaGCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg              >  2:1441126/1‑91 (MQ=37)
                                                                               aGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGt             <  1:1799349/91‑1 (MQ=255)
                                                                                gCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTg            >  2:100296/1‑91 (MQ=38)
                                                                                gCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTg            >  1:950310/1‑91 (MQ=38)
                                                                                 cAAAGTACAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc           <  2:1667479/91‑1 (MQ=255)
                                                                                 cAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc           >  1:187074/1‑91 (MQ=38)
                                                                                 cAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc           >  1:1537288/1‑91 (MQ=38)
                                                                                 cAAAGCTCAAATCTGAACTTCGCTTTTTTCATAGATGGAAATAAAATGATGATGAATAATTTTAATGCTAAATGAACTAATACGAAGGAGc           >  2:2526971/1‑91 (MQ=255)
                                                                                  aaaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCt          >  1:933727/1‑91 (MQ=255)
                                                                                  aaaGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCt          >  1:2805770/1‑91 (MQ=38)
                                                                                   aaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTa         >  1:3686696/1‑91 (MQ=255)
                                                                                    aGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTaa        >  1:601223/1‑91 (MQ=255)
                                                                                       tCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAg     <  2:2476475/91‑1 (MQ=255)
                                                                                        ccgATTTAAGCTTTGCTTCTCTCATAGATGGAACTTAAATGCTGAAGAATAATTTTAATGCTGAATGAACTAATACGAAGGTGCTAAAAGc    >  2:496140/3‑91 (MQ=255)
                                                                                        cAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGc    >  2:155524/1‑91 (MQ=255)
                                                                                         aGATTTGAGCTTTGCTTTTTTTTTTTTTTTTTATTAAATGCTGAAGAATAATTTTTATGCTGTTTGAACTAATACGAAGGTGCTAATAGCa   >  2:3785272/1‑91 (MQ=255)
                                                                                          gATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCaa  <  2:3846352/91‑1 (MQ=255)
                                                                                          gATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCaa  >  2:31778/1‑91 (MQ=255)
                                                                                          gATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCaa  >  2:2499498/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTGATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCAA  >  USA300TCH1516_ALE/228226‑228406

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