Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 536,863 T→A 7.9% 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 7.9% 439.4 / 21.1 127intergenic (+56/‑255)rplY/pth50S ribosomal protein L25/Peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (63/54);  new base A (8/2);  total (71/56)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTA  >  USA300TCH1516_ALE/536773‑536952
                                                                                          |                                                                                         
aaaGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTAtt                                                                                           <  1:793557/91‑1 (MQ=255)
  aGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTAtttt                                                                                         <  1:299346/91‑1 (MQ=255)
  aGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTAtttt                                                                                         <  1:1577351/91‑1 (MQ=255)
   gCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTa                                                                                        >  1:687255/1‑91 (MQ=255)
   gCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTa                                                                                        >  1:1793500/1‑91 (MQ=255)
     aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAat                                                                                      >  1:3958225/1‑91 (MQ=255)
     aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAat                                                                                      >  1:3958692/1‑91 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTGAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                     >  1:3301614/1‑91 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                     >  1:3408001/1‑91 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                     >  1:854304/1‑91 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                     >  1:1492939/1‑91 (MQ=255)
         aagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATaa                                                                                  >  1:789882/1‑91 (MQ=255)
         aagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGCTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATaa                                                                                  >  1:4140822/1‑91 (MQ=255)
          agaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAg                                                                                 >  1:2131768/1‑91 (MQ=255)
            aCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCa                                                                               <  1:1491328/91‑1 (MQ=255)
            aCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCa                                                                               <  1:1204366/91‑1 (MQ=255)
              gaagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                             <  1:1605561/91‑1 (MQ=255)
              gaagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAACTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                             >  1:3595460/1‑91 (MQ=255)
               aagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATgg                                                                            >  1:1380849/1‑91 (MQ=255)
                agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGt                                                                           <  1:2589894/91‑1 (MQ=255)
                agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCacagt                                                                           <  1:2471950/91‑5 (MQ=25)
                 gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                          <  1:1646866/91‑1 (MQ=255)
                  aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                         <  1:202018/91‑1 (MQ=255)
                     aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                      >  1:3014134/1‑91 (MQ=255)
                     aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                      >  1:3014675/1‑91 (MQ=255)
                     aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                      >  1:1853884/1‑91 (MQ=255)
                     aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                      >  1:269249/1‑91 (MQ=255)
                     aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                      <  1:868746/91‑1 (MQ=255)
                      aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                     <  1:2566856/91‑1 (MQ=255)
                      aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                     <  1:2317564/91‑1 (MQ=255)
                      aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                     <  1:1163210/91‑1 (MQ=255)
                       aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                    >  1:1377322/1‑91 (MQ=255)
                       aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                    >  1:1182640/1‑91 (MQ=255)
                        cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                   <  1:169388/91‑1 (MQ=255)
                        cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                   >  1:4117885/1‑91 (MQ=255)
                          gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                 >  1:2046189/1‑91 (MQ=255)
                           aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                >  1:4112711/1‑91 (MQ=255)
                            agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                               <  1:51634/91‑1 (MQ=255)
                            agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                               <  1:119198/91‑1 (MQ=255)
                             gagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtta                                                              >  1:2018899/1‑91 (MQ=255)
                              agTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttat                                                             <  1:687442/91‑1 (MQ=255)
                                tAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatta                                                           >  1:1947535/1‑91 (MQ=255)
                                tAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatta                                                           >  1:1844355/1‑91 (MQ=255)
                                 aaTTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTtattat                                                          <  1:2429295/91‑1 (MQ=255)
                                   ttttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaa                                                        <  1:2877748/91‑1 (MQ=255)
                                   ttttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaa                                                        <  1:2878797/91‑1 (MQ=255)
                                    tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                       >  1:3260700/1‑91 (MQ=255)
                                    tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                       >  1:374619/1‑91 (MQ=255)
                                       aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                    <  1:4001719/91‑1 (MQ=255)
                                        aTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtt                                                   <  1:2527354/91‑1 (MQ=255)
                                        aTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtt                                                   <  1:367487/91‑1 (MQ=255)
                                         tCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtta                                                  >  1:1615835/1‑91 (MQ=255)
                                         tCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtta                                                  >  1:500165/1‑91 (MQ=255)
                                           tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                <  1:1877373/91‑1 (MQ=255)
                                             ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                              <  1:3746850/91‑1 (MQ=255)
                                                   tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                        <  1:115980/91‑1 (MQ=255)
                                                   tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                        <  1:1708203/91‑1 (MQ=255)
                                                    aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                       <  1:1028991/91‑1 (MQ=255)
                                                    aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                       <  1:2736064/91‑1 (MQ=255)
                                                    aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                       <  1:4008572/91‑1 (MQ=255)
                                                    aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                       >  1:3647925/1‑91 (MQ=255)
                                                    aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                       >  1:1739314/1‑91 (MQ=255)
                                                        tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTAc                                   >  1:3804462/1‑91 (MQ=255)
                                                         ttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACt                                  >  1:770689/1‑91 (MQ=255)
                                                          tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                 <  1:3481445/91‑1 (MQ=255)
                                                             ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                              <  1:972601/91‑1 (MQ=255)
                                                             ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                              >  1:3169269/1‑91 (MQ=255)
                                                              taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                             >  1:2783409/1‑91 (MQ=255)
                                                                cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                           >  1:1344801/1‑91 (MQ=255)
                                                                cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                           >  1:2346977/1‑91 (MQ=255)
                                                                 tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                          >  1:353841/1‑91 (MQ=255)
                                                                   tGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaaa                        >  1:2342343/1‑91 (MQ=255)
                                                                     tttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGt                      <  1:2187784/91‑1 (MQ=255)
                                                                      ttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGtt                     >  1:1944938/1‑91 (MQ=255)
                                                                       tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGttt                    >  1:744307/1‑91 (MQ=255)
                                                                       tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGttt                    <  1:2792901/91‑1 (MQ=255)
                                                                        aaCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTAGTACAAAGTAATAAAGTtta                   <  1:4090881/91‑1 (MQ=25)
                                                                         aCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                  <  1:3425541/91‑1 (MQ=255)
                                                                          cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                 <  1:3791256/91‑1 (MQ=255)
                                                                           aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                >  1:16466/1‑91 (MQ=255)
                                                                           aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                >  1:17339/1‑91 (MQ=255)
                                                                             gCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaattaa              <  1:406619/91‑1 (MQ=255)
                                                                             gCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaattaa              >  1:1559965/1‑91 (MQ=255)
                                                                              cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat             >  1:164605/1‑91 (MQ=255)
                                                                               accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt            >  1:2585883/5‑91 (MQ=25)
                                                                               accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt            >  1:209919/5‑91 (MQ=25)
                                                                               accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt            >  1:1378983/5‑91 (MQ=25)
                                                                                ccatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt           >  1:1483079/4‑91 (MQ=37)
                                                                                ccatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt           >  1:541923/4‑91 (MQ=37)
                                                                                ccatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt           >  1:887862/4‑91 (MQ=37)
                                                                                 tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt          >  1:2983988/1‑91 (MQ=255)
                                                                                   tgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg        <  1:3671996/91‑1 (MQ=255)
                                                                                   tgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg        >  1:3915028/1‑91 (MQ=255)
                                                                                    gCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGt       >  1:2822263/1‑91 (MQ=255)
                                                                                    gCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGt       >  1:969383/1‑91 (MQ=255)
                                                                                     cTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTg      <  1:1724007/91‑1 (MQ=255)
                                                                                     cTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTg      <  1:990802/91‑1 (MQ=255)
                                                                                       tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg    <  1:1855566/91‑1 (MQ=255)
                                                                                       tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg    <  1:1855685/91‑1 (MQ=255)
                                                                                         ttttAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTa  <  1:3131642/91‑1 (MQ=255)
                                                                                          |                                                                                         
AAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTA  >  USA300TCH1516_ALE/536773‑536952

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