Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 917,327 T→G 5.1% intergenic (+86/+244) sufB_2 → / ← USA300HOU_RS04545 FeS cluster assembly protein SufB/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE917,3270TG 5.1% 921.6 / 26.9 277intergenic (+86/+244)sufB_2/USA300HOU_RS04545FeS cluster assembly protein SufB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (139/123);  new base G (7/8);  total (146/131)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

GATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAA  >  USA300TCH1516_ALE/917237‑917417
                                                                                          |                                                                                          
  tAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACAtca                                                                                          >  2:1561306/1‑89 (MQ=255)
  tAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTa                                                                                          <  1:1487110/91‑1 (MQ=255)
   aGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAAGGGCTATTGAGACCATTTTGATACTGAATATAACAAAAAGCCACATTAc                                                                                         <  2:3930020/91‑1 (MQ=255)
    gCTTTAAACCGCGTTGTTAAGCCATTCTTGTCTTCCGAAAATGGCTATTGATACCATTTTGATACTGATTATAACAAAAAGCCACATTACt                                                                                        >  2:3189323/1‑91 (MQ=255)
    gCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACt                                                                                        >  1:1513223/1‑91 (MQ=255)
    gCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACAGTaat                                                                                        <  1:2941542/91‑3 (MQ=255)
     cTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACtg                                                                                       <  2:2591518/91‑1 (MQ=255)
     cTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACtg                                                                                       <  1:2918359/91‑1 (MQ=255)
        tAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACAGTAATGTgg                                                                                    <  2:2063443/91‑1 (MQ=255)
         aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAACAAGCCACATTACTGTGGc                                                                                   >  2:89545/1‑91 (MQ=255)
         aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGc                                                                                   >  2:1843247/1‑91 (MQ=255)
         aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGc                                                                                   <  2:501458/91‑1 (MQ=255)
         aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGc                                                                                   >  1:3767358/1‑91 (MQ=255)
         aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGc                                                                                   >  2:2643752/1‑91 (MQ=255)
         aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACAGTAATGTGGc                                                                                   <  1:1523278/91‑1 (MQ=255)
          aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGACTATACCAAAAAGCCACATTACCGTGGCt                                                                                  >  2:2720041/1‑91 (MQ=255)
           aCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtt                                                                                 <  1:1642373/91‑1 (MQ=255)
           aCCGCGTTGTTAAGCCATGCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtt                                                                                 <  2:463505/91‑1 (MQ=255)
              gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttt                                                                              <  1:1795614/91‑1 (MQ=255)
              gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttt                                                                              >  2:716368/1‑91 (MQ=255)
              gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttt                                                                              <  1:534121/91‑1 (MQ=255)
              gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttt                                                                              >  1:2307472/1‑91 (MQ=255)
               cgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtttttt                                                                             >  2:330727/1‑91 (MQ=255)
               cgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtttttt                                                                             >  1:4707/1‑91 (MQ=255)
               cgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtttttt                                                                             >  2:3375173/1‑91 (MQ=255)
               cgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtttttt                                                                             >  1:2019265/1‑91 (MQ=255)
               cgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGAGACCATTTTGATACTGAATAGAACAAAAAGCCACAGTAATGTGGCtttttg                                                                             <  2:1957287/91‑2 (MQ=38)
                 ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTg                                                                           >  1:2752856/1‑91 (MQ=255)
                 ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTg                                                                           >  1:516521/1‑91 (MQ=255)
                  tgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGt                                                                          <  1:930227/91‑1 (MQ=255)
                  tgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGt                                                                          <  1:1287568/91‑1 (MQ=255)
                    ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGttt                                                                        <  2:1057336/91‑1 (MQ=255)
                    ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGttt                                                                        >  1:2190988/1‑91 (MQ=255)
                      aaGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTta                                                                      >  1:1463937/1‑91 (MQ=255)
                        gCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTtata                                                                    <  1:2498123/91‑1 (MQ=255)
                        gCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTtata                                                                    <  1:1666285/91‑1 (MQ=255)
                               ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             <  2:2919776/91‑1 (MQ=255)
                               ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             <  1:3546363/91‑1 (MQ=255)
                               ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             <  2:1104353/91‑1 (MQ=255)
                               ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             >  2:3569637/1‑91 (MQ=255)
                               ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             >  1:1091770/1‑91 (MQ=255)
                                tGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             <  2:2261468/90‑1 (MQ=255)
                                 gACTTCCGAAAATGGCTATTGATACCATTTTGTTACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                           <  1:3373368/91‑1 (MQ=255)
                                 gACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                           >  1:2789439/1‑91 (MQ=255)
                                  aCTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCgg                                                          >  1:3179599/1‑91 (MQ=255)
                                   cTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGa                                                         >  1:3599626/1‑91 (MQ=255)
                                    ttCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGAt                                                        >  2:1293580/1‑91 (MQ=255)
                                         aaaaTGGCTATTGATACAACTTTGATACTGAATATAACAAAAAGCCACATTACCGAGGCTTTTTTTACCTTATAACTAAATCGGCTTcata                                                   >  2:1123777/1‑91 (MQ=255)
                                          aaaTGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatag                                                  <  2:1893474/91‑1 (MQ=255)
                                          aaaTGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatag                                                  >  1:1125696/1‑91 (MQ=255)
                                          aaaTGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatag                                                  >  1:630134/1‑91 (MQ=255)
                                           aaTGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAACCCGATTgataga                                                 >  2:1653942/1‑91 (MQ=255)
                                            aTGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatagat                                                <  1:2177422/91‑1 (MQ=255)
                                            aTGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatagat                                                <  2:3103043/91‑1 (MQ=255)
                                             tGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatagata                                               <  1:2249111/91‑1 (MQ=255)
                                             tGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatagata                                               <  2:1233355/91‑1 (MQ=255)
                                             tGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGACTTTTTTTGTTTTATAACTAAATCGGATTgatagata                                               <  1:3321681/91‑1 (MQ=255)
                                              ggCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATaa                                              >  1:3270702/1‑91 (MQ=255)
                                               gCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAg                                             >  1:3460742/1‑91 (MQ=255)
                                                 tATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTCTAACTAAATCGGATTGATAGATAAGCt                                           >  2:2485862/1‑91 (MQ=255)
                                                 tATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCt                                           >  2:3428294/1‑91 (MQ=255)
                                                  aTTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                          <  1:1975454/91‑1 (MQ=255)
                                                   ttGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCttt                                         >  1:2316753/1‑91 (MQ=255)
                                                        aCCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                    >  1:760959/1‑91 (MQ=255)
                                                         ccATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTAttt                                   >  2:2669128/1‑91 (MQ=255)
                                                         ccATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtt                                   >  1:190456/1‑91 (MQ=255)
                                                          cATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtta                                  <  1:685857/91‑1 (MQ=255)
                                                           aTTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttat                                 >  2:1876471/1‑91 (MQ=255)
                                                            ttttGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttatt                                >  1:610684/1‑91 (MQ=255)
                                                              ttGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattta                              <  2:3262299/91‑1 (MQ=255)
                                                              ttGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttatttat                             >  2:1725000/1‑91 (MQ=255)
                                                               tGATACTGAATATAACAAAAAGCCCCATTACTGTGGCTTTTTTTGTTTTATAACTCAATCGGCTTGATAGATAAGCTTTGTACttatttat                             >  2:1033077/1‑91 (MQ=255)
                                                               tGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttatttat                             >  1:2557378/1‑91 (MQ=255)
                                                               tGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTAATAGATAAGCTTTGTACttatttat                             >  1:3667115/1‑91 (MQ=255)
                                                                gATACTGAATATAACAAAAAGCCACATTCCTGTGACTTTTTTTGTTTTATAACTCAATAGGATTGATAGGTAAGCTTTGTACTTCTCtata                            >  1:1686782/1‑91 (MQ=255)
                                                                 aTACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTatat                           >  2:3869539/1‑91 (MQ=255)
                                                                 aTACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTatat                           >  1:1422220/1‑91 (MQ=255)
                                                                  tACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATc                          <  2:2681928/91‑1 (MQ=255)
                                                                    cTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAg                        >  1:3013338/1‑91 (MQ=255)
                                                                    cTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAg                        <  2:28126/91‑1 (MQ=255)
                                                                    cTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAg                        <  1:1155147/91‑1 (MQ=255)
                                                                    cTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAg                        >  2:3140346/1‑91 (MQ=255)
                                                                     tGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       <  2:344941/91‑1 (MQ=255)
                                                                     tGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       >  1:2083376/1‑91 (MQ=255)
                                                                     tGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       <  1:1541802/91‑1 (MQ=255)
                                                                      gAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGtt                      <  2:2290471/91‑1 (MQ=255)
                                                                      gAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGtt                      <  2:2219995/91‑1 (MQ=255)
                                                                        atatAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCg                    <  2:572253/91‑1 (MQ=255)
                                                                           taACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                 <  2:2979001/91‑1 (MQ=255)
                                                                           taACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                 <  1:2322458/91‑1 (MQ=255)
                                                                            aaCAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttt                <  2:3700824/91‑1 (MQ=255)
                                                                             aCAAGAAGCCTCAGGACTGTGGCTTTTTTTGTTTTATAACTATATCGGATTGATCGGTAAGCGTTGTGCTTGGCTAGATCAGGTCGCtttt               <  2:1397458/91‑1 (MQ=255)
                                                                             aCAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt               >  2:1534657/1‑91 (MQ=255)
                                                                              cAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt              >  2:672709/1‑91 (MQ=255)
                                                                                  aaGCCACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt          >  1:949731/1‑91 (MQ=255)
                                                                                    gCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTg        <  1:3429360/91‑1 (MQ=255)
                                                                                        caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt    <  1:3294361/91‑1 (MQ=255)
                                                                                         aTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATACATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTa   <  1:3946579/91‑1 (MQ=255)
                                                                                          ttACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaa  <  2:1677054/91‑1 (MQ=255)
                                                                                          ttACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaa  <  2:947441/91‑1 (MQ=255)
                                                                                          |                                                                                          
GATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAA  >  USA300TCH1516_ALE/917237‑917417

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: