Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,869,042 G→A 6.5% intergenic (‑449/+21) ytnP ← / ← trmB putative quorum‑quenching lactonase YtnP/tRNA (guanine‑N(7)‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,869,0420GA 6.5% 661.2 / 21.5 200intergenic (‑449/+21)ytnP/trmBputative quorum‑quenching lactonase YtnP/tRNA (guanine‑N(7)‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (102/85);  new base A (9/4);  total (111/89)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 214 total aligned reads displayed.

AGTTAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCA  >  USA300TCH1516_ALE/1868952‑1869132
                                                                                          |                                                                                          
aGTTAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTAtgtg                                                                                            <  1:2800618/91‑1 (MQ=255)
   tAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTAc                                                                                         <  1:504411/91‑1 (MQ=255)
   tAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTAc                                                                                         <  1:505087/91‑1 (MQ=255)
    aaaaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACa                                                                                        <  1:2968738/91‑1 (MQ=255)
    aaaaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACa                                                                                        >  1:3085850/1‑91 (MQ=255)
    aaaaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACa                                                                                        >  1:3404008/1‑91 (MQ=255)
     aaaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACAt                                                                                       >  1:363530/1‑91 (MQ=255)
      aaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATa                                                                                      >  2:190760/1‑91 (MQ=255)
      aaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATa                                                                                      >  2:188403/1‑91 (MQ=255)
        aCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATAcc                                                                                    <  2:2465235/91‑1 (MQ=255)
          tttttCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCt                                                                                  >  1:1463821/1‑91 (MQ=255)
           ttttCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATAcaatc                                                                                 <  1:2953763/91‑5 (MQ=37)
               ccATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATaa                                                                             >  1:2940084/1‑91 (MQ=255)
               ccATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATaa                                                                             >  2:1770234/1‑91 (MQ=255)
                cATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATACAATCATAAg                                                                            <  1:3581673/91‑1 (MQ=25)
                 aTCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGt                                                                           <  2:567786/91‑1 (MQ=255)
                 aTCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATACAATCATAAGt                                                                           <  1:1762090/91‑1 (MQ=21)
                   ctgtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGc                                                                         <  2:1030305‑M2/80‑1 (MQ=255)
                   ctgtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGc                                                                         <  2:3589371‑M2/80‑1 (MQ=255)
                   cAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGGACGAAAAATCACTTATGATTGTATGTGTACATACCCTCATAAGTGc                                                                         <  1:3097286/91‑1 (MQ=255)
                    aGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGct                                                                        <  1:1506447/91‑1 (MQ=255)
                     gtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGctc                                                                       <  2:2426359‑M2/82‑1 (MQ=255)
                      tccgatttATTATAGAACCTTTATCTATAGAACCATTATAGGACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGctct                                                                      <  2:3410390‑M2/83‑1 (MQ=255)
                        cgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTAt                                                                    >  2:2616020‑M2/7‑91 (MQ=255)
                          atttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTAttt                                                                  >  1:3658245‑M2/5‑91 (MQ=255)
                              attatAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATAGTATGTGTACATACCCTCATAACTGCTCTATTATtgt                                                              >  2:1172538/1‑91 (MQ=255)
                               ttatAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTtgtg                                                             >  2:3080935/1‑91 (MQ=255)
                                tataGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGa                                                            <  2:2837147/91‑1 (MQ=255)
                                 ataGAACCTTTATCTGTAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGaa                                                           >  1:233791/1‑91 (MQ=255)
                                 ataGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGaa                                                           >  1:1058642/1‑91 (MQ=255)
                                 ataGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGaa                                                           >  1:1574265/1‑91 (MQ=255)
                                 ataGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGaa                                                           <  1:1242907/91‑1 (MQ=255)
                                   aGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATg                                                         <  1:2137093/91‑1 (MQ=255)
                                     aaCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGa                                                       <  2:248023/91‑1 (MQ=255)
                                     aaCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGa                                                       <  1:2166042/91‑1 (MQ=255)
                                      aCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGaa                                                      <  2:1579540/91‑1 (MQ=255)
                                       ccTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAt                                                     >  1:3143344/1‑91 (MQ=255)
                                       ccTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAt                                                     >  1:1823324/1‑91 (MQ=255)
                                         tttATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAttt                                                   >  1:3665336/1‑91 (MQ=255)
                                         tttATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAttt                                                   <  1:3841661/91‑1 (MQ=255)
                                           tATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCg                                                 >  1:1268391/1‑91 (MQ=255)
                                           tATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCg                                                 >  1:612651/1‑91 (MQ=255)
                                             tCTATAGAACCATTATAGTACGAAAAAACCCTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTTGCt                                               >  2:3385597/1‑91 (MQ=255)
                                              cTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCtt                                              <  2:136300/91‑1 (MQ=255)
                                                ataGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTcc                                            >  2:142386/1‑91 (MQ=255)
                                                 taGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCa                                           >  2:3745399/1‑91 (MQ=255)
                                                  aGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCAt                                          <  1:3741413/91‑1 (MQ=255)
                                                   gAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATa                                         >  2:1139161/1‑91 (MQ=255)
                                                        aTTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATa                                    >  2:108766/1‑91 (MQ=255)
                                                        aTTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATa                                    <  2:1542953/91‑1 (MQ=255)
                                                         ttATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATaa                                   >  1:1141032/1‑91 (MQ=255)
                                                         ttATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATaa                                   >  1:2426546/1‑91 (MQ=255)
                                                           ataGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTTTGAATGGAATTTCGCTTCCATACCATAAAt                                 >  2:573621/1‑91 (MQ=255)
                                                           ataGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAAt                                 >  2:3683211/1‑91 (MQ=255)
                                                           ataGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAAt                                 >  1:2580973/1‑91 (MQ=255)
                                                            taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                >  2:2887496/1‑91 (MQ=255)
                                                            taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                <  1:2546175/91‑1 (MQ=255)
                                                            taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                >  2:2253036/1‑91 (MQ=255)
                                                            taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                >  2:2184928/1‑91 (MQ=255)
                                                              gTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATACGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACg                              >  2:1696125/1‑91 (MQ=255)
                                                               tACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGt                             <  2:1067317/91‑1 (MQ=255)
                                                                aCGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTg                            >  1:1658448/1‑91 (MQ=255)
                                                                 cGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGa                           <  2:165060/91‑1 (MQ=255)
                                                                  gAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGAc                          <  1:3417390/91‑1 (MQ=255)
                                                                   aaaaaaCACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGAcc                         <  2:2026050/91‑1 (MQ=255)
                                                                      gagcacTTATGAGGGTATGTACACATACAATCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCtt                      >  2:184023/4‑91 (MQ=18)
                                                                       aaCACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCttt                     >  2:3778818/1‑91 (MQ=255)
                                                                        gcacTTATGAGGGTATGTACACATACAATCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                    >  1:1404192/2‑91 (MQ=21)
                                                                        acacTTATGCTTGTATGTGTACATACCCTCATAAGTGCTCTATTATTGTGAATGGAATTTCGCTTCCATACGATACATACGTGACCCTTTa                    >  2:27678/1‑91 (MQ=255)
                                                                        acacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                    >  2:849386/1‑91 (MQ=255)
                                                                        acacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                    >  1:1867894/1‑91 (MQ=255)
                                                                          acTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                  >  1:2400399/1‑91 (MQ=255)
                                                                          acTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                  >  1:706600/1‑91 (MQ=255)
                                                                          acTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                  >  1:1647371/1‑91 (MQ=255)
                                                                           cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                 >  2:3809674/1‑91 (MQ=255)
                                                                           cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                 >  2:156687/1‑91 (MQ=255)
                                                                           cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                 <  1:2474646/91‑1 (MQ=255)
                                                                           cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                 >  1:692519/1‑91 (MQ=255)
                                                                           cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                 >  2:2846874/1‑91 (MQ=255)
                                                                           cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                 <  1:419342/91‑1 (MQ=255)
                                                                           cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTACCTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTAccc                 >  2:2564291/1‑91 (MQ=255)
                                                                            ttATGATTGTATGTGCACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCg                <  2:1325868/91‑1 (MQ=255)
                                                                              aTGATTGTGTGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaa              <  2:2789358/91‑1 (MQ=255)
                                                                               tGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaa             <  2:1057864/91‑1 (MQ=255)
                                                                               tGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaa             <  2:1060713/91‑1 (MQ=255)
                                                                               tGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaa             <  1:3157606/91‑1 (MQ=255)
                                                                               tGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaa             <  1:2242651/91‑1 (MQ=255)
                                                                                gATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaaa            <  2:1203663/91‑1 (MQ=255)
                                                                                  ttGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAtt          >  2:3197863/1‑91 (MQ=255)
                                                                                   tGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAttt         >  2:2346228/1‑91 (MQ=255)
                                                                                   tGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAttt         <  1:2894020/91‑1 (MQ=255)
                                                                                    gTATGTACACATACCCTCATAAGTGCTCTATTTTTTTGAATGGTATTTCGCTTCCATACGATAAATACGTCACCCTTTATCCGAAAATTTc        >  1:87505/1‑91 (MQ=255)
                                                                                     tATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCt       <  2:3242293/91‑1 (MQ=255)
                                                                                     tATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCt       <  1:1965464/91‑1 (MQ=255)
                                                                                      aTGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCtt      <  2:2650947/91‑1 (MQ=255)
                                                                                       tgtgTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCttt     >  1:376105/1‑91 (MQ=255)
                                                                                       tgtgTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCtaa     <  2:1972358/91‑3 (MQ=255)
                                                                                          gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaaa            >  1:1219919/1‑81 (MQ=255)
                                                                                          gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCa  <  2:526122/91‑1 (MQ=255)
                                                                                          gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCa  >  1:1003778/1‑91 (MQ=255)
                                                                                          |                                                                                          
AGTTAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCA  >  USA300TCH1516_ALE/1868952‑1869132

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