Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,044,674 C→G 9.4% intergenic (+23/+23) 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,6740CG 9.4% 790.8 / 51.4 240intergenic (+23/+23)yfmC_1/USA300HOU_RS05165Fe(3+)‑citrate‑binding protein YfmC/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (119/97);  new base G (11/12);  total (131/109)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

CAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTA  >  USA300TCH1516_ALE/1044584‑1044764
                                                                                          |                                                                                          
cAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATc                                                                                            >  1:1781920/1‑91 (MQ=255)
  aGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCat                                                                                          >  1:1461648/1‑91 (MQ=255)
     gggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCAtata                                                                                       >  1:2683015/1‑91 (MQ=255)
      ggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATaa                                                                                      >  2:2750235/1‑91 (MQ=255)
      ggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATaa                                                                                      >  1:1272041/1‑91 (MQ=255)
        aTTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAaga                                                                                    <  2:882543/91‑1 (MQ=255)
          tATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTGAATCATATAAGAGt                                                                                  <  1:1713154/91‑1 (MQ=255)
          tATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGt                                                                                  >  1:91218/1‑91 (MQ=255)
            tCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTgg                                                                                <  2:2339811/91‑1 (MQ=255)
            tCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTgg                                                                                >  2:2208438/1‑91 (MQ=255)
             cGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgc                                                                               >  1:613546/1‑91 (MQ=255)
              gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                              <  1:858864/91‑1 (MQ=255)
              gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                              >  2:941040/1‑91 (MQ=255)
              gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                              <  2:2683015/91‑1 (MQ=255)
                aaatGAAAATAGAGCATAAGATTTAGAAAAAATTGAAGAAAAAAAAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                            <  2:3509642/88‑1 (MQ=255)
                aaGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                            <  2:3461636/91‑1 (MQ=255)
                 aGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTg                                                                           <  1:437326/91‑1 (MQ=255)
                  gTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGt                                                                          <  2:3096092/91‑1 (MQ=255)
                  gTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGt                                                                          <  2:122784/91‑1 (MQ=17)
                   tGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTa                                                                         <  1:2890280/91‑1 (MQ=255)
                    gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAt                                                                        <  1:390665/91‑1 (MQ=255)
                    gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAt                                                                        >  1:317563/1‑91 (MQ=255)
                    gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTAt                                                                        <  1:3201735/91‑1 (MQ=17)
                     aaaGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtt                                                                       >  1:3543564/1‑91 (MQ=255)
                       aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                                     >  1:1665485/1‑91 (MQ=255)
                       aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                                     >  1:151632/1‑91 (MQ=255)
                       aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTAtttt                                                                     <  1:3631805/91‑1 (MQ=17)
                        gTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAttttt                                                                    >  2:2103518/1‑91 (MQ=255)
                        gTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAttttt                                                                    <  1:797667/91‑1 (MQ=255)
                           tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                                 >  1:2021175/1‑91 (MQ=255)
                           tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                                 >  1:1928552/1‑91 (MQ=255)
                           tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTatt                                                                 <  1:358/91‑3 (MQ=14)
                              cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                              >  1:1576102/1‑91 (MQ=255)
                              cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                              <  2:1720304/91‑1 (MQ=255)
                              cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                              <  2:1860672/91‑1 (MQ=255)
                               agaagaGTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATg                                                             <  2:3264871/91‑1 (MQ=255)
                                gaagaTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGt                                                            >  2:3078521/1‑91 (MQ=255)
                                    aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                        <  1:3331038/91‑1 (MQ=255)
                                    aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                        >  2:3867494/1‑91 (MQ=255)
                                    aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                        >  2:3869249/1‑91 (MQ=255)
                                     tttAGAAAAAATTGCAGATAAAGCAGGTTAAAAATACAGTGCTACTCGTAAATCATAGAAGTGTGGCGCTGTAATTTTCATTATGTttatt                                                       >  2:1523630/1‑91 (MQ=255)
                                      ttAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatta                                                      >  1:3773094/1‑91 (MQ=255)
                                       tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                     <  1:999628/91‑1 (MQ=255)
                                       tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                     >  1:77490/1‑91 (MQ=255)
                                       tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                     >  1:1090456/1‑91 (MQ=255)
                                       tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                     <  1:50221/91‑1 (MQ=255)
                                       tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                                     >  2:280505/1‑91 (MQ=255)
                                         gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                                   <  2:2257837/91‑1 (MQ=255)
                                         gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                                   <  2:3738408/91‑1 (MQ=255)
                                          aaaaaaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTg                                                  <  1:1002222/91‑1 (MQ=255)
                                           aaaaaTTGCATAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGt                                                 <  2:1202913/91‑1 (MQ=255)
                                              aaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCgg                                              <  1:1321493/91‑1 (MQ=255)
                                                ttGCAAAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAt                                            <  2:1211906/91‑1 (MQ=255)
                                                 tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtg                                           <  1:1422922/91‑1 (MQ=255)
                                                 tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGCAtg                                           >  1:3496234/1‑91 (MQ=255)
                                                   cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                         <  1:115770/91‑1 (MQ=255)
                                                   cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                         >  1:1072183/1‑91 (MQ=255)
                                                   cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                         >  1:23701/1‑91 (MQ=255)
                                                   cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGAAGCGCTGTATTTTTAATTATGTTTATTATATGTCGGAtgtg                                         >  2:20330/1‑91 (MQ=255)
                                                    aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                        >  2:2481024/1‑91 (MQ=255)
                                                    aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                        >  2:834063/1‑91 (MQ=255)
                                                     gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                                       <  1:2313120/91‑1 (MQ=255)
                                                     gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                                       >  1:2634149/1‑91 (MQ=255)
                                                        aaaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTgatgat                                    <  1:2238579/91‑1 (MQ=255)
                                                         aaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGAtt                                   <  2:789093/91‑1 (MQ=255)
                                                           gCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTa                                 >  2:1495103/1‑91 (MQ=255)
                                                             aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                               <  2:2960378/91‑1 (MQ=255)
                                                             aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                               >  2:163974/1‑91 (MQ=255)
                                                             aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTCAcc                               >  1:3008298/1‑91 (MQ=255)
                                                              aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                              >  2:1920118/1‑91 (MQ=255)
                                                              aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                              >  2:3697686/1‑91 (MQ=255)
                                                               aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                             >  2:1746171/1‑91 (MQ=255)
                                                                ttaaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGa                            >  2:456600/2‑91 (MQ=17)
                                                                 taaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaa                           <  1:2410113/91‑1 (MQ=255)
                                                                  aaaaataCAGCGTTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                          <  2:2567048/91‑1 (MQ=255)
                                                                  aaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                          <  2:3230649/91‑1 (MQ=255)
                                                                  aaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                          <  1:982298/91‑1 (MQ=255)
                                                                   aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                         <  1:2876193/91‑1 (MQ=255)
                                                                   aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                         <  2:1998588/91‑1 (MQ=255)
                                                                   aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg                         <  1:536380/91‑1 (MQ=255)
                                                                    aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCAGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt                        >  2:2052248/1‑91 (MQ=255)
                                                                     aataCAGCGCTACTCGTAAATCATGTAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGtt                       >  1:2767664/1‑91 (MQ=255)
                                                                     aataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGtt                       <  1:2951011/91‑1 (MQ=255)
                                                                      ataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGttt                      >  2:491083/1‑91 (MQ=255)
                                                                      ataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGttt                      <  1:1776503/91‑1 (MQ=255)
                                                                       taCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTa                     >  1:1824279/1‑91 (MQ=17)
                                                                       taCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATATTTAATTATGTTTATTATCTGTCGGGTGTGATGATTTACCTGAAAGTTTa                     >  2:3597761/1‑91 (MQ=17)
                                                                          aGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAttt                  <  1:3069106/91‑1 (MQ=255)
                                                                           gcgcCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTc                 >  1:1117954/1‑91 (MQ=18)
                                                                            cgcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCg                >  2:3062353/1‑91 (MQ=255)
                                                                             gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa               <  2:1576102/91‑1 (MQ=255)
                                                                                aCTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATa            <  2:554125/91‑1 (MQ=255)
                                                                                aCTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATa            <  1:15473/91‑1 (MQ=255)
                                                                                   cGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAt         >  2:1706234/1‑91 (MQ=255)
                                                                                   cGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAt         <  2:2878686/91‑1 (MQ=255)
                                                                                     tAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAttt       >  1:3439925/1‑91 (MQ=255)
                                                                                     tAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAttt       >  2:1015883/1‑91 (MQ=255)
                                                                                      aaaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTa      >  1:3569351/1‑91 (MQ=255)
                                                                                       aaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTaa     <  1:3402182/91‑1 (MQ=255)
                                                                                       aaTCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATGTaa     >  1:3487254/1‑91 (MQ=255)
                                                                                          |                                                                                          
CAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAATTTAATTA  >  USA300TCH1516_ALE/1044584‑1044764

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: