Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,044,680 A→C 9.5% intergenic (+29/+17) yfmC_1 → / ← USA300HOU_RS05165 Fe(3+)‑citrate‑binding protein YfmC/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,044,6800AC 9.5% 769.4 / 53.1 238intergenic (+29/+17)yfmC_1/USA300HOU_RS05165Fe(3+)‑citrate‑binding protein YfmC/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (116/99);  new base C (11/12);  total (127/111)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

GGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTACATAAC  >  USA300TCH1516_ALE/1044590‑1044770
                                                                                          |                                                                                          
ggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATaa                                                                                            >  1:1272041/1‑91 (MQ=255)
  aTTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAaga                                                                                          <  2:882543/91‑1 (MQ=255)
     aTCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGAAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTg                                                                                       <  2:2454386/91‑1 (MQ=255)
      tCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTgg                                                                                      >  2:2208438/1‑91 (MQ=255)
      tCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTgg                                                                                      <  2:2339811/91‑1 (MQ=255)
       cGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgc                                                                                     >  1:614609/1‑91 (MQ=255)
        gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                                    >  2:941040/1‑91 (MQ=255)
        gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                                    <  2:2683015/91‑1 (MQ=255)
          aaatGAAAATAGAGCATAAGATTTAGAAAAAATTGAAGAAAAAAAAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                                  <  2:3509642/88‑1 (MQ=255)
          aaGTGAAAGTATGGCAGAAGATTTAGAAACAATGGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                                  <  1:483634/91‑1 (MQ=255)
          aaGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                                  <  1:3526952/91‑1 (MQ=255)
           aGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATCAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTg                                                                                 <  1:3365655/91‑1 (MQ=255)
           aGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTg                                                                                 >  2:2970441/1‑91 (MQ=255)
           aGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTg                                                                                 >  2:3431937/1‑91 (MQ=255)
            gTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGt                                                                                <  2:3096092/91‑1 (MQ=255)
            gTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGt                                                                                <  2:122784/91‑1 (MQ=17)
            ctttaaGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGt                                                                                >  1:519468/5‑91 (MQ=255)
             tGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTa                                                                               <  1:894054/91‑1 (MQ=255)
              gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAt                                                                              <  1:390665/91‑1 (MQ=255)
              gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAt                                                                              >  1:317578/1‑91 (MQ=255)
                 aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTAtttt                                                                           <  2:2076859/91‑1 (MQ=17)
                     tGGCAGAAGATTTAGAAAAAATTGCGGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                                       <  2:3829762/91‑1 (MQ=255)
                     tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                                       >  1:2021175/1‑91 (MQ=255)
                     tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                                       >  1:1928552/1‑91 (MQ=255)
                     tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTatt                                                                       <  1:2395924/91‑3 (MQ=14)
                        cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                                    >  1:1576101/1‑91 (MQ=255)
                             gaTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTtta                                                               <  1:1460790/91‑1 (MQ=255)
                              aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                              >  2:3869249/1‑91 (MQ=255)
                               tttAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatt                                                             >  2:1523477/1‑91 (MQ=255)
                                ttAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatta                                                            <  1:813329/91‑1 (MQ=255)
                                 tAGAAAAAATTGCAGAAAAAGCAAACTAAAAATCCAGCGCTACTCGTAAATCATACAAGAGTGGCGCTGTCTTTCTAATTATGTCtatgat                                                           >  2:3277029/1‑91 (MQ=255)
                                 tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                           <  1:3056194/91‑1 (MQ=255)
                                 tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                           >  1:2752620/1‑91 (MQ=255)
                                 tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                           >  1:1090456/1‑91 (MQ=255)
                                 tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                           <  1:999321/91‑1 (MQ=255)
                                  aGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATc                                                          <  2:1019987/91‑1 (MQ=255)
                                  aGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATc                                                          >  2:1009729/1‑91 (MQ=255)
                                   gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                                         <  2:3738408/91‑1 (MQ=255)
                                   gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                                         >  1:3559951/1‑91 (MQ=255)
                                    aaaaaaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTg                                                        <  1:1002222/91‑1 (MQ=255)
                                        aaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCgg                                                    <  1:1321493/91‑1 (MQ=255)
                                         aTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGa                                                   <  1:3583245/91‑1 (MQ=255)
                                         aTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGa                                                   <  1:1726087/91‑1 (MQ=255)
                                           tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtg                                                 >  1:1699720/1‑91 (MQ=255)
                                           tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCATAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtg                                                 <  2:2239220/91‑1 (MQ=255)
                                            gCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgt                                                >  2:418082/1‑91 (MQ=255)
                                            gCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgt                                                >  1:2265886/1‑91 (MQ=255)
                                             cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                               >  2:1706683/1‑91 (MQ=255)
                                             cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                               >  1:1072183/1‑91 (MQ=255)
                                             cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                               <  1:3656924/91‑1 (MQ=255)
                                             cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                               <  1:115770/91‑1 (MQ=255)
                                             cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                               >  1:3450503/1‑91 (MQ=255)
                                             cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGAAGCGCTGTATTTTTAATTATGTTTATTATATGTCGGAtgtg                                               >  2:20330/1‑91 (MQ=255)
                                              aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                              >  2:834063/1‑91 (MQ=255)
                                              aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                              >  1:1601271/1‑91 (MQ=255)
                                               gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                                             >  1:2634149/1‑91 (MQ=255)
                                               gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                                             <  1:2313120/91‑1 (MQ=255)
                                                  aaaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTAAATTATGTTTATTATCTGTCGGATGTgatgat                                          >  1:3319406/1‑91 (MQ=255)
                                                   aaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGAtt                                         <  2:789093/91‑1 (MQ=255)
                                                     gCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTa                                       <  2:2752620/91‑1 (MQ=255)
                                                      cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                                      >  1:1535134/1‑91 (MQ=255)
                                                      cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                                      >  1:37880/1‑91 (MQ=255)
                                                       aaaataaaaataCAGTGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTGTTATCTGTCGGATGTGATGATTTAcc                                     >  1:1487927/1‑91 (MQ=255)
                                                        aattaaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                                    >  1:2957015/4‑91 (MQ=11)
                                                        aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                                    >  2:3697686/1‑91 (MQ=255)
                                                        aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                                    >  2:1920118/1‑91 (MQ=255)
                                                         aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                                   >  2:2894173/1‑91 (MQ=255)
                                                         aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                                   >  2:1746171/1‑91 (MQ=255)
                                                         aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                                   >  2:1955016/1‑91 (MQ=255)
                                                         aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                                   >  1:560801/1‑91 (MQ=255)
                                                           taaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaa                                 <  1:2410113/91‑1 (MQ=255)
                                                            aaaaataCAGCGTTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                                <  2:2567048/91‑1 (MQ=255)
                                                            aaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                                <  2:3230649/91‑1 (MQ=255)
                                                            aaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                                >  2:2860290/1‑91 (MQ=14)
                                                             aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                               <  1:3086149/91‑1 (MQ=255)
                                                             aaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                               >  2:1839212/1‑91 (MQ=17)
                                                              aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt                              <  1:779190/91‑1 (MQ=255)
                                                              aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCAGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt                              >  2:2052248/1‑91 (MQ=255)
                                                               aataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGtt                             >  1:2767901/1‑91 (MQ=255)
                                                                ataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGttt                            <  1:1781392/91‑1 (MQ=255)
                                                                ataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGttt                            >  2:491083/1‑91 (MQ=255)
                                                                 taCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTa                           >  1:1824279/1‑91 (MQ=17)
                                                                 taCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATATTTAATTATGTTTATTATCTGTCGGGTGTGATGATTTACCTGAAAGTTTa                           >  2:3597761/1‑91 (MQ=17)
                                                                  aCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAt                          >  1:2884257/1‑91 (MQ=255)
                                                                     gcgcCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTc                       >  1:1117954/1‑91 (MQ=18)
                                                                      cgcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCg                      >  2:3062353/1‑91 (MQ=255)
                                                                       gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa                     <  2:1576101/91‑1 (MQ=255)
                                                                       gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa                     <  1:923581/91‑1 (MQ=255)
                                                                         tACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAAt                   <  1:2894173/91‑1 (MQ=255)
                                                                         tACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAAt                   >  2:1380495/1‑91 (MQ=255)
                                                                          aCTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATa                  <  2:554125/91‑1 (MQ=255)
                                                                            tCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATaaa                >  2:1683994/1‑91 (MQ=255)
                                                                                aaaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTa            >  1:3569653/1‑91 (MQ=255)
                                                                                 aaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTaa           <  1:3402182/91‑1 (MQ=255)
                                                                                 aaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTaa           <  2:1211267/91‑1 (MQ=255)
                                                                                 aaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATGTaa           >  1:3487254/1‑91 (MQ=255)
                                                                                     atatAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTAc       <  2:1607942/91‑1 (MQ=255)
                                                                                     atatAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTAc       <  2:1173596/91‑1 (MQ=255)
                                                                                         aaGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTACATaa   >  2:1586425/1‑91 (MQ=255)
                                                                                          agagTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTACATAAc  >  2:1407299/1‑91 (MQ=255)
                                                                                          |                                                                                          
GGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTACATAAC  >  USA300TCH1516_ALE/1044590‑1044770

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: