Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 598,492 A→G 10.8% intergenic (+41/+241) tuf → / ← yxeP_2 Elongation factor Tu/putative hydrolase YxeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE598,4920AG10.8% 655.7 / 57.4 241intergenic (+41/+241)tuf/yxeP_2Elongation factor Tu/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base A (118/96);  new base G (8/18);  total (126/115)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

Only 100 of 247 total aligned reads displayed.

AAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTT  >  USA300TCH1516_ALE/598402‑598582
                                                                                          |                                                                                            
 aGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCaa                                                                                             >  1:1543379/1‑91 (MQ=255)
  ggTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAAt                                                                                            <  1:1507126/91‑1 (MQ=255)
    tGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCg                                                                                          >  2:1740024/1‑91 (MQ=255)
      gACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTAAATCGAg                                                                                        >  1:1962540/1‑91 (MQ=255)
       aCGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgg                                                                                       <  2:2853428/91‑1 (MQ=255)
       aCGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGACCCCTCGATTGAgg                                                                                       <  2:2663437/91‑1 (MQ=255)
        cGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAggg                                                                                      <  1:919608/91‑1 (MQ=255)
         gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                     <  1:1118163/91‑1 (MQ=255)
         gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                     <  1:1280133/91‑1 (MQ=255)
            cTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTATTTCCTTAGATTTTATATAAACAGAAGCTCCCTCAATCGAGGGGTCt                                                                                 >  2:2284724/1‑91 (MQ=255)
            cTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCt                                                                                 <  2:2377422/91‑1 (MQ=255)
             tGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtt                                                                                >  1:1176416/1‑91 (MQ=255)
              gTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGACCCCTCGATTGAGGGATCttc                                                                               <  1:627314/91‑2 (MQ=38)
               tAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtttt                                                                             >  1:3648582/1‑91 (MQ=255)
                  gATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCttttttt                                                                          >  1:2919314/1‑91 (MQ=255)
                  gATCAGGCGTTGTTACTGAAATCATTAAAGAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCttctttt                                                                          <  1:3173383/91‑1 (MQ=25)
                     cAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAt                                                                       >  1:1353313/1‑91 (MQ=255)
                     cAGGCGTTGTTACTGAAATAATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTata                                                                       <  2:1122329/91‑3 (MQ=18)
                       ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                     <  1:2264499/91‑1 (MQ=255)
                        gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                    >  2:165597/1‑91 (MQ=255)
                        gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                    >  1:2785544/1‑91 (MQ=255)
                        gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatata                                                                     <  1:2281300/91‑5 (MQ=18)
                           ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                                 <  2:2685199/91‑1 (MQ=255)
                           ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                                 >  2:1190560/1‑91 (MQ=255)
                            tgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaaa                                                                <  1:511837/91‑1 (MQ=255)
                               tACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAttt                                                             <  2:3405993/91‑1 (MQ=255)
                                aCTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtttt                                                            >  2:3294694/1‑91 (MQ=255)
                                 tagAAATTATTTAATGATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTg                                                           >  2:2222816/3‑91 (MQ=255)
                                 cTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTACATTTTg                                                           >  2:1304006/1‑91 (MQ=255)
                                 cTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTg                                                           >  2:3078569/1‑91 (MQ=255)
                                  tggaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGt                                                          <  2:1283409/88‑1 (MQ=255)
                                  tGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGt                                                          <  1:2978291/91‑1 (MQ=255)
                                    aaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTaa                                                        <  2:2907769/91‑1 (MQ=255)
                                    aaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTaa                                                        >  2:729834/1‑91 (MQ=255)
                                       tCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTACATTTTGTAATgg                                                     >  2:408434/1‑91 (MQ=255)
                                       tCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATgg                                                     >  2:1610309/1‑91 (MQ=255)
                                         aTTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCt                                                   <  2:2499380/91‑1 (MQ=255)
                                          ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTTAATGTGTAAATTTTGTAATGGCt                                                   <  1:1323799/91‑1 (MQ=255)
                                          ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTa                                                  >  1:323258/1‑91 (MQ=255)
                                          ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTa                                                  >  2:3630308/1‑91 (MQ=255)
                                           tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                                 >  1:2454468/1‑91 (MQ=255)
                                           tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                                 <  1:3832777/91‑1 (MQ=255)
                                             aataatTTCTCATTTCTTAGATTTTATATAATAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTc                                               <  2:2605280/91‑1 (MQ=255)
                                             aataatTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTc                                               >  1:3943022/1‑91 (MQ=255)
                                             aataatTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTc                                               >  1:3311209/1‑91 (MQ=255)
                                              ataatTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCg                                              >  1:1694693/1‑91 (MQ=255)
                                               taatttctaatttcTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGa                                             >  2:2724492/1‑91 (MQ=255)
                                               taatttctaatttcTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGa                                             >  1:839042/1‑91 (MQ=255)
                                                  tttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGAttt                                          >  1:2126103/1‑91 (MQ=255)
                                                     ctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTACCTTTTGTAATCGCTATTAGATTTaga                                       >  2:3427207/1‑91 (MQ=255)
                                                      taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                      <  1:989452/91‑1 (MQ=255)
                                                      taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                      >  2:536138/1‑91 (MQ=255)
                                                      taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                      >  1:1648630/1‑91 (MQ=255)
                                                       aatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGTGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaag                                     >  2:1414603/1‑91 (MQ=255)
                                                         tttctTAGATTTTATATAAAAAGAAGGTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTGAATTTTGTAGTGGCTATTCGATTTAgaagaa                                   <  2:1464062/91‑1 (MQ=255)
                                                             ttAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                               >  1:2109839/1‑91 (MQ=255)
                                                             ttAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                               >  1:1463517/1‑91 (MQ=255)
                                                              tAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaata                              <  2:3805799/91‑1 (MQ=255)
                                                              tAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaata                              >  2:646625/1‑91 (MQ=255)
                                                               aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                             >  1:420385/1‑91 (MQ=255)
                                                               aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                             >  2:2636889/1‑91 (MQ=255)
                                                                gcttttATATAAACAGAAGCCCCCTCAACCGAGGGGTCCTTTTTTAATGTGTAAATTTTGTAAGGGCTATTCGATTTAGACGAACaataat                            >  2:3651771/3‑91 (MQ=255)
                                                                     tatataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatg                       <  1:369813/91‑1 (MQ=255)
                                                                     tatataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatg                       <  2:169497/91‑1 (MQ=255)
                                                                      atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                      <  2:2468789/91‑1 (MQ=255)
                                                                      atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                      >  1:1870711/1‑91 (MQ=255)
                                                                      atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                      >  1:2867861/1‑91 (MQ=255)
                                                                       tataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaa                     >  2:2038095/1‑91 (MQ=255)
                                                                        ataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaaa                    >  1:1018340/1‑91 (MQ=255)
                                                                          aaaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAga                  <  2:2484053/91‑1 (MQ=255)
                                                                          aaaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAga                  <  2:1216438/91‑1 (MQ=255)
                                                                           aaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgac                 <  1:2589366/91‑1 (MQ=255)
                                                                           aaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgac                 <  2:1582081/91‑1 (MQ=255)
                                                                            aaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgact                >  2:3502222/1‑91 (MQ=255)
                                                                            aaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgact                <  1:312429/91‑1 (MQ=255)
                                                                              agaagaTCCCTCAATCGAGGGGTCTTTTTTTACTGTGTCCATTTTGTACTGGCTATTCGCTTTCGAAGAACACTAATTGACGAACgactga              >  2:456300/1‑91 (MQ=255)
                                                                              agaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga              <  1:2702459/91‑1 (MQ=255)
                                                                              agaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga              <  1:2198970/91‑1 (MQ=255)
                                                                              agaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga              <  1:241134/91‑1 (MQ=255)
                                                                               gaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac             >  2:743794/1‑91 (MQ=255)
                                                                               gaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac             <  2:1847744/91‑1 (MQ=255)
                                                                                 agaccccTCGATTGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta           >  1:714759/5‑91 (MQ=255)
                                                                                 agaccccTCGATTGAAGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATATAGAAGAACAATAATTGATGAAAGACTGACta           >  1:3366106/5‑91 (MQ=255)
                                                                                 agaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta           >  1:139557/1‑91 (MQ=255)
                                                                                  gaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGACTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa          >  2:3314615/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat         >  1:3524402/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat         >  2:873764/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat         >  1:3734842/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat         <  2:3109029/91‑1 (MQ=255)
                                                                                    tCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata        >  2:3273171/1‑91 (MQ=255)
                                                                                    tCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata        >  1:1704519/1‑91 (MQ=255)
                                                                                    tCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATCCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata        >  2:2895903/1‑91 (MQ=255)
                                                                                     cccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa       >  1:3439256/1‑91 (MQ=255)
                                                                                     cccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa       >  1:773625/1‑91 (MQ=255)
                                                                                     cccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa       >  2:2080624/1‑91 (MQ=255)
                                                                                      ccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa      >  2:3548375/1‑91 (MQ=255)
                                                                                      ccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa      <  2:392818/91‑1 (MQ=255)
                                                                                        tCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAc    <  2:2834162/91‑1 (MQ=255)
                                                                                         cAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACt   <  1:3612792/91‑1 (MQ=255)
                                                                                          aaTCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACtt  <  2:814664/91‑1 (MQ=255)
                                                                                          |                                                                                            
AAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTT  >  USA300TCH1516_ALE/598402‑598582

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: