Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 228,319 T→C 11.7% intergenic (+165/+29) 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,3190TC11.7% 555.4 / 52.2 207intergenic (+165/+29)ybbH_1/USA300TCH1516_00200putative HTH‑type transcriptional regulator YbbH/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (97/84);  new base C (15/9);  total (113/94)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 212 total aligned reads displayed.

ATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCAA  >  USA300TCH1516_ALE/228229‑228406
                                                                                          |                                                                                       
   aaaaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGa                                                                                      <  2:2623760/91‑1 (MQ=255)
      aaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgct                                                                                   >  1:871221/1‑91 (MQ=255)
      aaTTGTTGCTACAAGCTTGTAAAACTTGCATAGAGGCTATTTTGAAAATTATAAAGGCGTGACAAATTAAAGCAACGTTCAGATTTGAgct                                                                                   >  2:3469563/1‑91 (MQ=255)
       aTTGTTGGTACAAGTTTGGAGAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctt                                                                                  <  2:1299429/91‑1 (MQ=255)
       aTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctt                                                                                  <  2:2091657/91‑1 (MQ=255)
       aTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactt                                                                                  <  1:468104/91‑1 (MQ=38)
       aTTGTTGCTACAAGTTGGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactt                                                                                  <  2:2294880/91‑1 (MQ=38)
        ttgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgcttt                                                                                 <  1:1003991/91‑1 (MQ=255)
        ttgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgcttt                                                                                 >  1:2393361/1‑91 (MQ=255)
         tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttg                                                                                >  1:1411439/1‑91 (MQ=255)
         tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttg                                                                                <  2:3659894/91‑1 (MQ=255)
         tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttg                                                                                <  2:226184/91‑1 (MQ=38)
          gttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgc                                                                               >  2:3320724/1‑91 (MQ=255)
           ttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAATGCGTGACAAATTAAATCAAAGTTCAGATTTGAgcttttct                                                                              >  2:820603/1‑91 (MQ=255)
            tgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgctt                                                                             >  1:771472/1‑91 (MQ=255)
             gCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgcttt                                                                            >  1:3810956/1‑91 (MQ=255)
              cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCtttt                                                                           >  2:605790/1‑91 (MQ=255)
               tACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgctttaa                                                                          <  2:3359784/91‑3 (MQ=37)
                   aGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCAt                                                                      >  1:2089036/1‑91 (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)
                       tGTAAAACTGACATAGAGGCTATTCTGAAAATGTGAAAGGCGTGACAAATTAAAGCAGAGTACAGATATGAGCTTTGCTTTTTTCATAggt                                                                  >  2:3143194/1‑89 (MQ=255)
                        gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                                 >  1:35227/1‑91 (MQ=255)
                        gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                                 >  1:2764234/1‑91 (MQ=255)
                         tAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATgg                                                                >  2:2728368/1‑91 (MQ=255)
                         tAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATgg                                                                <  2:40608/91‑1 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                               >  2:1520411/1‑91 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                               >  2:2141611/1‑91 (MQ=255)
                             aCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAAt                                                            <  2:182255/91‑1 (MQ=255)
                             aCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAAt                                                            <  1:3844544/91‑1 (MQ=255)
                               tGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTa                                                          <  2:1849921/91‑1 (MQ=255)
                                ggCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTaa                                                         <  1:3433117/91‑1 (MQ=255)
                                 gCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTaaa                                                        >  2:219599/1‑91 (MQ=255)
                                  cATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAAt                                                       >  1:3051443/1‑91 (MQ=255)
                                    tAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGc                                                     >  1:3938578/1‑91 (MQ=255)
                                    tAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGc                                                     >  1:1880348/1‑91 (MQ=255)
                                     agagGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCt                                                    <  1:3598602/91‑1 (MQ=255)
                                      gagGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTg                                                   >  1:3475213/1‑91 (MQ=255)
                                       agGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgg                                                  >  2:1599575/1‑90 (MQ=255)
                                        ggCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaa                                                 <  1:1993107/91‑1 (MQ=255)
                                          cTATTTTGAAGATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTCTTTTCATAGATGGAAATTAAATGCTgaaga                                               >  2:3532347/1‑91 (MQ=255)
                                          cTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaaga                                               <  2:3110368/91‑1 (MQ=255)
                                           tATTTTGCAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaagaa                                              <  2:3516496/91‑1 (MQ=255)
                                            aTTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                             >  2:1134185/1‑91 (MQ=255)
                                             ttttGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaata                                            >  1:2384681/1‑91 (MQ=255)
                                             ttttGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaata                                            <  1:2208545/91‑1 (MQ=255)
                                              tttGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaataa                                           >  2:683949/1‑91 (MQ=255)
                                               ttGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaataat                                          >  1:2484987/1‑91 (MQ=255)
                                                tGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtt                                         <  1:2610807/91‑1 (MQ=255)
                                                 gAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAttt                                        >  1:876204/1‑91 (MQ=255)
                                                  aaaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtttt                                       >  2:2397267/1‑91 (MQ=255)
                                                  aaaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtttt                                       <  2:469478/91‑1 (MQ=255)
                                                   aaaTGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTa                                      <  2:1047170/91‑1 (MQ=255)
                                                      tGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATg                                   <  1:2688980/91‑1 (MQ=255)
                                                       gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                                  <  1:1780676/91‑1 (MQ=255)
                                                       gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                                  <  2:2054120/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)
                                                           aaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaaa                              <  2:1775132/91‑1 (MQ=255)
                                                           aaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaaa                              >  2:2587204/1‑91 (MQ=255)
                                                           aaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaaa                              >  2:2320283/1‑91 (MQ=255)
                                                            aGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTTTTATTTGGTTATTAAATGCTGAAGAATAATTTTAATGCTAAAt                             >  2:2704157/1‑91 (MQ=255)
                                                            aGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAAt                             >  1:1243698/1‑91 (MQ=255)
                                                             ggCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATg                            <  1:3224443/91‑1 (MQ=255)
                                                              gCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGa                           <  2:3578514/91‑1 (MQ=255)
                                                                gTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAAc                         >  2:49779/1‑91 (MQ=255)
                                                                gTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAAc                         <  2:761239/91‑1 (MQ=255)
                                                                 tGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACt                        <  2:2556204/91‑1 (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                 <  1:317037/91‑1 (MQ=255)
                                                                        ttAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGa                 <  1:880400/91‑1 (MQ=255)
                                                                           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            >  1:950310/1‑91 (MQ=38)
                                                                             gCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTg            >  2:100296/1‑91 (MQ=38)
                                                                              cAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc           <  1:2433816/91‑1 (MQ=255)
                                                                              cAAAGTACAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc           <  2:1667479/91‑1 (MQ=255)
                                                                              cAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc           >  1:1537288/1‑91 (MQ=38)
                                                                              cAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc           >  1:187074/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:31778/1‑91 (MQ=255)
                                                                                       gATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCaa  <  2:3846352/91‑1 (MQ=255)
                                                                                       gATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCaa  >  2:2499498/1‑91 (MQ=255)
                                                                                          |                                                                                       
ATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCAA  >  USA300TCH1516_ALE/228229‑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: