Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 536,877 G→A 5.2% intergenic (+70/‑241) 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,8770GA 5.2% 755.8 / 17.5 230intergenic (+70/‑241)rplY/pth50S ribosomal protein L25/Peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base G (116/101);  new base A (9/3);  total (126/104)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 237 total aligned reads displayed.

GAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACTC  >  USA300TCH1516_ALE/536787‑536967
                                                                                          |                                                                                          
  agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGt                                                                                          >  1:3223153/1‑91 (MQ=255)
  agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGt                                                                                          >  1:160435/1‑91 (MQ=255)
    aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                                        >  2:955786/1‑91 (MQ=255)
    aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                                        <  2:3703195/91‑1 (MQ=255)
     aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                                       <  1:1892314/91‑1 (MQ=255)
      aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                                      >  1:917327/1‑91 (MQ=255)
      aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                                      >  2:2035301/1‑91 (MQ=255)
      aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCACAGTGCtt                                                                                      <  2:2204140/91‑1 (MQ=21)
       aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                                     >  2:1583149/1‑91 (MQ=255)
       aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTGAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                                     <  2:3186025/91‑1 (MQ=255)
        aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                                    >  1:762854/1‑91 (MQ=255)
         aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                                   >  1:1641654/1‑91 (MQ=255)
         aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                                   <  2:1032139/91‑1 (MQ=255)
          cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                                  <  1:2684606/91‑1 (MQ=255)
          cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCACAGTGCttgtta                                                                                  <  2:3751939/91‑5 (MQ=18)
           tGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgt                                                                                 >  2:2018183/1‑91 (MQ=255)
           tGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATCGTGCTTTTtgt                                                                                 >  1:3518111/1‑91 (MQ=255)
              agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                              >  2:707323/1‑91 (MQ=255)
                  tAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatta                                                                          <  1:3557345/91‑1 (MQ=255)
                    aTTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATtata                                                                        >  1:3749902/1‑91 (MQ=255)
                    aTTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTCTTGTGTTATtata                                                                        >  2:2553830/1‑91 (MQ=255)
                      tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                                      >  2:2594824/1‑91 (MQ=255)
                      tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                                      <  1:1085037/91‑1 (MQ=255)
                      tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                                      >  1:3704909/1‑91 (MQ=255)
                        tAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                                    <  1:1551907/91‑1 (MQ=255)
                        tAATCTGTTACATTAAAGTTTTTATACTTTGTGTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                                    <  1:3313872/91‑1 (MQ=255)
                         aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                           <  1:2628867/83‑1 (MQ=255)
                         aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                           >  2:2628867/1‑83 (MQ=255)
                              gTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGTttaatt                                                              <  2:3800139/91‑1 (MQ=255)
                               ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                                             >  1:2515765/1‑91 (MQ=255)
                                tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                            >  2:1223387/1‑91 (MQ=255)
                                tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                            >  2:648975/1‑91 (MQ=255)
                                      aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                                      <  2:2285346/91‑1 (MQ=255)
                                       aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                                     <  2:1533933/91‑1 (MQ=255)
                                       aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                                     <  2:1266503/91‑1 (MQ=255)
                                       aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                                     >  2:2648790/1‑91 (MQ=255)
                                       aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                                     >  2:405459/1‑91 (MQ=255)
                                       aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttag                                                     <  1:953403/91‑2 (MQ=255)
                                        aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatt                                                    <  1:2872594/91‑1 (MQ=255)
                                        aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatt                                                    <  1:285555/91‑1 (MQ=255)
                                        aGTTTTTATACTGTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatt                                                    >  1:2973833/1‑91 (MQ=255)
                                         gTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatta                                                   <  1:1050197/91‑1 (MQ=255)
                                          tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                                            <  2:1195792/49‑1 (MQ=255)
                                          tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                                            >  1:1195792/1‑49 (MQ=255)
                                           ttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACt                                                 >  1:1907928/1‑91 (MQ=255)
                                            tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                                <  2:3209058/91‑1 (MQ=255)
                                            tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                                >  2:2863887/1‑91 (MQ=255)
                                            tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                                <  2:2840213/91‑1 (MQ=255)
                                              tataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTg                                              <  1:202338/91‑1 (MQ=255)
                                               ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAGAGCTTAATTAAACTTTATTACTTTGt                                             <  2:2218806/91‑1 (MQ=255)
                                               ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                                             <  2:2038342/91‑1 (MQ=255)
                                                taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                                            >  2:1432297/1‑91 (MQ=255)
                                                 aCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTAc                                           <  2:3096476/91‑1 (MQ=255)
                                                  cTTTGTTTAACGAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGGTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                          <  2:3623763/91‑1 (MQ=255)
                                                  cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                          >  1:311250/1‑91 (MQ=255)
                                                  cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                          >  1:1453327/1‑91 (MQ=255)
                                                   tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                                         >  1:2481988/1‑91 (MQ=255)
                                                   tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                                         >  1:418465/1‑91 (MQ=255)
                                                   tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                                         <  1:3028013/91‑1 (MQ=255)
                                                    ttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaa                                        >  1:3492174/1‑91 (MQ=255)
                                                       tttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGt                                     >  2:3402720/1‑91 (MQ=255)
                                                        ttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGtt                                    <  2:1625779/91‑1 (MQ=255)
                                                         tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGttt                                   <  1:2594824/91‑1 (MQ=255)
                                                         tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTAGTACAAAGTAATAAAGttt                                   <  2:2763081/91‑1 (MQ=25)
                                                           aaaaGCACCATGCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                                 >  2:795440/3‑91 (MQ=21)
                                                           aCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                                 <  1:1432297/91‑1 (MQ=255)
                                                            cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                                <  1:466757/91‑1 (MQ=255)
                                                            cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                                >  1:3681920/1‑91 (MQ=255)
                                                            cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                                >  1:2577730/1‑91 (MQ=255)
                                                             aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                               <  2:1448611/91‑1 (MQ=255)
                                                              aGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatta                              >  2:2550121/1‑91 (MQ=255)
                                                                cacagtgCTTGTTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat                            >  1:1678178/5‑91 (MQ=255)
                                                                cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat                            <  1:1826906/91‑1 (MQ=255)
                                                                 accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt                           >  1:1943303/5‑91 (MQ=25)
                                                                 aCTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt                           >  1:8474/1‑91 (MQ=255)
                                                                  ccatgCTTATATTAAAATAAGCACAGCGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt                          >  1:1890876/4‑91 (MQ=37)
                                                                  cTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt                          >  1:2076035/1‑91 (MQ=255)
                                                                   tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt                         >  2:243229/1‑91 (MQ=255)
                                                                   tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt                         >  2:1676353/1‑91 (MQ=255)
                                                                   tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt                         >  1:583188/1‑91 (MQ=255)
                                                                    gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa                        >  2:2372175/1‑91 (MQ=255)
                                                                    gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa                        >  1:2693325/1‑91 (MQ=255)
                                                                        ttATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGa                    <  2:2868749/91‑1 (MQ=255)
                                                                         tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg                   >  1:2761704/1‑91 (MQ=255)
                                                                            attAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaa                >  1:1066151/2‑91 (MQ=255)
                                                                             ttAATATAAGCATGGTGCTTTTTGTGTTTTTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaa               >  2:3350444/1‑91 (MQ=255)
                                                                             ttAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaa               <  2:2122823/91‑1 (MQ=255)
                                                                              tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa              <  1:3806007/91‑1 (MQ=255)
                                                                              tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa              >  1:3691537/1‑91 (MQ=255)
                                                                              tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa              >  2:494394/1‑91 (MQ=255)
                                                                              tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa              <  1:2171068/91‑1 (MQ=255)
                                                                                   taAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACAt         <  2:1901153/91‑1 (MQ=255)
                                                                                   taAGCATGGTGCTTGTTGTGTTATTATAGAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAcgt         <  2:1076044/91‑3 (MQ=255)
                                                                                    aaGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACAtt        <  1:3781688/91‑1 (MQ=255)
                                                                                       cATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTaaa     <  2:3504978/91‑1 (MQ=255)
                                                                                        aTGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAAc    >  1:647569/1‑91 (MQ=255)
                                                                                         tGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACt   >  2:3512772/1‑91 (MQ=255)
                                                                                         tGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACt   >  2:1026613/1‑91 (MQ=255)
                                                                                          ggTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACTc  >  1:2334426/1‑91 (MQ=255)
                                                                                          ggTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACTc  >  1:2094574/1‑91 (MQ=255)
                                                                                          |                                                                                          
GAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGACATTAAACTC  >  USA300TCH1516_ALE/536787‑536967

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: