Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 419,965 A→T 7.0% intergenic (+183/+83) rpsR → / ← USA300HOU_RS01950 30S ribosomal protein S18/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE419,9650AT 7.0% 374.7 / 13.1 112intergenic (+183/+83)rpsR/USA300HOU_RS0195030S ribosomal protein S18/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (47/55);  new base T (4/4);  total (51/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

TGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACT  >  USA300TCH1516_ALE/419875‑420054
                                                                                          |                                                                                         
tGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaa                                                                                           <  1:2088828/91‑1 (MQ=255)
 gTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaat                                                                                          <  1:3792561/91‑1 (MQ=255)
  ttttGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaata                                                                                         >  1:3599855/1‑91 (MQ=255)
   tttGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaataa                                                                                        <  1:2435117/91‑1 (MQ=255)
   tttGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaataa                                                                                        <  1:3165047/91‑1 (MQ=255)
   tttGTGCTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaataa                                                                                        >  1:3223294/1‑91 (MQ=255)
    ttGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaataaa                                                                                       >  1:3069971/1‑91 (MQ=255)
     tgtgTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATaaaa                                                                                      >  1:3241515/1‑91 (MQ=255)
      gtgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAAt                                                                                     >  1:219043/1‑91 (MQ=255)
       tgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATg                                                                                    >  1:3998236/1‑91 (MQ=255)
       tgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATg                                                                                    >  1:2140459/1‑91 (MQ=255)
       tgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATg                                                                                    >  1:313979/1‑91 (MQ=255)
        gtTGGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTATTTAAAATGc                                                                                   <  1:3739479/91‑1 (MQ=255)
        gtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGc                                                                                   >  1:2648868/1‑91 (MQ=255)
           cGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAGTAAAATGCATa                                                                                <  1:3037062/91‑1 (MQ=255)
            ggAATTTGCACAAAGATATGTTTATATTGCAAGAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAg                                                                               <  1:3532853/91‑1 (MQ=255)
            ggAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAg                                                                               <  1:3436322/91‑1 (MQ=255)
            ggAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAg                                                                               >  1:1794684/1‑91 (MQ=255)
            ggAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAg                                                                               <  1:3041199/91‑1 (MQ=255)
             gAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGc                                                                              >  1:2572035/1‑91 (MQ=255)
             gAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGc                                                                              <  1:3972105/91‑1 (MQ=255)
              aaTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTATGCATTTTATTTAAAATGCGTAgtt                                                                             <  1:2051329/91‑3 (MQ=37)
               aTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTATTTAAAATGCTTAgttc                                                                            <  1:3222864/91‑4 (MQ=38)
               aTTTGCACAAAGATATGCTTATATTGCAAAAAAAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGctc                                                                            <  1:3554034/91‑1 (MQ=255)
                tttGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGctct                                                                           <  1:3051083/91‑1 (MQ=255)
                  tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttc                                                                         >  1:972783/1‑91 (MQ=255)
                  tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttc                                                                         >  1:3774747/1‑91 (MQ=255)
                  tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttc                                                                         >  1:534429/1‑91 (MQ=255)
                  tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttc                                                                         >  1:938021/1‑91 (MQ=255)
                   gCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttct                                                                        <  1:1225027/91‑1 (MQ=255)
                      caAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCtttt                                                                     <  1:726170/91‑1 (MQ=255)
                      caAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCtttt                                                                     <  1:2552511/91‑1 (MQ=255)
                      caAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCtttt                                                                     <  1:2552781/91‑1 (MQ=255)
                      caAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCtttt                                                                     <  1:2552823/91‑1 (MQ=255)
                      caAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCtttt                                                                     <  1:2552939/91‑1 (MQ=255)
                       aaaGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCttttt                                                                    <  1:2907633/91‑1 (MQ=255)
                       aaaGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCttttt                                                                    <  1:2907519/91‑1 (MQ=255)
                          gATATGTTTATATTGCAAAAATAATATGAATTTGGATGCATAAAAGAAGAGCTATGCATTTTATTTAAAATGCGTAGTTCTttttttatgc                                                                 <  1:633380/91‑5 (MQ=2)
                            tatGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGc                                                               >  1:2214409/1‑91 (MQ=255)
                             atGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCa                                                              <  1:629041/91‑1 (MQ=255)
                             atGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCa                                                              >  1:3884480/1‑91 (MQ=255)
                              tGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCAt                                                             <  1:1218000/91‑1 (MQ=255)
                              tGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCAt                                                             <  1:1661173/91‑1 (MQ=255)
                                 ttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACg                                                          >  1:106523/1‑91 (MQ=255)
                                   atatTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGaa                                                        >  1:945944/1‑91 (MQ=255)
                                   atatTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGaa                                                        >  1:3355209/1‑91 (MQ=255)
                                    tatTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAAt                                                       <  1:3667112/91‑1 (MQ=255)
                                    tatTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAAt                                                       >  1:3655686/1‑91 (MQ=255)
                                        gCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaa                                                   >  1:2993469/1‑91 (MQ=255)
                                        gCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaa                                                   >  1:1513920/1‑91 (MQ=255)
                                         cAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaaa                                                  <  1:360785/91‑1 (MQ=255)
                                          aaaaaTAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAAt                                                 <  1:596555/91‑1 (MQ=255)
                                          aaaaaTAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAAt                                                 <  1:2209757/91‑1 (MQ=255)
                                             aataatATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAAc                                              <  1:342489/91‑1 (MQ=255)
                                                aataTGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcg                                           >  1:2634/1‑91 (MQ=255)
                                                 atatGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgc                                          <  1:1913883/91‑1 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        >  1:1197664/1‑91 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        <  1:841640/91‑1 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        <  1:3416693/91‑1 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        <  1:3604838/91‑1 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        >  1:1510709/1‑91 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        >  1:1510355/1‑91 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        >  1:1510249/1‑91 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                        <  1:2886128/91‑1 (MQ=255)
                                                   atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTAGCATACGAATTAAAATAACTCGCga                                        >  1:3857876/1‑91 (MQ=255)
                                                     gAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                                      >  1:3269210/1‑91 (MQ=255)
                                                     gAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                                      >  1:3269250/1‑91 (MQ=255)
                                                     gAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                                      <  1:708228/91‑1 (MQ=255)
                                                             aTGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCCCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAg                              >  1:1291704/1‑91 (MQ=255)
                                                              tGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGt                             >  1:2216199/1‑91 (MQ=255)
                                                               gCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtc                            <  1:4084906/91‑1 (MQ=255)
                                                                cATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtct                           <  1:2956669/91‑1 (MQ=255)
                                                                 aTAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtctc                          <  1:207002/91‑1 (MQ=255)
                                                                  tAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTctct                         <  1:1268201/91‑1 (MQ=255)
                                                                  tAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTctct                         <  1:1761253/91‑1 (MQ=255)
                                                                  tAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTctct                         <  1:1761594/91‑1 (MQ=255)
                                                                    aaaaaaGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCttt                       >  1:4050294/1‑91 (MQ=255)
                                                                      agaaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTcc                     >  1:3952187/3‑91 (MQ=18)
                                                                      aaaaGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTcc                     >  1:3644181/1‑91 (MQ=255)
                                                                      aaaaGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTcc                     >  1:485074/1‑91 (MQ=255)
                                                                       gaaGAGCTATGCATTTTATTTAAAATGCGTAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCt                    >  1:1993006/2‑91 (MQ=25)
                                                                         aGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCa                  >  1:1783716/1‑91 (MQ=255)
                                                                         aGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCa                  <  1:3880582/91‑1 (MQ=255)
                                                                          gAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCAc                 >  1:2444397/1‑91 (MQ=255)
                                                                          gAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCAc                 <  1:429205/91‑1 (MQ=255)
                                                                           aaCTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACt                >  1:1317410/1‑91 (MQ=255)
                                                                           aaCTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACt                <  1:3417028/91‑1 (MQ=255)
                                                                            aCTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTa               >  1:1903865/1‑91 (MQ=255)
                                                                            aCTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTa               >  1:2464254/1‑91 (MQ=255)
                                                                              tACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGa             <  1:1210673/91‑1 (MQ=255)
                                                                              tACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTATCTTTCCTCACTAGa             <  1:3540265/91‑1 (MQ=255)
                                                                               atgCATTTTATTTAAAATGCGTAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGAt            >  1:2003271/3‑91 (MQ=38)
                                                                                   aTTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGtt        <  1:2871079/91‑1 (MQ=255)
                                                                                   aTTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGtt        <  1:3607079/91‑1 (MQ=255)
                                                                                    ttttAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGttt       <  1:2336697/91‑1 (MQ=255)
                                                                                    ttttAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGttt       <  1:2907618/91‑1 (MQ=255)
                                                                                      ttAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTtat     <  1:1363862/91‑1 (MQ=255)
                                                                                         aataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTGACTAGATAGTTTATACt  >  1:3554101/1‑91 (MQ=255)
                                                                                         aataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACt  <  1:1951541/91‑1 (MQ=255)
                                                                                         aataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCAATAAGTCTCTTTCCTCACTAGATAGTTTATACt  >  1:3399970/1‑91 (MQ=255)
                                                                                          |                                                                                         
TGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACT  >  USA300TCH1516_ALE/419875‑420054

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