Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,537,428 A→T 6.5% I94I (ATT→ATA aroB ← 3‑dehydroquinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,537,4280AT 6.5% 622.1 / 23.6 185I94I (ATT→ATAaroB3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base A (89/84);  new base T (9/3);  total (98/87)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 209 total aligned reads displayed.

CGTTGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGCTGG  >  USA300TCH1516_ALE/1537338‑1537517
                                                                                          |                                                                                         
cGTTGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCaata                                                                                           <  2:3308553/91‑1 (MQ=255)
cGTTGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCaata                                                                                           >  2:1551900/1‑91 (MQ=255)
 gTTGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCaataa                                                                                          <  1:500275/91‑1 (MQ=255)
  ttGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCaataat                                                                                         >  1:3478715/1‑91 (MQ=255)
  ttGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCaataat                                                                                         <  2:2349786/91‑1 (MQ=255)
   tGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAattatt                                                                                        <  1:1801942/91‑1 (MQ=255)
    ggCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATTATTg                                                                                       <  2:2455549/91‑1 (MQ=255)
    ggCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTg                                                                                       <  1:1922038/91‑1 (MQ=255)
       aCTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATTATTGCTg                                                                                    <  1:2361596/91‑1 (MQ=255)
        cTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGt                                                                                   <  2:1693244/91‑1 (MQ=255)
           gTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTAtt                                                                                >  1:325995/1‑91 (MQ=255)
           gTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTCGCACCACCACCAACAGCAATAATTGCTGTAtt                                                                                <  2:149036/91‑1 (MQ=255)
             ataAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTAc                                                                              <  2:649619/91‑1 (MQ=255)
                aaGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAg                                                                           >  2:1613101/1‑91 (MQ=255)
                 aGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGt                                                                          >  1:2117162/1‑91 (MQ=255)
                 aGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGt                                                                          >  2:1056424/1‑91 (MQ=255)
                    gAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTAc                                                                       <  1:3742361/91‑1 (MQ=255)
                    gAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTAc                                                                       <  1:1958900/91‑1 (MQ=255)
                       cGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACatg                                                                    >  1:3048768/1‑91 (MQ=255)
                       cGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACatg                                                                    <  1:60848/91‑1 (MQ=255)
                         ccTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACatgat                                                                  <  2:1773725/91‑1 (MQ=255)
                         ccTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACatgat                                                                  >  1:3590941/1‑91 (MQ=255)
                         ccTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACatgat                                                                  >  2:3235972/1‑91 (MQ=255)
                            cGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATggg                                                               >  2:790912/1‑91 (MQ=255)
                             gTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGa                                                              <  1:318647/91‑1 (MQ=255)
                             gTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGa                                                              >  1:511062/1‑91 (MQ=255)
                             gTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGa                                                              <  1:318480/91‑1 (MQ=255)
                              tAAAGGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGAt                                                             <  1:3235972/91‑1 (MQ=255)
                              tAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGAt                                                             >  2:259640/1‑91 (MQ=255)
                               aaaaGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATa                                                            <  1:2474725/91‑1 (MQ=255)
                               aaaaGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATa                                                            >  1:3121234/1‑91 (MQ=255)
                                 aaGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATaaa                                                          >  1:1829879/1‑91 (MQ=255)
                                  aGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATaaaa                                                         >  1:1209954/1‑91 (MQ=255)
                                  aGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATaaaa                                                         >  2:1472763/1‑91 (MQ=255)
                                   gtgtCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGt                                                                          <  2:1337414/73‑1 (MQ=255)
                                   gtgtCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGt                                                                          >  1:1337414/1‑73 (MQ=255)
                                   gtgtCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAat                                                        >  2:3780731/1‑91 (MQ=255)
                                   gtgtCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAat                                                        >  2:2485942/1‑91 (MQ=255)
                                    tgtCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAata                                                       <  1:175887/91‑1 (MQ=255)
                                       cgccgcAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATAc                                                    <  2:1458838/91‑1 (MQ=255)
                                        gccgcAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATAct                                                   <  1:1062562/91‑1 (MQ=255)
                                        gccgcAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATAct                                                   <  1:768102/91‑1 (MQ=255)
                                        gccgcAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATAct                                                   >  1:1116031/1‑91 (MQ=255)
                                         ccgcAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATActc                                                  >  1:1150260/1‑91 (MQ=255)
                                          cgcAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATActct                                                 >  2:430143/1‑91 (MQ=255)
                                           gcAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTa                                                <  2:3121234/91‑1 (MQ=255)
                                            cAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTaa                                               >  2:3105499/1‑91 (MQ=255)
                                                 aaCCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGttt                                          <  2:3630257/91‑1 (MQ=255)
                                                  aCCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATACTTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTc                                         >  2:671591/1‑91 (MQ=255)
                                                  aCCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATAAACTCTAATGTTTc                                         >  1:795429/1‑91 (MQ=255)
                                                   cccAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCt                                        >  1:1867233/1‑91 (MQ=255)
                                                    ccAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCtt                                       <  2:1680011/91‑1 (MQ=255)
                                                    ccAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCtt                                       <  1:3194360/91‑1 (MQ=255)
                                                     cAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTg                                      >  1:640178/1‑91 (MQ=255)
                                                     cAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTg                                      <  1:2216216/91‑1 (MQ=255)
                                                     cAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTg                                      <  1:2213127/91‑1 (MQ=255)
                                                      aGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGa                                     >  1:2527911/1‑91 (MQ=255)
                                                      aGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGa                                     <  2:3874410/91‑1 (MQ=255)
                                                      aGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGa                                     <  1:2082062/91‑1 (MQ=255)
                                                             tCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCt                              >  1:552265/1‑91 (MQ=255)
                                                               aCCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCa                            <  1:1613101/91‑1 (MQ=255)
                                                               aCCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCa                            <  1:2062232/91‑1 (MQ=255)
                                                               aCCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCa                            <  2:2637735/91‑1 (MQ=255)
                                                                ccAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCaa                           >  2:2083429/1‑91 (MQ=255)
                                                                 cAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCaaa                          <  1:888187/91‑1 (MQ=255)
                                                                 cAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCaaa                          <  1:2547450/91‑1 (MQ=255)
                                                                 cAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCaaa                          >  1:1073008/1‑91 (MQ=255)
                                                                    ttGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAACATATACTCTAATGTTTCTTGATATTGCTCAAATGt                       >  2:2965664/1‑91 (MQ=255)
                                                                      gCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACGCTAATGTTTCTTGATATTGCTCAAATGttt                     >  1:2382746/1‑91 (MQ=255)
                                                                       caccaccacCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGtttt                    <  2:3513271/91‑1 (MQ=255)
                                                                       caccaccacCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGtttt                    >  1:678799/1‑91 (MQ=255)
                                                                         ccaccaccaACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCg                  >  1:683186/1‑91 (MQ=255)
                                                                         ccaccaccaACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCg                  >  1:1022794/1‑91 (MQ=255)
                                                                          caccaccaACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGt                 <  2:183723/91‑1 (MQ=255)
                                                                          caccaccaACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGt                 <  1:1505519/91‑1 (MQ=255)
                                                                           accaccaACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGtttt                    >  1:970286/1‑87 (MQ=255)
                                                                            ccaccaACAGCAATAATTGCTGTATTACGAGTTACATGCTGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCt               >  1:172890/1‑91 (MQ=255)
                                                                               ccaACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCtttt            >  1:3631897/1‑91 (MQ=255)
                                                                               ccaACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGATTCTTGATATTGCTCAAATGTTTTCGTCtttt            >  1:3324260/1‑91 (MQ=255)
                                                                                caACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTc           <  2:2184056/91‑1 (MQ=255)
                                                                                 taCAGCAATTATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTCTCGTCTTTTCa          >  1:3524480/2‑91 (MQ=255)
                                                                                 aaCAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCa          <  1:1320537/91‑1 (MQ=255)
                                                                                  aCAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCGAATGTTTTCGTCTTTTCAc         <  1:2273371/91‑1 (MQ=255)
                                                                                  aCAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCAc         <  2:880118/91‑1 (MQ=255)
                                                                                  aCAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCAc         <  1:1472763/91‑1 (MQ=255)
                                                                                  aCAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCAc         <  1:1737143/91‑1 (MQ=255)
                                                                                   cAGCAATTATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCAcc        >  2:3081308/1‑91 (MQ=255)
                                                                                   cAGCAATTATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCAcc        >  1:2853394/1‑91 (MQ=255)
                                                                                   cAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCAcc        >  2:2155565/1‑91 (MQ=255)
                                                                                   cAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCAcc        >  2:3724363/1‑91 (MQ=255)
                                                                                    aGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCa       >  1:2496346/1‑91 (MQ=255)
                                                                                     gCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAAt                                              <  2:1182223/51‑1 (MQ=255)
                                                                                      cAATTATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGc     >  2:2274481/1‑91 (MQ=255)
                                                                                      cAATTATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGc     >  2:2222462/1‑91 (MQ=255)
                                                                                      cAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGc     >  2:2025142/1‑91 (MQ=255)
                                                                                      cAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGc     >  2:534145/1‑91 (MQ=255)
                                                                                      cAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGc     <  1:2452044/91‑1 (MQ=255)
                                                                                        ataatTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGCTg   <  2:1064561/91‑1 (MQ=255)
                                                                                         taatTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGCTgg  <  1:1259071/91‑1 (MQ=255)
                                                                                         taatTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGCTgg  >  2:187252/1‑91 (MQ=255)
                                                                                          |                                                                                         
CGTTGGCACTTGTATAAAGTGAACGCCTCGTAAAAGTGTCGCCGCAATAAACCCAGCAAAATCACCAGTTGCACCACCACCAACAGCAATAATTGCTGTATTACGAGTTACATGATGGGATAAAATATACTCTAATGTTTCTTGATATTGCTCAAATGTTTTCGTCTTTTCACCAGCTGG  >  USA300TCH1516_ALE/1537338‑1537517

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: