Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,413 C→G 22.1% 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,4130CG22.1% 625.4 / 174.3 271intergenic (+164/+28)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (104/106);  new base G (19/41);  total (124/147)
Fisher's exact test for biased strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

Only 100 of 276 total aligned reads displayed.

TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503
                                                                                          |                                                                                          
  tCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAt                                                                                          <  1:1106435/91‑1 (MQ=255)
   cGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAtt                                                                                         <  2:468183/91‑1 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTg                                                                                        <  1:795974/91‑1 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTg                                                                                        <  2:2528914/91‑1 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTg                                                                                        >  2:406834/1‑91 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTg                                                                                        >  1:2501657/1‑91 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTgg                                                                                       <  1:818151/91‑1 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggt                                                                                      <  2:3626348/91‑2 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggt                                                                                      <  1:2238519/91‑2 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                      <  1:2774719/91‑1 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                      <  1:3820410/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  1:2076921/91‑1 (MQ=255)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                   >  1:1541196/1‑91 (MQ=255)
          ttAATGACAAACATAATATCGCTTTATGTTATTTCGCTTCGTATATTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                  <  2:1279032/91‑1 (MQ=255)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  2:2994197/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  2:1747190/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  1:978500/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  1:3315182/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  1:495917/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  2:2475670/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  2:3025415/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCACTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  1:572904/91‑1 (MQ=37)
            aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAGGAAAAAAACGACTGCCAATCATTGGCAGTCGt                                                                                <  2:2234557/91‑1 (MQ=255)
            aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGt                                                                                <  1:2924929/91‑1 (MQ=38)
             aaGACATTCATAATATCGCTTTATGTTATTTCGCTTTTTATTTTATCCTTACTTCAAAGAAAAAAACGTCTACCAATCTTTTGCAGTCGtt                                                                               >  1:121671/1‑91 (MQ=255)
             aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGtt                                                                               <  1:2561409/91‑1 (MQ=38)
               gACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTGCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttt                                                                             <  2:1934511/91‑1 (MQ=255)
               gACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttt                                                                             <  2:2097539/91‑1 (MQ=255)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGttttt                                                                            <  1:3631853/91‑1 (MQ=37)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                            <  1:77260/91‑1 (MQ=255)
                 cATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttttt                                                                           >  1:1906816/1‑91 (MQ=255)
                 cATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttttt                                                                           >  1:410449/1‑91 (MQ=255)
                   ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                         >  1:1274173/1‑91 (MQ=255)
                   ttCATAATATCGCTTTATGTTATATCGCTTCTTATTTGATCCTTACCTCAAAGAAAAAACCGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                         >  1:3154576/1‑91 (MQ=255)
                     cATTATATCGGTTTATGGTGTTTGGCTTCTTATTTTAGCGTTCCTTCAAAGAAAAAAACAACAACCAATGCGTGGTAGTCGTTTTTTgtgt                                                                       >  2:3735787/1‑91 (MQ=255)
                     cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                       >  2:277571/1‑91 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      >  2:1991198/1‑91 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  1:860756/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  1:521121/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  1:2265984/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  1:2292741/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  2:3653670/91‑1 (MQ=255)
                       taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                     <  2:670581/91‑1 (MQ=255)
                        aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                                    >  1:914689/1‑91 (MQ=255)
                        aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                                    <  1:2187025/91‑1 (MQ=255)
                          tatCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGa                                                                  <  1:3357322/91‑1 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                 <  1:2344050/91‑1 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                 >  2:2977479/1‑91 (MQ=255)
                             cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTCTTTTGTGTCTGGAttt                                                               >  2:3790456/1‑91 (MQ=255)
                              gCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTAGATTTa                                                              <  1:3080757/91‑1 (MQ=255)
                               cTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTaa                                                             >  1:3934951/1‑91 (MQ=255)
                                 ttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAcc                                                           <  2:3442025/91‑1 (MQ=255)
                                    tGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATg                                                        >  2:2342871/1‑91 (MQ=255)
                                     gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGa                                                       <  1:3884151/91‑1 (MQ=255)
                                       tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                     <  1:2008019/91‑1 (MQ=255)
                                        aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                                    >  1:1036596/1‑91 (MQ=255)
                                        aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                                    >  2:1194641/1‑91 (MQ=255)
                                          ttCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACaaa                                                  >  1:606750/1‑91 (MQ=255)
                                           tCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAat                                                 >  1:3000488/1‑91 (MQ=255)
                                              cttcttATTTTATCCATACTTTAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTCTTTTGTGTCTGGATATAACCATCACGAAATATg                                              >  2:1359359/1‑91 (MQ=255)
                                               ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                             <  2:1592732/91‑1 (MQ=255)
                                                  ttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATTATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGt                                          >  2:718123/1‑91 (MQ=255)
                                                   tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                         <  1:3501837/91‑1 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                        >  2:2834334/1‑91 (MQ=255)
                                                      tttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAg                                      >  2:3220928/1‑91 (MQ=255)
                                                         aTCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaat                                   <  2:3791637/91‑1 (MQ=255)
                                                          tCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                  <  2:54721/91‑1 (MQ=255)
                                                          tCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                  >  1:213390/1‑91 (MQ=255)
                                                               aCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAg                             <  1:1150042/91‑1 (MQ=255)
                                                                cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                            <  1:1752035/91‑1 (MQ=255)
                                                                    aaaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCa                        <  1:3232401/91‑1 (MQ=255)
                                                                    aaaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCa                        <  1:3811819/91‑1 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACACATATGAAGTACAGAATGAATAGTACCac                       >  2:1040257/1‑90 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      >  1:2155986/1‑91 (MQ=255)
                                                                        aaaaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTaa                    >  2:186036/1‑91 (MQ=255)
                                                                         aaaaaaCGACTACCAATCATTGGCAGTCGTTTTTTTTGCCTGGATTTAACCATGAACCCATATGAAGTACAGAATGAATAGTACTATTAat                   >  2:797500/1‑91 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                >  1:3276268/1‑91 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               >  2:961944/1‑91 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTGTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               <  2:2795183/91‑1 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              <  1:2679098/91‑1 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              >  2:2120116/1‑91 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             <  1:3442766/91‑1 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             >  2:3339308/1‑91 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             <  1:3449696/91‑1 (MQ=255)
                                                                                 actgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacat           >  2:188983/5‑91 (MQ=38)
                                                                                  ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  1:584330/4‑91 (MQ=39)
                                                                                  ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  1:2459081/4‑91 (MQ=39)
                                                                                  ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  2:3928282/4‑91 (MQ=39)
                                                                                  ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  2:2697119/4‑91 (MQ=39)
                                                                                  cTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  2:2363627/1‑91 (MQ=255)
                                                                                    gccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGCATAGTACCATTAATACATACATTAt        >  2:2588186/2‑91 (MQ=39)
                                                                                    aCCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt        <  1:1450956/91‑1 (MQ=255)
                                                                                    aCCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt        >  1:1788762/1‑91 (MQ=255)
                                                                                     ccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATc       >  2:1379510/1‑91 (MQ=255)
                                                                                     ccAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATc       >  1:758573/1‑91 (MQ=255)
                                                                                      cAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg      <  1:1295776/91‑1 (MQ=255)
                                                                                       aaTGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGCATGAATAGTACCATTAATACATACATTATCgg     >  2:3876120/1‑91 (MQ=255)
                                                                                       aaTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAACTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg     >  2:89442/1‑91 (MQ=255)
                                                                                         tCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAt   <  2:3959128/91‑1 (MQ=255)
                                                                                          cATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATg  <  1:3700157/91‑1 (MQ=255)
                                                                                          |                                                                                          
TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503

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