Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,044,677 A→T 9.6% intergenic (+26/+20) 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,6770AT 9.6% 761.7 / 50.7 236intergenic (+26/+20)yfmC_1/USA300HOU_RS05165Fe(3+)‑citrate‑binding protein YfmC/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (116/96);  new base T (11/12);  total (128/108)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

GGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTACAT  >  USA300TCH1516_ALE/1044588‑1044767
                                                                                         |                                                                                          
  ggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATaa                                                                                         >  1:1272041/1‑91 (MQ=255)
    aTTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAaga                                                                                       <  2:882543/91‑1 (MQ=255)
         cGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgc                                                                                  >  1:614609/1‑91 (MQ=255)
         cGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgc                                                                                  >  1:1780156/1‑91 (MQ=255)
          gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                                 >  2:941040/1‑91 (MQ=255)
            aaatGAAAATAGAGCATAAGATTTAGAAAAAATTGAAGAAAAAAAAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                               <  2:3509642/88‑1 (MQ=255)
            aaGTGAAAGTATGGCAGAAGATTTAGAAACAATGGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                               <  1:483634/91‑1 (MQ=255)
            aaGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                               <  2:3461636/91‑1 (MQ=255)
            aaGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                               <  1:3526952/91‑1 (MQ=255)
             aGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATCAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTg                                                                              <  1:3365655/91‑1 (MQ=255)
              gTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGt                                                                             <  2:3096092/91‑1 (MQ=255)
              ctttaaGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGt                                                                             >  1:519468/5‑91 (MQ=255)
               tGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTa                                                                            <  1:894054/91‑1 (MQ=255)
                gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTAt                                                                           <  1:3201735/91‑1 (MQ=17)
                   aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                                        >  2:2695378/1‑91 (MQ=255)
                    gTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAttttt                                                                       >  2:2225589/1‑91 (MQ=255)
                    gTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAttttt                                                                       >  2:2103518/1‑91 (MQ=255)
                       tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                                    >  1:2021175/1‑91 (MQ=255)
                       tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTatt                                                                    <  1:2396622/91‑3 (MQ=14)
                          cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                                 <  2:1860672/91‑1 (MQ=255)
                          cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                                 <  2:1720304/91‑1 (MQ=255)
                          cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                                 >  1:1576101/1‑91 (MQ=255)
                          cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                                 >  2:2433623/1‑91 (MQ=255)
                           agaagaTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGCAAATCATATAAGAGTGGCGCTGTATTTTTAATTATg                                                                >  2:3063344/1‑91 (MQ=255)
                           agaagaGTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATg                                                                <  2:3264871/91‑1 (MQ=255)
                                aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                           >  2:3869249/1‑91 (MQ=255)
                                aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                           >  2:3867494/1‑91 (MQ=255)
                                aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                           <  1:3331038/91‑1 (MQ=255)
                                 tttAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatt                                                          >  2:1523477/1‑91 (MQ=255)
                                  ttAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatta                                                         <  1:813329/91‑1 (MQ=255)
                                   tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                        >  1:1090456/1‑91 (MQ=255)
                                   tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                        <  1:999321/91‑1 (MQ=255)
                                    aGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATc                                                       >  2:1675034/1‑91 (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)
                                     gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                                      <  2:2257837/91‑1 (MQ=255)
                                      aaaaaaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTg                                                     >  1:1851125/1‑91 (MQ=255)
                                      aaaaaaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTg                                                     <  1:1002222/91‑1 (MQ=255)
                                          aaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCgg                                                 <  1:1321493/91‑1 (MQ=255)
                                           aTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGa                                                <  1:1726087/91‑1 (MQ=255)
                                           aTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGa                                                <  1:3583245/91‑1 (MQ=255)
                                            ttGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAt                                               >  1:1523611/1‑91 (MQ=255)
                                             tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtg                                              >  1:1699720/1‑91 (MQ=255)
                                             tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtg                                              <  1:1422922/91‑1 (MQ=255)
                                              gCAGAAAAAGCAAAATAAAAATACAGCGCTACACGTAAATCATATAAGAGTGGCGCTGTATTTCTAATTATGTATATTATCTGTCAGAtgt                                             >  2:1226836/1‑91 (MQ=39)
                                               cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                            <  1:115770/91‑1 (MQ=255)
                                               cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                            >  1:1072183/1‑91 (MQ=255)
                                               cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                            >  1:3450503/1‑91 (MQ=255)
                                               cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                            >  1:23701/1‑91 (MQ=255)
                                               cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                            <  1:3656924/91‑1 (MQ=255)
                                               cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGAAGCGCTGTATTTTTAATTATGTTTATTATATGTCGGAtgtg                                            >  2:20330/1‑91 (MQ=255)
                                                aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                           >  2:834063/1‑91 (MQ=255)
                                                aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                           >  1:1601271/1‑91 (MQ=255)
                                                aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                           >  2:2481024/1‑91 (MQ=255)
                                                 gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                                          >  1:2634149/1‑91 (MQ=255)
                                                 gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                                          >  2:1588685/1‑91 (MQ=255)
                                                    aaaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTgatgat                                       <  1:2238579/91‑1 (MQ=255)
                                                     aaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGAtt                                      <  2:789093/91‑1 (MQ=255)
                                                       gCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTa                                    <  2:2752620/91‑1 (MQ=255)
                                                        cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                                   >  1:37880/1‑91 (MQ=255)
                                                        cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                                   >  2:1837850/1‑91 (MQ=255)
                                                         aaaataaaaataCAGTGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTGTTATCTGTCGGATGTGATGATTTAcc                                  >  1:1487927/1‑91 (MQ=255)
                                                         aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                                  >  1:2655127/1‑91 (MQ=255)
                                                          aattaaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                                 >  1:2957015/4‑91 (MQ=11)
                                                          aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                                 >  2:3697686/1‑91 (MQ=255)
                                                          aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                                 >  2:1920460/1‑91 (MQ=255)
                                                           aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                                >  2:2894173/1‑91 (MQ=255)
                                                           aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                                >  1:560801/1‑91 (MQ=255)
                                                            ttaaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGa                               >  2:456600/2‑91 (MQ=17)
                                                              aaaaataCAGCGTTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                             <  2:2567048/91‑1 (MQ=255)
                                                              aaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                             <  2:3230649/91‑1 (MQ=255)
                                                               aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                            <  2:1998588/91‑1 (MQ=255)
                                                               aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                            <  1:3086149/91‑1 (MQ=255)
                                                               aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                            <  1:2876193/91‑1 (MQ=255)
                                                                aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt                           <  1:779190/91‑1 (MQ=255)
                                                                aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt                           >  2:2864852/1‑91 (MQ=255)
                                                                aaataCAGCGCTACTCGTAAATCATATAAGAGCGGCGCTGTACCTCTCATTATGTTTATTCCCTGTCGGATGTGATGATTTACCTGAAAgc                           >  2:2052286/1‑90 (MQ=255)
                                                                 aataCAGCGCTACTCGTAAATCATGTAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGtt                          >  1:2767664/1‑91 (MQ=255)
                                                                 aataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGtt                          <  1:2951011/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)
                                                                    aCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAt                       >  1:401063/1‑91 (MQ=255)
                                                                      aGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAttt                     <  1:3069106/91‑1 (MQ=255)
                                                                       gcgcCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTc                    >  1:1117954/1‑91 (MQ=18)
                                                                        cgcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATGTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCg                   >  2:3061183/1‑91 (MQ=255)
                                                                         gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa                  <  1:923581/91‑1 (MQ=255)
                                                                         gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa                  <  2:1576101/91‑1 (MQ=255)
                                                                           tACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAAt                >  2:1380495/1‑91 (MQ=255)
                                                                            aCTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATa               <  2:554125/91‑1 (MQ=255)
                                                                               cGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAt            >  2:1706234/1‑91 (MQ=255)
                                                                                 tAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAttt          <  1:228938/91‑1 (MQ=255)
                                                                                 tAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAttt          >  1:3439906/1‑91 (MQ=255)
                                                                                 tAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAttt          >  1:348880/1‑91 (MQ=255)
                                                                                  aaaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTa         >  1:3569653/1‑91 (MQ=255)
                                                                                   aaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTaa        <  2:1211267/91‑1 (MQ=255)
                                                                                       atatAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTAc    <  2:1173596/91‑1 (MQ=255)
                                                                                         |                                                                                          
GGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTACAT  >  USA300TCH1516_ALE/1044588‑1044767

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: