Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,413 C→G 20.7% intergenic (+164/+28) USA300TCH1516_02262 → / ← pbuG hypothetical protein/Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,372,4130CG20.7% 304.6 / 74.0 125intergenic (+164/+28)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (44/55);  new base G (7/19);  total (51/74)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 127 total aligned reads displayed.

TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAT  >  USA300TCH1516_ALE/2372323‑2372502
                                                                                          |                                                                                         
tatCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATc                                                                                           <  1:952821/91‑1 (MQ=255)
tatCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATc                                                                                           >  1:199362/1‑91 (MQ=255)
  tCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAt                                                                                         >  1:759573/1‑91 (MQ=255)
   cGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGAtt                                                                                        <  1:2027104/91‑1 (MQ=255)
   cGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGAtt                                                                                        <  1:2027239/91‑1 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTg                                                                                       <  1:2531411/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTgg                                                                                      <  1:918686/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTgg                                                                                      <  1:1471665/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTgg                                                                                      >  1:1761336/1‑91 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggt                                                                                     <  1:1727902/91‑2 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                     <  1:261891/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                   <  1:3191991/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                   <  1:2477427/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                   <  1:1240006/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                   <  1:4025121/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTCACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                   <  1:2392525/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTAAGTTATTTCGCTTCGGATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                   <  1:1240024/91‑1 (MQ=255)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGt                                                                                  <  1:3319852/91‑1 (MQ=38)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                  >  1:2661097/1‑91 (MQ=255)
          ttAAAGACATTCATAATATCGCTTTCTGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTc                                                                                 <  1:3856226/91‑1 (MQ=38)
          ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTc                                                                                 <  1:3132479/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTATTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                <  1:1896286/91‑1 (MQ=37)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                <  1:2023298/91‑1 (MQ=255)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                <  1:938857/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                <  1:733610/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCg                                                                                >  1:1878631/1‑91 (MQ=255)
            aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGt                                                                               <  1:889931/91‑1 (MQ=255)
             aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGtt                                                                              <  1:2303269/91‑1 (MQ=38)
             aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtt                                                                              >  1:2996340/1‑91 (MQ=255)
                  aTTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtg                                                                         <  1:1392129/91‑1 (MQ=255)
                    tCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgtg                                                                       <  1:2787899/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                     <  1:1076575/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                     <  1:4047671/91‑1 (MQ=255)
                       taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                    <  1:824504/91‑1 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                >  1:2290965/1‑91 (MQ=255)
                              gCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTa                                                             >  1:2688217/1‑91 (MQ=255)
                              gCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTa                                                             >  1:6290/1‑91 (MQ=255)
                                  tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                         <  1:1522283/91‑1 (MQ=255)
                                  tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATGTAACCa                                                         <  1:3229399/91‑1 (MQ=255)
                                  tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATCTAACCa                                                         <  1:3226961/91‑1 (MQ=255)
                                    tGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATg                                                       <  1:1932151/91‑1 (MQ=255)
                                      ttATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGaa                                                     >  1:3789576/1‑91 (MQ=255)
                                      ttATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGaa                                                     <  1:4005443/91‑1 (MQ=255)
                                       tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                    <  1:935656/91‑1 (MQ=255)
                                       tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                    <  1:2655584/91‑1 (MQ=255)
                                           tCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAat                                                <  1:2550397/91‑1 (MQ=255)
                                           tCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAat                                                <  1:2550411/91‑1 (MQ=255)
                                           tCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAat                                                <  1:518888/91‑1 (MQ=255)
                                            cGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAata                                               >  1:1675351/1‑91 (MQ=255)
                                              cttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATCGGCAGTCGTTTTTTGTGTCTGGATTTAACGATGAACAAATATg                                             >  1:4148062/1‑91 (MQ=255)
                                               ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                            <  1:2866640/91‑1 (MQ=255)
                                               ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                            <  1:1314409/91‑1 (MQ=255)
                                                 cttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                          >  1:710425/1‑91 (MQ=255)
                                                 cttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                          >  1:710282/1‑91 (MQ=255)
                                                  ttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGt                                         <  1:2233898/91‑1 (MQ=255)
                                                  ttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGt                                         >  1:2635301/1‑91 (MQ=255)
                                                  ttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGt                                         <  1:2233897/91‑1 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                       <  1:3724782/91‑1 (MQ=255)
                                                     ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACa                                      <  1:2397149/91‑1 (MQ=255)
                                                     ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACa                                      <  1:2397212/91‑1 (MQ=255)
                                                      tttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAg                                     >  1:1206997/1‑91 (MQ=255)
                                                        tATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaa                                   <  1:1521357/91‑1 (MQ=255)
                                                            cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                               <  1:3575524/91‑1 (MQ=255)
                                                               aCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAg                            >  1:328712/1‑91 (MQ=255)
                                                                cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                           >  1:628978/1‑91 (MQ=255)
                                                                 ttCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCCTGCACACATATGAAGTACAGAATGGATAGTa                          >  1:1363285/1‑91 (MQ=255)
                                                                 ttCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTa                          <  1:1827294/91‑1 (MQ=255)
                                                                  tCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAc                         >  1:113880/1‑91 (MQ=255)
                                                                   cAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                        <  1:2266807/91‑1 (MQ=255)
                                                                    aaaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCa                       <  1:1117850/91‑1 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                      >  1:1756436/1‑91 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                      >  1:945867/1‑91 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                      >  1:2604734/1‑91 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                      >  1:2604750/1‑91 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                      >  1:2663785/1‑91 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                      >  1:3951819/1‑91 (MQ=255)
                                                                          aaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAata                 >  1:467333/1‑91 (MQ=255)
                                                                          aaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAata                 >  1:3350130/1‑91 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGGACCATTAataca               <  1:3194797/91‑1 (MQ=255)
                                                                             aaCGACTGCCAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat              >  1:1379748/1‑91 (MQ=38)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat              <  1:2054719/91‑1 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata             >  1:2185444/1‑91 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata             >  1:618390/1‑91 (MQ=255)
                                                                               cGACTGCCAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac            >  1:3622879/1‑91 (MQ=38)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac            >  1:2091535/1‑91 (MQ=255)
                                                                                gACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca           <  1:2763801/91‑1 (MQ=255)
                                                                                   tACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa        <  1:2366048/91‑1 (MQ=255)
                                                                                   tACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa        <  1:2366033/91‑1 (MQ=255)
                                                                                    gccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt       >  1:3079110/2‑91 (MQ=255)
                                                                                    aCCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt       <  1:1054613/91‑1 (MQ=255)
                                                                                     ccAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTGACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATc      >  1:2960519/1‑91 (MQ=255)
                                                                                      cAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg     >  1:4026647/1‑91 (MQ=255)
                                                                                      cAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg     >  1:2753452/1‑91 (MQ=255)
                                                                                      cAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg     >  1:3886619/1‑91 (MQ=255)
                                                                                      cAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg     <  1:622351/91‑1 (MQ=255)
                                                                                      cAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg     <  1:1222867/91‑1 (MQ=255)
                                                                                       aaTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg    >  1:2871700/1‑91 (MQ=255)
                                                                                       aaTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg    >  1:2870466/1‑91 (MQ=255)
                                                                                        aTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGa   <  1:3419610/91‑1 (MQ=255)
                                                                                         tCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAt  >  1:162659/1‑91 (MQ=255)
                                                                                          |                                                                                         
TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAT  >  USA300TCH1516_ALE/2372323‑2372502

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: