Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,413 C→G 5.5% 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,4130CG 5.5% 1605.9 / 61.7 440intergenic (+164/+28)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (199/216);  new base G (13/11);  total (212/228)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 449 total aligned reads displayed.

TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503
                                                                                          |                                                                                          
tatCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATc                                                                                            >  2:2830620/1‑91 (MQ=255)
tatCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATc                                                                                            >  2:1368829/1‑91 (MQ=255)
  tCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAt                                                                                          <  1:1351597/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTgg                                                                                       <  1:1198674/91‑1 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                      <  2:3561499/91‑1 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                      >  1:5720673/1‑91 (MQ=255)
        atgtAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  2:6183604/88‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    >  2:4402603/1‑91 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  1:1729507/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  2:5519014/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  2:3871415/91‑1 (MQ=255)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                   >  1:5110373/1‑91 (MQ=255)
          ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                  <  1:2222228/91‑1 (MQ=255)
          ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                  <  1:3082523/91‑1 (MQ=255)
          ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                  <  2:1961722/91‑1 (MQ=255)
          ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                  <  1:1403773/91‑1 (MQ=255)
            aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGt                                                                                <  1:2685072/91‑1 (MQ=255)
             aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAATACGCCTACCATTCATCTGCTGTCGtt                                                                               >  2:6302731/1‑91 (MQ=255)
             aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtt                                                                               >  2:5451465/1‑91 (MQ=255)
             aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtt                                                                               <  1:4752721/91‑1 (MQ=255)
              aGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGttt                                                                              <  1:2782551/91‑1 (MQ=38)
               gACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttt                                                                             >  2:2092324/1‑91 (MQ=255)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                            <  2:4812841/91‑1 (MQ=255)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                            <  1:6784678/91‑1 (MQ=255)
                   ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                         >  2:762934/1‑91 (MQ=255)
                    tCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgtg                                                                        >  1:1602124/1‑91 (MQ=255)
                     cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                       >  2:2219448/1‑91 (MQ=255)
                     cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                       >  2:3646945/1‑91 (MQ=255)
                     cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                       >  1:4534008/1‑91 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  2:3487173/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  2:3047182/91‑1 (MQ=255)
                       taataTCGCTTTATGTTATTTCGCTTCTTATTTTTCCCTTACTTCAAAGAAAAAAACGACTACCATTCATTGGCAGTCGTTTTTTgtgttt                                                                     <  2:5868899/91‑3 (MQ=255)
                         atatCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTgg                                                                   >  1:5247375/1‑91 (MQ=255)
                         atatCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTgg                                                                   >  1:2609896/1‑91 (MQ=255)
                         atatCGCCTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTgg                                                                   <  1:5573057/91‑1 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                 >  2:5235576/1‑91 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                 >  2:4721033/1‑91 (MQ=255)
                            tCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAtt                                                                >  1:2288895/1‑91 (MQ=255)
                             cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                               <  1:3943259/91‑1 (MQ=255)
                                  tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                          <  2:3744764/91‑1 (MQ=255)
                                    tGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATg                                                        <  2:4998245/91‑1 (MQ=255)
                                     gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGa                                                       <  2:5720673/91‑1 (MQ=255)
                                      ttATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGaa                                                      >  1:4335128/1‑91 (MQ=255)
                                       tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                     <  2:3437030/91‑1 (MQ=255)
                                        aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                                    >  1:3352544/1‑91 (MQ=255)
                                             gCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAatat                                               >  1:5462380/1‑91 (MQ=255)
                                             gCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAatat                                               >  2:458849/1‑91 (MQ=255)
                                               ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                             <  1:6097786/91‑1 (MQ=255)
                                               ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                             <  1:5031565/91‑1 (MQ=255)
                                               ctcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACTAATCATTGTCAGTCGTTTTTTGTTTTTGGATTTAACCCTGAACACATATGa                                             >  2:3204210/2‑91 (MQ=255)
                                                 cttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                           <  1:719412/91‑1 (MQ=255)
                                                  ttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGt                                          >  1:4063627/1‑91 (MQ=255)
                                                   tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                         <  2:6933907/91‑1 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                        <  2:4215425/91‑1 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                        <  2:5362307/91‑1 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                        >  2:6564514/1‑91 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                        >  2:2358747/1‑91 (MQ=255)
                                                     ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACa                                       >  1:1543970/1‑91 (MQ=255)
                                                      tttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTGTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAg                                      <  2:212986/91‑1 (MQ=255)
                                                       ttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAga                                     >  1:2401671/1‑91 (MQ=255)
                                                        tATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaa                                    >  2:6013751/1‑91 (MQ=255)
                                                        tATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaa                                    <  2:367287/91‑1 (MQ=255)
                                                        tATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaa                                    >  2:6742885/1‑91 (MQ=255)
                                                         aTCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaat                                   <  1:3014288/91‑1 (MQ=255)
                                                         aTCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaat                                   <  1:6392293/91‑1 (MQ=255)
                                                          tCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                  <  2:2942920/91‑1 (MQ=255)
                                                           ccTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatga                                 >  2:1113299/1‑91 (MQ=255)
                                                           ccTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatga                                 <  2:414629/91‑1 (MQ=255)
                                                            cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                                <  2:1713597/91‑1 (MQ=255)
                                                             ttACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaat                               >  1:1005035/1‑91 (MQ=255)
                                                             ttACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaat                               <  1:2053529/91‑1 (MQ=255)
                                                                cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                            <  1:3763406/91‑1 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                       <  2:2480369/91‑1 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                       >  1:2444841/1‑91 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAATTACAGAATGAATAGTACCAtt                      <  2:3267142/91‑1 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      >  2:1327978/1‑91 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      >  1:6125627/1‑91 (MQ=255)
                                                                       gAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTa                     >  1:5962435/1‑91 (MQ=255)
                                                                          aaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAata                  >  1:425041/1‑91 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                >  2:6085346/1‑91 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               >  2:1868563/1‑91 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               >  1:486882/1‑91 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               <  1:1868563/91‑1 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              >  2:4541464/1‑91 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              >  1:3534375/1‑91 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             <  2:6987988/91‑1 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             <  2:843861/91‑1 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             <  2:4028893/91‑1 (MQ=255)
                                                                                gACTGCCAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca            >  2:4333940/1‑91 (MQ=38)
                                                                                gACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca            <  2:2709155/91‑1 (MQ=255)
                                                                                 aCTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacat           >  1:58195/1‑91 (MQ=255)
                                                                                 aCTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacat           >  1:4614392/1‑91 (MQ=255)
                                                                                  cTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  1:6239690/1‑91 (MQ=255)
                                                                                   tgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa         >  1:289143/3‑91 (MQ=255)
                                                                                   tgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa         >  1:1071731/3‑91 (MQ=39)
                                                                                   tACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa         <  1:802422/91‑1 (MQ=255)
                                                                                     ccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATc       >  1:4170069/1‑91 (MQ=255)
                                                                                        aTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGa    <  2:4666472/91‑1 (MQ=255)
                                                                                        aTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGa    <  1:3239079/91‑1 (MQ=255)
                                                                                         tCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAt   >  1:5286629/1‑91 (MQ=255)
                                                                                          |                                                                                          
TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503

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