Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,869,043 T→C 6.7% intergenic (‑450/+20) 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,0430TC 6.7% 649.5 / 24.0 204intergenic (‑450/+20)ytnP/trmBputative quorum‑quenching lactonase YtnP/tRNA (guanine‑N(7)‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (103/87);  new base C (9/5);  total (112/92)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TTAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCAT  >  USA300TCH1516_ALE/1868954‑1869133
                                                                                         |                                                                                          
 tAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTAc                                                                                          <  1:504986/91‑1 (MQ=255)
 tAAAAACTTTTTCCATCAGTCCGATTTATTAGAGAACCTTTATCTATAGAACCTTTATAGTACGAAAAAACACTTATGATTGTATGTGTAc                                                                                          <  1:3103704/91‑1 (MQ=255)
 gtttttaacttttcatctgtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTAc                                                                                          >  1:2962357‑M2/28‑91 (MQ=255)
  aaaaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACa                                                                                         <  2:1184094/91‑1 (MQ=255)
   aaaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACAt                                                                                        >  1:365179/1‑91 (MQ=255)
    tttaacttttcatctgtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATa                                                                                       <  2:2596281‑M2/67‑1 (MQ=255)
    aaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATa                                                                                       >  2:190760/1‑91 (MQ=255)
    aaaCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATa                                                                                       >  2:188403/1‑91 (MQ=255)
      aCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATAcc                                                                                     <  2:3055553/91‑1 (MQ=255)
      aCTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATAcc                                                                                     <  2:2465235/91‑1 (MQ=255)
         ttttCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATAcaatc                                                                                  <  1:2953763/91‑5 (MQ=37)
               aTCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGt                                                                            <  2:567786/91‑1 (MQ=255)
               aTCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGAGGGTATGTACACATACAATCATAAGt                                                                            <  1:1762090/91‑1 (MQ=21)
                tCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTg                                                                           <  2:1204927/91‑1 (MQ=255)
                tCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTg                                                                           <  1:2928617/91‑1 (MQ=255)
                 ctgtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGc                                                                          <  2:1030305‑M2/80‑1 (MQ=255)
                 ctgtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGc                                                                          <  2:3589371‑M2/80‑1 (MQ=255)
                  aGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGct                                                                         <  1:1506447/91‑1 (MQ=255)
                   gtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGctc                                                                        <  2:2426359‑M2/82‑1 (MQ=255)
                   gtccgatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGctc                                                                        <  1:314651‑M2/82‑1 (MQ=255)
                    tccgatttATTATAGAACCTTTATCTATAGAACCATTATAGGACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGctct                                                                       <  2:3410390‑M2/83‑1 (MQ=255)
                       gatttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTAtt                                                                    >  2:1699765‑M2/6‑91 (MQ=255)
                        atttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTAttt                                                                   >  1:97793‑M2/5‑91 (MQ=255)
                        atttATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTAttt                                                                   >  1:3658307‑M2/5‑91 (MQ=255)
                            attatAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTtgt                                                               >  1:583472/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:1242907/91‑1 (MQ=255)
                               ataGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGaa                                                            >  1:1574265/1‑91 (MQ=255)
                                 aGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATg                                                          <  1:2137093/91‑1 (MQ=255)
                                  gAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATgg                                                         <  2:3745174/91‑1 (MQ=255)
                                   aaCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGa                                                        <  1:2166042/91‑1 (MQ=255)
                                   aaCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGa                                                        <  2:248023/91‑1 (MQ=255)
                                    aCCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGaa                                                       <  2:1579540/91‑1 (MQ=255)
                                     ccTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAt                                                      >  1:1823324/1‑91 (MQ=255)
                                       tttATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAAttt                                                    <  1:3841661/91‑1 (MQ=255)
                                         tATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCg                                                  >  2:1105620/1‑91 (MQ=255)
                                         tATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCg                                                  >  1:1268391/1‑91 (MQ=255)
                                           tCTATAGAACCATTATAGTACGAAAAAACCCTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTTGCt                                                >  2:3385597/1‑91 (MQ=255)
                                            cTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCtt                                               <  2:136300/91‑1 (MQ=255)
                                            cTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCtt                                               >  2:221258/1‑91 (MQ=255)
                                               taGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCa                                            >  1:309065/1‑91 (MQ=255)
                                                     cATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGAt                                      >  1:305629/1‑91 (MQ=255)
                                                      aTTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATa                                     <  1:664262/91‑1 (MQ=255)
                                                      aTTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATa                                     >  1:2820848/1‑91 (MQ=255)
                                                      aTTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATa                                     <  2:1542953/91‑1 (MQ=255)
                                                       ttATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATaa                                    >  1:498070/1‑91 (MQ=255)
                                                       ttATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATaa                                    >  1:1141032/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                                 >  2:2184735/1‑91 (MQ=255)
                                                          taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                 >  1:3159872/1‑91 (MQ=255)
                                                          taGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATa                                 <  1:2546175/91‑1 (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:3479294/1‑91 (MQ=255)
                                                                      acacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                     >  1:1466569/1‑91 (MQ=255)
                                                                      acacTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTa                     >  1:1867957/1‑91 (MQ=255)
                                                                        acTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                   >  1:1647371/1‑91 (MQ=255)
                                                                        acTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                   >  1:706600/1‑91 (MQ=255)
                                                                        acTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATc                   >  1:2403414/1‑91 (MQ=255)
                                                                         cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATcc                  >  2:2846874/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:3733341/1‑91 (MQ=255)
                                                                         cTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTACCTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTAccc                  >  2:2564291/1‑91 (MQ=255)
                                                                          ttATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCg                 >  2:1671930/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              <  1:2242651/91‑1 (MQ=255)
                                                                             tGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaa              <  2:1060713/91‑1 (MQ=255)
                                                                               aTTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAt            <  1:2444703/91‑1 (MQ=255)
                                                                                ttGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAtt           >  2:3197863/1‑91 (MQ=255)
                                                                                 tGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAAttt          >  2:2346228/1‑91 (MQ=255)
                                                                                  gTATGTACACATACCCTCATAAGTGCTCTATTTTTTTGAATGGTATTTCGCTTCCATACGATAAATACGTCACCCTTTATCCGAAAATTTc         >  1:87505/1‑91 (MQ=255)
                                                                                   tATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCt        <  1:1965464/91‑1 (MQ=255)
                                                                                   tATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCt        <  2:3242293/91‑1 (MQ=255)
                                                                                   tATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCt        <  1:3588607/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)
                                                                                      gtgtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCtttt     <  2:1141030/91‑1 (MQ=255)
                                                                                       tgtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTCCCATACGATAAATACGTGACCCTTTATCCTAAAATTTCTTTTc    >  1:3885004/1‑91 (MQ=255)
                                                                                        gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGaaaa             >  1:1219919/1‑81 (MQ=255)
                                                                                        gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCa   <  2:526122/91‑1 (MQ=255)
                                                                                        gtACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCa   >  1:1003778/1‑91 (MQ=255)
                                                                                         tACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTcgt  <  2:1452449/91‑3 (MQ=255)
                                                                                         |                                                                                          
TTAAAAACTTTTTCCATCAGTCCGATTTATTATAGAACCTTTATCTATAGAACCATTATAGTACGAAAAAACACTTATGATTGTATGTGTACATACCCTCATAAGTGCTCTATTTTTGTGAATGGAATTTCGCTTCCATACGATAAATACGTGACCCTTTATCCGAAAATTTCTTTTCAT  >  USA300TCH1516_ALE/1868954‑1869133

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: