Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 536,876 T→C 5.2% intergenic (+69/‑242) rplY → / → pth 50S ribosomal protein L25/Peptidyl‑tRNA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE536,8760TC 5.2% 735.1 / 19.7 228intergenic (+69/‑242)rplY/pth50S ribosomal protein L25/Peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (114/102);  new base C (9/3);  total (123/105)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 235 total aligned reads displayed.

CGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACT  >  USA300TCH1516_ALE/536786‑536966
                                                                                          |                                                                                          
cGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCAt                                                                                            >  1:242212/1‑91 (MQ=255)
   agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGt                                                                                         >  1:160435/1‑91 (MQ=255)
    gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                                        <  1:471207/91‑1 (MQ=255)
     aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                                       >  2:955786/1‑91 (MQ=255)
     aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                                       <  2:3703195/91‑1 (MQ=255)
      aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                                      <  1:1892314/91‑1 (MQ=255)
      aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                                      <  1:3513473/91‑1 (MQ=255)
       aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                                     >  2:2065200/1‑91 (MQ=255)
        aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                                    <  1:2948012/91‑1 (MQ=255)
        aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTGAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                                    <  2:3186025/91‑1 (MQ=255)
          aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                                  >  1:1641654/1‑91 (MQ=255)
          aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                                  >  2:1648984/1‑91 (MQ=255)
          aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                                  <  2:1550054/91‑1 (MQ=255)
           cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCACAGTGCttgtta                                                                                 <  2:3751939/91‑5 (MQ=18)
            tGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgt                                                                                >  2:2018183/1‑91 (MQ=255)
             gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                               >  2:256031/1‑91 (MQ=255)
              aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                              >  1:932267/1‑91 (MQ=255)
              aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                              >  1:2197994/1‑91 (MQ=255)
               agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                             >  2:1457000/1‑91 (MQ=255)
               agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                             >  2:707323/1‑91 (MQ=255)
               agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                             <  2:2244191/91‑1 (MQ=255)
                 agTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttat                                                                           <  2:1225903/91‑1 (MQ=255)
                  gTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                          >  1:638671/1‑91 (MQ=255)
                  gTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAACCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                          >  1:3078602/1‑91 (MQ=255)
                     aTTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATtata                                                                       >  1:3749180/1‑91 (MQ=255)
                       tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                                     >  2:2594824/1‑91 (MQ=255)
                       tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                                     <  1:1085037/91‑1 (MQ=255)
                       tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                                     >  1:2868942/1‑91 (MQ=255)
                         tAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                                   <  1:1551907/91‑1 (MQ=255)
                          aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                          >  2:2628867/1‑83 (MQ=255)
                          aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                          <  1:2628867/83‑1 (MQ=255)
                          aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                                  >  2:1435378/1‑91 (MQ=255)
                          aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                                  <  2:3538774/91‑1 (MQ=255)
                          aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                                  <  1:2850635/91‑1 (MQ=255)
                           aTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtt                                                                 <  2:1217954/91‑1 (MQ=255)
                            tCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtta                                                                >  2:1590832/1‑91 (MQ=255)
                              tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                              >  1:1917202/1‑91 (MQ=255)
                               gTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGTttaatt                                                             <  2:3800139/91‑1 (MQ=255)
                                ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                                            >  1:2515765/1‑91 (MQ=255)
                                 tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                           >  2:648975/1‑91 (MQ=255)
                                      tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                                      <  1:460211/91‑1 (MQ=255)
                                      tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                                      >  1:3141576/1‑91 (MQ=255)
                                      tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTCAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                                      >  1:3360596/1‑91 (MQ=255)
                                       aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                                     <  2:2285346/91‑1 (MQ=255)
                                        aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                                    >  2:405459/1‑91 (MQ=255)
                                        aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                                    <  2:1266503/91‑1 (MQ=255)
                                        aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                                    >  2:3365603/1‑91 (MQ=255)
                                          gTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatta                                                  <  1:1050197/91‑1 (MQ=255)
                                           tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                                           >  1:1195792/1‑49 (MQ=255)
                                            ttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACt                                                >  1:2532138/1‑91 (MQ=255)
                                            ttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACt                                                >  1:1907928/1‑91 (MQ=255)
                                             tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                               >  2:2863606/1‑91 (MQ=255)
                                             tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                               >  2:2863888/1‑91 (MQ=255)
                                             tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                               <  1:955786/91‑1 (MQ=255)
                                             tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                               <  2:3209058/91‑1 (MQ=255)
                                               tataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTg                                             <  1:202338/91‑1 (MQ=255)
                                               tataCTTTGTTTAACAAGCACTATGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTg                                             <  2:3121959/91‑1 (MQ=255)
                                                ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                                            <  2:2038226/91‑1 (MQ=255)
                                                ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                                            <  2:1731610/91‑1 (MQ=255)
                                                 taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                                           <  1:3681623/91‑1 (MQ=255)
                                                   cTTTGTTTAACGAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGGTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                         <  2:3623763/91‑1 (MQ=255)
                                                   cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                         >  1:1453327/1‑91 (MQ=255)
                                                   cTTTGGTTAACAAGCACTGTGCTTATGTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                         <  1:3766062/91‑1 (MQ=255)
                                                    tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                                        <  2:2671489/91‑1 (MQ=255)
                                                    tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                                        >  1:3802781/1‑91 (MQ=255)
                                                      tGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaaa                                      >  1:3250442/1‑91 (MQ=255)
                                                        tttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGt                                    >  2:3402720/1‑91 (MQ=255)
                                                          tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTAGTACAAAGTAATAAAGttt                                  <  2:2763081/91‑1 (MQ=25)
                                                            aaaaGCACCATGCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                                >  2:795440/3‑91 (MQ=21)
                                                            aCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                                <  1:1432297/91‑1 (MQ=255)
                                                            aCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                                >  1:2170579/1‑91 (MQ=255)
                                                             cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                               <  1:660099/91‑1 (MQ=255)
                                                               aGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatta                             >  2:2550121/1‑91 (MQ=255)
                                                                 cacagtgCTTGTTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat                           >  1:1678178/5‑91 (MQ=255)
                                                                 cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTACAGTTtaattaat                           <  1:2986159/91‑1 (MQ=255)
                                                                 cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat                           <  1:1826906/91‑1 (MQ=255)
                                                                  accatgCTTATATTAAAATAAGCACAGTGCTTTTTGCGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt                          >  2:237221/5‑91 (MQ=25)
                                                                   ccatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt                         >  1:1890819/4‑91 (MQ=37)
                                                                   cTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt                         >  1:2076035/1‑91 (MQ=255)
                                                                   cTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTAGTACAAAGTAATAAAGTttaattaagct                         <  1:3843260/91‑4 (MQ=18)
                                                                    tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt                        >  2:243229/1‑91 (MQ=255)
                                                                     gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa                       >  1:2693309/1‑91 (MQ=255)
                                                                      tgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg                      >  1:901108/1‑91 (MQ=255)
                                                                      tgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg                      >  1:3545520/1‑91 (MQ=255)
                                                                       gCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGt                     >  1:2647/1‑91 (MQ=255)
                                                                        cTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTg                    >  2:1973834/1‑91 (MQ=255)
                                                                         ttATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGa                   <  2:2868749/91‑1 (MQ=255)
                                                                          tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg                  >  1:2583562/1‑91 (MQ=255)
                                                                          tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg                  >  1:2761704/1‑91 (MQ=255)
                                                                           aTTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGt                 <  1:775002/91‑1 (MQ=255)
                                                                             attAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaa               >  1:1066151/2‑91 (MQ=255)
                                                                              ttAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaa              <  2:2122823/91‑1 (MQ=255)
                                                                               tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa             >  1:3691599/1‑91 (MQ=255)
                                                                               tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa             >  2:494394/1‑91 (MQ=255)
                                                                                    taAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACAt        <  1:3684674/91‑1 (MQ=255)
                                                                                    taAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTCATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACAt        >  2:1095018/1‑91 (MQ=255)
                                                                                     aaGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACAtt       <  2:1037750/91‑1 (MQ=255)
                                                                                        cATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTaaa    <  2:3504978/91‑1 (MQ=255)
                                                                                         aTGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAAc   >  2:2208742/1‑91 (MQ=255)
                                                                                          tGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACt  >  2:1026696/1‑91 (MQ=255)
                                                                                          |                                                                                          
CGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACT  >  USA300TCH1516_ALE/536786‑536966

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: