Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 584,121 G→T intergenic (+205/‑67) rplA → / → rplJ 50S ribosomal protein L1/50S ribosomal protein L10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE584,1210GT98.6% 499.8 / ‑1.7 144intergenic (+205/‑67)rplA/rplJ50S ribosomal protein L1/50S ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (75/67);  total (76/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 144 total aligned reads displayed.

AGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAG  >  USA300TCH1516_ALE/584033‑584206
                                                                                        |                                                                                      
aGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTg                                                                                      >  2:623264/1‑91 (MQ=255)
 ggAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTgg                                                                                     >  2:769799/1‑91 (MQ=255)
 ggAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTAGCCTTgg                                                                                     >  1:1545795/1‑91 (MQ=255)
  gAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTggg                                                                                    >  2:2013541/1‑91 (MQ=255)
   aGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTggga                                                                                   >  2:890623/1‑90 (MQ=255)
   aGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGt                                                                                   >  2:890554/1‑91 (MQ=255)
    gTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTa                                                                                  >  1:1453514/1‑91 (MQ=255)
        tttatAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGaaa                                                                              >  2:2792156/1‑91 (MQ=255)
         ttatAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAg                                                                             <  2:510718/91‑1 (MQ=255)
         ttatAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAg                                                                             <  1:283657/91‑1 (MQ=255)
          tataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGt                                                                            >  1:408075/1‑91 (MQ=255)
          tataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGt                                                                            >  1:2713645/1‑91 (MQ=255)
           ataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                           <  2:2459177/91‑1 (MQ=255)
           ataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                           >  1:1106388/1‑91 (MQ=255)
           ataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                           >  2:1073299/1‑91 (MQ=255)
           ataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                           <  1:933958/91‑1 (MQ=255)
           ataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTg                                                                           >  1:418552/1‑91 (MQ=255)
              aCTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCtt                                                                        >  1:1089825/1‑91 (MQ=255)
                ttAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCtttt                                                                      >  1:949329/1‑91 (MQ=255)
                ttAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCtttt                                                                      <  1:2080260/91‑1 (MQ=255)
                   aaaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCttttttt                                                                   >  2:807568/1‑91 (MQ=255)
                    aaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTtta                                                                  >  1:2412307/1‑91 (MQ=255)
                     aTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTttat                                                                 <  1:1717297/91‑1 (MQ=255)
                     aTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTttat                                                                 <  1:1716701/91‑1 (MQ=255)
                      tttATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTttatt                                                                >  1:2413318/1‑91 (MQ=255)
                      tttATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTttatt                                                                >  2:357852/1‑91 (MQ=255)
                        tATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTaa                                                              >  2:1918273/1‑91 (MQ=255)
                        tATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTaa                                                              >  2:2013110/1‑91 (MQ=255)
                         aTCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAt                                                             >  2:813535/1‑91 (MQ=255)
                         aTCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAt                                                             <  1:483807/91‑1 (MQ=255)
                         aTCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAt                                                             >  2:24122/1‑91 (MQ=255)
                          tCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtt                                                            >  1:3103174/1‑91 (MQ=255)
                          tCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtt                                                            >  2:1856496/1‑91 (MQ=255)
                          tCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtt                                                            >  2:2626085/1‑91 (MQ=255)
                          tCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtt                                                            >  2:2626086/1‑91 (MQ=255)
                          tCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtt                                                            >  2:2626232/1‑91 (MQ=255)
                          tCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtt                                                            <  2:2870568/91‑1 (MQ=255)
                           ccTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAttt                                                          <  2:2902693/91‑1 (MQ=255)
                           ccTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAttt                                                          <  2:3095513/91‑1 (MQ=255)
                           ccTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAttt                                                          <  2:1468638/91‑1 (MQ=255)
                           ccTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAttt                                                          >  1:759713/1‑91 (MQ=255)
                            cTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtttt                                                         >  1:675232/1‑91 (MQ=255)
                                  aGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTaaaaaa                                                   >  2:1699958/1‑91 (MQ=255)
                                  aGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTaaaaaa                                                   <  1:2770810/91‑1 (MQ=255)
                                   ggCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTaaaaaaa                                                  <  1:2459236/91‑1 (MQ=255)
                                   ggCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTGATTAATTTTaaaaaaa                                                  >  2:2381443/1‑91 (MQ=255)
                                      tAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCa                                               <  2:93657/91‑1 (MQ=255)
                                      tAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCa                                               <  2:1493196/91‑1 (MQ=255)
                                      tAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCa                                               <  1:2399423/91‑1 (MQ=255)
                                      tAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCa                                               >  2:263168/1‑91 (MQ=255)
                                         aaTTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCa                                            <  1:1369746/91‑1 (MQ=255)
                                           ttGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAGGTGCTTTTTTTAAAATTAATAAAAAAAGCActttc                                          <  1:428201/91‑5 (MQ=18)
                                             gACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAACAAAGCACCaaaaa                                        >  2:1128646/1‑91 (MQ=255)
                                              aCTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAt                                       <  2:538079/91‑1 (MQ=255)
                                               cTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAATCACCAAAAAtt                                      >  1:2809689/1‑91 (MQ=255)
                                               cTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAtt                                      >  1:982565/1‑91 (MQ=255)
                                               cTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAtt                                      <  2:2136296/91‑1 (MQ=255)
                                               cTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAtt                                      <  1:666156/91‑1 (MQ=255)
                                               cTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAtt                                      <  2:2134032/91‑1 (MQ=255)
                                                   aaCGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTaaa                                  <  2:152573/91‑1 (MQ=255)
                                                   aaCGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTaaa                                  <  2:428269/91‑1 (MQ=255)
                                                    aCGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAAt                                 >  2:2803195/1‑91 (MQ=255)
                                                    aCGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAAt                                 <  1:3014550/91‑1 (MQ=255)
                                                    aCGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAAt                                 <  2:945551/91‑1 (MQ=255)
                                                       tgatgaTCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGa                              >  2:263424/1‑91 (MQ=255)
                                                       tgatgaTCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGa                              >  1:2840722/1‑91 (MQ=255)
                                                        gatgatCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAg                             <  1:401120/91‑1 (MQ=255)
                                                        gatgatCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTGTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAg                             <  1:90160/91‑1 (MQ=255)
                                                         atgatctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAgg                            <  2:2503101/91‑1 (MQ=255)
                                                           gatctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGgtg                          >  1:1869286/1‑91 (MQ=255)
                                                           gatctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGgtg                          >  2:1157956/1‑91 (MQ=255)
                                                           gatctatgatctTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGgtg                          <  2:1274655/91‑1 (MQ=255)
                                                              ctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCt                       <  1:1810381/91‑1 (MQ=255)
                                                              ctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCt                       <  1:420991/91‑1 (MQ=255)
                                                               tatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTg                      <  1:495392/91‑1 (MQ=255)
                                                                atgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGCGTCTGa                     >  2:243242/1‑91 (MQ=255)
                                                                  gatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAAt                   >  1:1172818/1‑91 (MQ=255)
                                                                    tctTTCAAGCACTTTTTTCCTTGGGTATAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGt                 >  2:2212565/1‑91 (MQ=255)
                                                                    tctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGt                 >  2:2212478/1‑91 (MQ=255)
                                                                     ctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTc                >  1:634278/1‑91 (MQ=255)
                                                                     ctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTc                <  1:270221/91‑1 (MQ=255)
                                                                     ctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTc                >  1:637340/1‑91 (MQ=255)
                                                                      tttCAAGCGCTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCTAAAATTTAAATGGAGGTGTCTGAATGTCt               <  2:2697353/91‑1 (MQ=255)
                                                                       ttCAAGCACTTTTTGCCTTGGGGAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTg              <  2:2496903/91‑1 (MQ=255)
                                                                        tCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGc             <  1:678970/91‑1 (MQ=255)
                                                                        tCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGc             >  1:3053924/1‑91 (MQ=255)
                                                                        tCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGc             <  1:1580441/91‑1 (MQ=255)
                                                                         cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAATATTTAAATGGAGGTGTCTGAATGTCTGCt            >  1:2488674/1‑91 (MQ=255)
                                                                         cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCt            >  2:1612391/1‑91 (MQ=255)
                                                                         cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCt            >  1:2521067/1‑91 (MQ=255)
                                                                         cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCt            >  1:981436/1‑91 (MQ=255)
                                                                          aaGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTa           <  1:1473459/91‑1 (MQ=255)
                                                                           aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTTATTTTAAAAAAAGAACCAAAAATTTAAATCGTGGTGTCTGAATGTCTGCTAt          >  2:269728/1‑91 (MQ=255)
                                                                           aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTAt          <  2:1886397/91‑1 (MQ=255)
                                                                           aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTAt          <  1:2035139/91‑1 (MQ=255)
                                                                           aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTAt          >  2:2252802/1‑91 (MQ=255)
                                                                            gCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATc         >  1:1050597/1‑91 (MQ=255)
                                                                            gCACTGTTTGCCTTGGGTAGAAAGTTCTTTTTTTATTAAGTTTAAAAAAAGCAGCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATc         <  1:2886296/91‑1 (MQ=255)
                                                                                 ttttGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGa    >  1:157440/1‑91 (MQ=255)
                                                                                   ttGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAg  >  2:722239/1‑91 (MQ=255)
                                                                                        |                                                                                      
AGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAG  >  USA300TCH1516_ALE/584033‑584206

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