Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 536,863 T→A 9.0% intergenic (+56/‑255) 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,8630TA 9.0% 753.6 / 47.2 233intergenic (+56/‑255)rplY/pth50S ribosomal protein L25/Peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (110/101);  new base A (9/12);  total (119/114)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 238 total aligned reads displayed.

AAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTA  >  USA300TCH1516_ALE/536773‑536952
                                                                                          |                                                                                         
aaaGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTAtt                                                                                           <  1:844082/91‑1 (MQ=255)
aaaGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTAtt                                                                                           <  2:2823311/91‑1 (MQ=255)
     aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAat                                                                                      <  1:1614218/91‑1 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                     <  1:2552655/91‑1 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                     >  2:2171214/1‑91 (MQ=255)
       agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                    >  1:1210767/1‑91 (MQ=255)
        gaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAAtata                                                                                   >  2:2119120/1‑91 (MQ=255)
         aagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATaa                                                                                  <  2:2359605/91‑1 (MQ=38)
          agaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAg                                                                                 <  2:2508377/91‑1 (MQ=255)
          agaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAg                                                                                 <  2:2915724/91‑1 (MQ=38)
           gaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGc                                                                                <  1:42956/91‑1 (MQ=255)
             cGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCAt                                                                              >  1:242212/1‑91 (MQ=255)
             cGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCac                                                                              <  1:775094/91‑2 (MQ=37)
             cGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCac                                                                              <  2:591991/91‑2 (MQ=37)
              gaagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCaca                                                                             <  2:3821165/91‑3 (MQ=37)
                agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAAtata                                                                                   >  2:3633014/1‑83 (MQ=255)
                agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCacagt                                                                           <  2:731996/91‑5 (MQ=25)
                 gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                          <  1:1490433/91‑1 (MQ=255)
                 gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                          <  1:1623441/91‑1 (MQ=255)
                 gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                          <  1:471207/91‑1 (MQ=255)
                  aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                         >  2:955786/1‑91 (MQ=255)
                   aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                        <  1:1892314/91‑1 (MQ=255)
                   aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                        <  1:3513473/91‑1 (MQ=255)
                    aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                       >  2:2035301/1‑91 (MQ=255)
                    aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                       >  2:2065200/1‑91 (MQ=255)
                     aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                      <  1:2948012/91‑1 (MQ=255)
                     aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTGAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                      <  2:3186025/91‑1 (MQ=255)
                      aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                     >  1:762854/1‑91 (MQ=255)
                       aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                    >  2:1649562/1‑91 (MQ=255)
                       aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                    <  2:1550054/91‑1 (MQ=255)
                       aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                    <  2:1032139/91‑1 (MQ=255)
                        cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                   <  1:2684606/91‑1 (MQ=255)
                         tGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgt                                                                  >  2:2018183/1‑91 (MQ=255)
                          gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                 >  2:3772481/1‑91 (MQ=255)
                            agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                               >  2:707323/1‑91 (MQ=255)
                              agTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttat                                                             <  2:1225903/91‑1 (MQ=255)
                               gTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAACCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                            >  1:3078602/1‑91 (MQ=255)
                                  aTTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATtata                                                         >  1:3749180/1‑91 (MQ=255)
                                    tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                       >  1:2868942/1‑91 (MQ=255)
                                      tAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                     <  1:1551907/91‑1 (MQ=255)
                                       aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                    >  2:1435378/1‑91 (MQ=255)
                                       aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                    >  2:415334/1‑91 (MQ=255)
                                       aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                    <  1:2850635/91‑1 (MQ=255)
                                       aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                    <  2:3538774/91‑1 (MQ=255)
                                         tCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtta                                                  <  1:2612839/91‑1 (MQ=255)
                                           tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                >  1:2641524/1‑91 (MQ=255)
                                             ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                              <  1:1649562/91‑1 (MQ=255)
                                             ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                              >  1:2515773/1‑91 (MQ=255)
                                              tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                             >  2:649090/1‑91 (MQ=255)
                                                cATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAAc                                           >  2:1466626/1‑91 (MQ=255)
                                                   tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                        >  1:3141576/1‑91 (MQ=255)
                                                   tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTCAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                        >  1:3360596/1‑91 (MQ=255)
                                                     aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttag                                      <  1:953403/91‑2 (MQ=255)
                                                       gTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatta                                    <  1:1050197/91‑1 (MQ=255)
                                                        tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                             <  2:1195792/49‑1 (MQ=255)
                                                        tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTAc                                   >  1:1914217/1‑91 (MQ=255)
                                                          tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                 >  2:2863887/1‑91 (MQ=255)
                                                          tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                 <  2:3209058/91‑1 (MQ=255)
                                                           ttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACttt                                <  2:3714216/91‑1 (MQ=255)
                                                            tataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTg                               <  1:202338/91‑1 (MQ=255)
                                                            tataCTTTGTTTAACAAGCACTATGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTg                               <  2:3121959/91‑1 (MQ=255)
                                                             ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAGAGCTTAATTAAACTTTATTACTTTGt                              <  2:2218806/91‑1 (MQ=255)
                                                             ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                              <  2:1731610/91‑1 (MQ=255)
                                                              taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                             <  1:3681623/91‑1 (MQ=255)
                                                                cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                           <  2:3623963/91‑1 (MQ=255)
                                                                 tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                          <  1:1097773/91‑1 (MQ=255)
                                                                 tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                          >  1:3802781/1‑91 (MQ=255)
                                                                 tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                          >  1:2481988/1‑91 (MQ=255)
                                                                   tGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaaa                        >  1:3250442/1‑91 (MQ=255)
                                                                     tttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGt                      >  2:3402720/1‑91 (MQ=255)
                                                                      ttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGtt                     <  2:2646996/91‑1 (MQ=255)
                                                                      ttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGtt                     <  2:1625779/91‑1 (MQ=255)
                                                                       tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTAGTACAAAGTAATAAAGttt                    <  2:2763081/91‑1 (MQ=25)
                                                                         aaaaGCACCATGCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                  >  2:795440/3‑91 (MQ=21)
                                                                         aCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                  >  1:2170579/1‑91 (MQ=255)
                                                                          cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGttatac                 >  1:1434514/1‑87 (MQ=255)
                                                                          cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                 <  2:747407/91‑1 (MQ=255)
                                                                          cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                 <  1:660099/91‑1 (MQ=255)
                                                                          cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                 >  1:3681920/1‑91 (MQ=255)
                                                                           aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                <  2:1390381/91‑1 (MQ=255)
                                                                              cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTACAGTTtaattaat             <  1:2986159/91‑1 (MQ=255)
                                                                              cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat             >  2:2553406/1‑91 (MQ=255)
                                                                               accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt            >  1:1943303/5‑91 (MQ=25)
                                                                               accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt            >  2:2284791/5‑91 (MQ=25)
                                                                               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)
                                                                                 tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt          >  1:583188/1‑91 (MQ=255)
                                                                                  gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa         >  1:1752011/1‑91 (MQ=255)
                                                                                  gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa         >  1:2693325/1‑91 (MQ=255)
                                                                                   tgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg        >  1:901108/1‑91 (MQ=255)
                                                                                   tgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg        >  1:3545520/1‑91 (MQ=255)
                                                                                   tgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg        >  2:2592269/1‑91 (MQ=38)
                                                                                     cTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTg      >  2:1973834/1‑91 (MQ=255)
                                                                                      ttATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGa     <  2:2868368/91‑1 (MQ=255)
                                                                                       tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg    >  1:2761704/1‑91 (MQ=255)
                                                                                       tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg    >  1:2583562/1‑91 (MQ=255)
                                                                                         ttttAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTa  >  2:2604225/1‑91 (MQ=255)
                                                                                          |                                                                                         
AAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTA  >  USA300TCH1516_ALE/536773‑536952

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