Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,869,051 C→A 5.5% intergenic (‑458/+12) 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,0510CA 5.5% 701.4 / 24.0 201intergenic (‑458/+12)ytnP/trmBputative quorum‑quenching lactonase YtnP/tRNA (guanine‑N(7)‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (105/85);  new base A (8/3);  total (113/88)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 211 total aligned reads displayed.

CTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCATATTCTGTT  >  USA300TCH1516_ALE/1868961‑1869141
                                                                                          |                                                                                          
cTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATAccc                                                                                            <  2:3776974/91‑1 (MQ=255)
 tttttCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCt                                                                                           >  1:1463208/1‑91 (MQ=255)
      ccATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATaa                                                                                      >  1:2940482/1‑91 (MQ=255)
       cATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATACAATCATAAg                                                                                     <  1:3581673/91‑1 (MQ=25)
        aTCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATACAATCATAAGt                                                                                    <  1:1762090/91‑1 (MQ=21)
         tCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTg                                                                                   <  1:255907/91‑1 (MQ=255)
         tCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTg                                                                                   <  1:2928617/91‑1 (MQ=255)
         tCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTg                                                                                   <  2:1204927/91‑1 (MQ=255)
         tCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATACAATCATAAGTg                                                                                   <  2:2590063/91‑1 (MQ=21)
          ctgtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGc                                                                                  <  2:1031563‑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)
             tccgatttATTATAGAACCTTTATCTATAGAACCATTATAGGACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGctct                                                                               <  2:3410390‑M2/83‑1 (MQ=255)
                gatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTAtt                                                                            >  2:1699765‑M2/6‑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                                                                     <  1:1583356/91‑1 (MQ=255)
                       tataGAACCGTTATCTATAGAGCCATTATCGTACGAAAAAACACGTATGATGGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGa                                                                     <  2:2837146/91‑1 (MQ=255)
                        ataGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGaa                                                                    >  1:1058642/1‑91 (MQ=255)
                        ataGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGaa                                                                    <  1:1242907/91‑1 (MQ=255)
                          aGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATg                                                                  <  1:2137093/91‑1 (MQ=255)
                           gAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATgg                                                                 <  2:3745174/91‑1 (MQ=255)
                            aaCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGa                                                                <  2:1085974/91‑1 (MQ=255)
                            aaCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGa                                                                <  2:248023/91‑1 (MQ=255)
                            aaCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGa                                                                <  1:2166042/91‑1 (MQ=255)
                              ccTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAt                                                              >  1:1823324/1‑91 (MQ=255)
                              ccTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAt                                                              >  1:3143344/1‑91 (MQ=255)
                                tttATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAttt                                                            >  1:3665336/1‑91 (MQ=255)
                                tttATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAttt                                                            <  1:3842457/91‑1 (MQ=255)
                                  tATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCg                                                          >  1:612651/1‑91 (MQ=255)
                                  tATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCg                                                          <  1:3411192/91‑1 (MQ=255)
                                  tATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCg                                                          >  2:1105620/1‑91 (MQ=255)
                                   aTCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGc                                                         >  1:1915984/1‑91 (MQ=255)
                                    tCTATAGAACCATTATAGTACGAAAAAACCCTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTTGCt                                                        >  2:3385597/1‑91 (MQ=255)
                                    tCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTTGCt                                                        >  2:1061346/1‑91 (MQ=255)
                                     cTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCtt                                                       >  2:2360436/1‑91 (MQ=255)
                                     cTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCtt                                                       <  2:136352/91‑1 (MQ=255)
                                        taGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCa                                                    <  2:814999/91‑1 (MQ=255)
                                         aGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCAt                                                   <  1:3741413/91‑1 (MQ=255)
                                              cATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGAt                                              >  1:305629/1‑91 (MQ=255)
                                               aTTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATa                                             <  2:1542953/91‑1 (MQ=255)
                                                ttATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATaa                                            >  1:1141032/1‑91 (MQ=255)
                                                  ataGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAAt                                          >  2:3683211/1‑91 (MQ=255)
                                                  ataGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAAt                                          >  2:573555/1‑91 (MQ=255)
                                                   taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                         >  2:2184735/1‑91 (MQ=255)
                                                   taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                         >  2:2887496/1‑91 (MQ=255)
                                                    aGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATAc                                        >  2:409296/1‑91 (MQ=255)
                                                    aGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATAc                                        >  1:2602933/1‑91 (MQ=255)
                                                      tACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGt                                      <  1:2363880/91‑1 (MQ=255)
                                                       aCGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTg                                     >  1:1657779/1‑91 (MQ=255)
                                                        cGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGa                                    <  2:165060/91‑1 (MQ=255)
                                                         gAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGAc                                   >  1:992424/1‑91 (MQ=255)
                                                          aaaaaaCACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGAcc                                  <  2:2026050/91‑1 (MQ=255)
                                                            agagcacTTATGAGGGTATGTACACATACAATCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCt                                >  2:2776165/5‑91 (MQ=18)
                                                             gagcacTTATGAGGGTATGTACACATACAATCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCtt                               >  1:3449640/4‑91 (MQ=18)
                                                             gagcacTTATGAGGGTATGTACACATACAATCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCtt                               >  2:184023/4‑91 (MQ=18)
                                                             aaaCACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCtt                               <  1:1970391/91‑1 (MQ=255)
                                                              agcacTTATGAGGGTATGTACACATACAATCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCttt                              >  1:2459215/3‑91 (MQ=21)
                                                               acacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGTATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                             >  2:27659/1‑91 (MQ=255)
                                                               acacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                             >  1:1867894/1‑91 (MQ=255)
                                                               acacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                             >  1:1466569/1‑91 (MQ=255)
                                                                cacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTAt                            >  1:251621/1‑91 (MQ=255)
                                                                 acTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                           >  1:2403414/1‑91 (MQ=255)
                                                                 acTTATGAGGGTATGTACACATACAATCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                           >  1:2304430/1‑91 (MQ=21)
                                                                  cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                          >  1:733480/1‑91 (MQ=255)
                                                                  cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                          >  2:2846874/1‑91 (MQ=255)
                                                                  cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                          <  1:419342/91‑1 (MQ=255)
                                                                  cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                          >  2:156687/1‑91 (MQ=255)
                                                                  cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                          >  1:692519/1‑91 (MQ=255)
                                                                   ttATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCg                         >  2:1671930/1‑91 (MQ=255)
                                                                   ttATGATTGTATGTGCACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCg                         <  2:1325868/91‑1 (MQ=255)
                                                                      tGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaa                      <  2:1058642/91‑1 (MQ=255)
                                                                      tGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaa                      <  1:3157606/91‑1 (MQ=255)
                                                                        aTTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAt                    <  1:2444703/91‑1 (MQ=255)
                                                                         ttGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAtt                   >  2:3197863/1‑91 (MQ=255)
                                                                          tGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAttt                  <  2:2343369/91‑1 (MQ=255)
                                                                           gTATGTACACATACCCTCATAAGTGCTCTATTTTTTTGAATGGTATTTCGCTTCCATACGATAAATACGTCACCCTTTATCCGAAAATTTc                 >  1:87505/1‑91 (MQ=255)
                                                                            tATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCt                <  2:3242293/91‑1 (MQ=255)
                                                                             aTGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCtt               <  1:1327898/91‑1 (MQ=255)
                                                                             aTGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCtt               <  2:2650947/91‑1 (MQ=255)
                                                                              tgtgTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCttt              <  2:2460473/91‑1 (MQ=255)
                                                                              tgtgTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCttt              >  1:376105/1‑91 (MQ=255)
                                                                              tgtgTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCtaa              <  2:1972358/91‑3 (MQ=255)
                                                                                tgtACATACCCTCATAAGTGCTCTATTTTTGTGACTGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTc            >  2:3791759/1‑91 (MQ=255)
                                                                                 gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaaa                     >  1:1219919/1‑81 (MQ=255)
                                                                                 gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCa           >  1:1003778/1‑91 (MQ=255)
                                                                                 gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCa           <  2:526272/91‑1 (MQ=255)
                                                                                 acaCATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCa           >  1:2808839/3‑91 (MQ=255)
                                                                                  tACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTcgt          <  2:1452449/91‑3 (MQ=255)
                                                                                  tACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCat          >  1:966502/1‑91 (MQ=255)
                                                                                   aCATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCata         <  2:1577748/91‑1 (MQ=255)
                                                                                   aCATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCata         <  2:1329849/91‑1 (MQ=255)
                                                                                   aCATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCata         >  1:500596/1‑91 (MQ=255)
                                                                                   aCATACCCTCAGAAGTGCTCTATCTTTGTGAATGGAATTTCGCTTCCATACGATAAATAGGTGACCCTTTATCCGAAAATTTCTTTTCata         <  2:2965365/91‑1 (MQ=255)
                                                                                      tACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCATATTc      >  2:909028/1‑91 (MQ=255)
                                                                                       aCCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCATATTCt     >  1:1287380/1‑91 (MQ=255)
                                                                                        cccTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCATATTCTg    <  2:220793/91‑1 (MQ=255)
                                                                                        cccTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCATATTCTg    >  2:1794510/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCATATTCTGTT  >  USA300TCH1516_ALE/1868961‑1869141

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: