Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,419 C→T 5.3% intergenic (+170/+22) 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,4190CT 5.3% 1507.7 / 47.7 442intergenic (+170/+22)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (203/215);  new base T (20/4);  total (223/219)
Fisher's exact test for biased strand distribution p-value = 1.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 454 total aligned reads displayed.

ATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTC  >  USA300TCH1516_ALE/2372329‑2372509
                                                                                          |                                                                                           
atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                             >  1:5277138/1‑91 (MQ=255)
  atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                           <  1:5210732/91‑1 (MQ=255)
  atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                           <  1:1955740/91‑1 (MQ=255)
   tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                          >  1:5110373/1‑91 (MQ=255)
    ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                         <  1:3908608/91‑1 (MQ=255)
    ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                         <  2:4006149/91‑1 (MQ=255)
    ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                         <  1:5883963/91‑1 (MQ=255)
     tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                        <  2:6086147/91‑1 (MQ=38)
     tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTCATTTTATCCTTGCTTCAAAGAAAAAAACGACTACCGATCATTGTCAGTCg                                                                                        >  1:5811408/1‑91 (MQ=255)
      aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑t                                                                                      <  1:459450/91‑1 (MQ=255)
       aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑tt                                                                                     >  2:5450619/1‑91 (MQ=255)
        aGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑ttt                                                                                    >  2:3510909/1‑91 (MQ=255)
         gACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑tttt                                                                                   <  1:5434057/91‑1 (MQ=255)
          aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑ttttt                                                                                  <  2:4812841/91‑1 (MQ=255)
          aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑ttttt                                                                                  <  1:6784678/91‑1 (MQ=255)
                ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTc                                                                            <  2:4353633/91‑1 (MQ=255)
                ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTc                                                                            <  2:29761/91‑1 (MQ=255)
                ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTc                                                                            <  2:188442/91‑1 (MQ=255)
                 taataTCGCTTTATGTTATTTCGCTTCTTATTTTTCCCTTACTTCAAAGAAAAAAACGACTACCATTCATTGGCAGTCG‑TTTTTTgtgttt                                                                           <  2:5868899/91‑3 (MQ=255)
                 taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCt                                                                           <  1:3250100/91‑1 (MQ=255)
                  aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCTATCATTGGCTGTCG‑TTTTTtgtgtgtg                                                                          >  2:6840454/1‑91 (MQ=255)
                     atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGAt                                                                       >  1:2619418/1‑91 (MQ=255)
                     atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGAt                                                                       >  2:4721033/1‑91 (MQ=255)
                      tCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGAtt                                                                      >  1:2288895/1‑91 (MQ=255)
                       cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGAttt                                                                     >  1:1590508/1‑91 (MQ=255)
                          tttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAAc                                                                  <  2:5332246/91‑1 (MQ=255)
                           ttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAAcc                                                                 <  2:2002716/91‑1 (MQ=255)
                              tGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTCGATTTAACCATg                                                              >  2:4469835/1‑91 (MQ=255)
                               gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGa                                                             <  2:5720673/91‑1 (MQ=255)
                                ttATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGaa                                                            <  2:2209968/91‑1 (MQ=255)
                                 tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAAc                                                           <  2:5478784/91‑1 (MQ=255)
                                  aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACa                                                          >  2:4978345/1‑91 (MQ=255)
                                  aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACa                                                          >  1:3352544/1‑91 (MQ=255)
                                  aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACa                                                          <  1:6222332/91‑1 (MQ=255)
                                      cGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAata                                                      <  2:3456881/91‑1 (MQ=255)
                                       gCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAatat                                                     >  2:458849/1‑91 (MQ=255)
                                         ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                                   >  2:2096426/1‑91 (MQ=255)
                                         ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                                   <  2:6203418/91‑1 (MQ=255)
                                           cttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                                 <  1:719412/91‑1 (MQ=255)
                                             tagtttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                               <  1:5692714/88‑1 (MQ=255)
                                             tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                               >  1:4333932/1‑91 (MQ=255)
                                             tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                               <  2:5163328/91‑1 (MQ=255)
                                             tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                               <  2:6933907/91‑1 (MQ=255)
                                             tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                               >  2:3372621/1‑91 (MQ=255)
                                              aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                              <  2:4215425/91‑1 (MQ=255)
                                               ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACa                                             >  1:4229731/1‑91 (MQ=255)
                                                 ttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAga                                           <  1:4124252/91‑1 (MQ=255)
                                                  tATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaa                                          >  2:6013751/1‑91 (MQ=255)
                                                    tCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                        >  2:2869911/1‑91 (MQ=255)
                                                    tCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                        >  1:4028893/1‑91 (MQ=255)
                                                      cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                                      <  1:5567676/91‑1 (MQ=255)
                                                      cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                                      >  1:3107585/1‑91 (MQ=255)
                                                      cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                                      <  1:3621692/91‑1 (MQ=255)
                                                      cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                                      <  1:4402603/91‑1 (MQ=255)
                                                       ttACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaat                                     >  1:1005035/1‑91 (MQ=255)
                                                       ttACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaat                                     <  1:2053529/91‑1 (MQ=255)
                                                        tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                                    >  1:1205213/1‑91 (MQ=255)
                                                        tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                                    >  1:6107910/1‑91 (MQ=255)
                                                        tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                                    <  2:3231293/91‑1 (MQ=255)
                                                         aCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAg                                   <  1:3036465/91‑1 (MQ=255)
                                                         aCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAg                                   >  2:366735/1‑91 (MQ=255)
                                                         aCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAg                                   <  1:1770624/91‑1 (MQ=255)
                                                          cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                                  <  1:6000651/91‑1 (MQ=255)
                                                            tCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAc                                >  1:2752304/1‑91 (MQ=255)
                                                             cAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                               <  1:2904583/91‑1 (MQ=255)
                                                             cAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                               >  1:24052/1‑91 (MQ=255)
                                                              aaaGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCa                              >  2:718151/1‑91 (MQ=255)
                                                               aaGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                             <  2:6603507/91‑1 (MQ=255)
                                                               aaGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                             >  2:1787570/1‑91 (MQ=255)
                                                               aaGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                             <  2:2480369/91‑1 (MQ=255)
                                                                aGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                            >  2:1327978/1‑91 (MQ=255)
                                                                    aaaaaCGACTGCCAATGATTGGTAGTCGTTTTTTTGTGTCTGGATTTAACCATGAACAAATATTAAGTACAGAATGAATAGTACCATTAat                         >  2:3630507/1‑91 (MQ=25)
                                                                     aaaaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatac                       <  2:3750922/91‑1 (MQ=255)
                                                                     aaaaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatac                       <  1:6391252/91‑1 (MQ=255)
                                                                      aaaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                      <  2:234450/91‑1 (MQ=255)
                                                                       aaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat                     >  2:4121357/1‑91 (MQ=255)
                                                                       aaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat                     <  2:3833002/91‑1 (MQ=255)
                                                                       aaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat                     >  1:486882/1‑91 (MQ=255)
                                                                       aaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat                     >  1:4944399/1‑91 (MQ=255)
                                                                       aaCGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat                     >  2:3138593/1‑91 (MQ=255)
                                                                          gACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca                  >  1:2474635/1‑91 (MQ=255)
                                                                           aCTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacat                 <  2:83012/91‑1 (MQ=255)
                                                                            ctgccAATGATTGGTAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt                >  1:1367519/4‑91 (MQ=39)
                                                                             tgccAATGATTGGTAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa               >  1:1071731/3‑91 (MQ=39)
                                                                             tACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa               <  1:802422/91‑1 (MQ=255)
                                                                              aCCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt              <  1:2225876/91‑1 (MQ=255)
                                                                               ccAATGATTGGTAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATc             >  1:4715231/1‑91 (MQ=255)
                                                                                 aatgaTTGGTAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg           >  1:1682428/5‑91 (MQ=255)
                                                                                 aaTCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGATCAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg           <  2:6983292/91‑1 (MQ=255)
                                                                                  aTCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGa          <  2:4666472/91‑1 (MQ=255)
                                                                                   tgaTTGGTAGTCG‑TTTTTTGTGTCTGGAATTAAACATGAACAAATATGAAGTAAAGAATGAATAGTACCATTAATAAATACATTATCGGaa         >  2:903641/3‑90 (MQ=255)
                                                                                   tCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAt         >  2:2659098/1‑91 (MQ=255)
                                                                                    cATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATg        <  2:1570087/91‑1 (MQ=255)
                                                                                      ttGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTa      >  1:1450830/1‑91 (MQ=255)
                                                                                        ggCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTAct    >  2:4314360/1‑91 (MQ=255)
                                                                                        ggCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTAct    >  2:1083769/1‑91 (MQ=255)
                                                                                          cAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTActtc  >  2:5610611/1‑91 (MQ=255)
                                                                                          cAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTActtc  >  2:6311034/1‑91 (MQ=255)
                                                                                          cAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTActtc  >  1:3420720/1‑91 (MQ=255)
                                                                                          cAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTActtc  >  1:27875/1‑91 (MQ=255)
                                                                                          |                                                                                           
ATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCG‑TTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTC  >  USA300TCH1516_ALE/2372329‑2372509

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: ‑