Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 16.5% intergenic (+711/‑6052) recR → / → speA_2 Recombination protein RecR/Arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE512,6430AG16.5% 582.9 / 100.3 243intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (120/83);  new base G (22/18);  total (142/101)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 265 total aligned reads displayed.

ATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAC  >  USA300TCH1516_ALE/512553‑512733
                                                                                           |                                                                                          
  ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                          <  1:558001/91‑1 (MQ=255)
    gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                        <  1:3311644/91‑1 (MQ=255)
    gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                        >  1:422510/1‑91 (MQ=255)
        atgCAAACAAAACAACTCAACTTAAGGAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATg                                                                                    <  2:1974594/89‑1 (MQ=255)
        aGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATg                                                                                    >  1:1366181/1‑91 (MQ=255)
        aGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATg                                                                                    <  1:803995/91‑1 (MQ=255)
         ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTGAATGTTAATAAAATGt                                                                                   >  2:3781977/1‑91 (MQ=255)
         ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGt                                                                                   >  2:3841761/1‑91 (MQ=255)
           caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat                                                                                 >  1:3434016/1‑91 (MQ=255)
             aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                               >  2:2555034/1‑91 (MQ=255)
             aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                               >  2:3497907/1‑91 (MQ=255)
             aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                               <  2:2105076/91‑1 (MQ=255)
              acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTCTTGACTTAAATGTTAATAAAATGTAtaat                                                                              >  1:3177621/1‑91 (MQ=255)
              acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaat                                                                              >  1:1855987/1‑91 (MQ=255)
              acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaat                                                                              >  1:625352/1‑91 (MQ=255)
               caaaacaaCTCAACTTAAGAAATATCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                             <  2:3685702/91‑1 (MQ=255)
               caaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                             <  1:823054/91‑1 (MQ=255)
                 aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                           >  2:3857216/1‑91 (MQ=255)
                 aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                           >  1:2255880/1‑91 (MQ=255)
                 aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                           >  2:2815581/1‑91 (MQ=255)
                  aacaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaat                                                                          <  1:2748944/91‑1 (MQ=255)
                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                         >  2:3799882/1‑91 (MQ=255)
                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                         <  2:1072221/91‑1 (MQ=255)
                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAACTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                         >  1:1194460/1‑91 (MQ=255)
                     aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCt                                                                       >  2:2591070/1‑91 (MQ=255)
                      actcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTCAATAAAATGTACAATTAATTCtt                                                                      >  2:3016171/1‑91 (MQ=255)
                       ctcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg                                                                     >  2:676841/1‑91 (MQ=255)
                        tcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGt                                                                    <  2:1399980/91‑1 (MQ=255)
                          aactTAAGAAATGACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCg                                                                  <  2:3288016/91‑1 (MQ=255)
                            ctTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                <  2:2318866/91‑1 (MQ=255)
                            ctTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                >  1:1021469/1‑91 (MQ=255)
                             ttAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa                                                               <  1:1797096/91‑1 (MQ=255)
                               aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                             <  2:2148972/91‑1 (MQ=255)
                               aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                             >  2:126945/1‑91 (MQ=255)
                                aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                            >  1:455017/1‑91 (MQ=255)
                                aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                            <  2:1526517/91‑1 (MQ=255)
                                 gAAATAACTTGAATTACTAACGAAAATTAATTATAAAAAAGTTATTGACTTAAATGTTAATAAACTGTATAATTAATTCTTGTCGGTAAGa                                                            <  2:2241869/91‑1 (MQ=255)
                                 gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTACGaa                                                           >  1:1913721/1‑91 (MQ=255)
                                 gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                           >  1:2944331/1‑91 (MQ=255)
                                 gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                           >  1:195708/1‑91 (MQ=255)
                                      aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg                                                      <  2:858582/91‑1 (MQ=255)
                                       aCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGa                                                     <  2:2690946/91‑1 (MQ=255)
                                       aCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGa                                                     >  2:1709999/1‑91 (MQ=255)
                                          tGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                  <  1:1578052/91‑1 (MQ=255)
                                          tGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                  <  2:812328/91‑1 (MQ=255)
                                          tGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                  <  2:1469287/91‑1 (MQ=255)
                                              ttACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGa                                              >  2:2610625/1‑91 (MQ=255)
                                               tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaa                                             >  1:1143107/1‑91 (MQ=255)
                                                aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaa                                            >  2:3154578/1‑91 (MQ=255)
                                                 cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                           >  1:2468511/1‑91 (MQ=255)
                                                 cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                           <  1:2629547/91‑1 (MQ=255)
                                                  tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                          >  1:3046018/1‑91 (MQ=255)
                                                   aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                         >  1:1645983/1‑91 (MQ=255)
                                                    aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                        <  2:3941381/91‑1 (MQ=255)
                                                    aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                        <  1:112592/91‑1 (MQ=255)
                                                       aaaaTTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt                                     >  1:3769064/1‑91 (MQ=255)
                                                        aaaTTTATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAGCTGAATg                                    <  2:575161/91‑1 (MQ=255)
                                                        aaaTTAATTTTAAAAAGTTATTGACTTAAATGTTAATACAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                    >  2:1640940/1‑91 (MQ=255)
                                                          attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAc                                  >  1:1286974/1‑91 (MQ=255)
                                                             aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAat                               <  2:1691048/91‑1 (MQ=255)
                                                             aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAat                               >  1:1447863/1‑91 (MQ=255)
                                                              attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                              <  2:937042/91‑1 (MQ=255)
                                                              attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                              >  2:3459298/1‑91 (MQ=255)
                                                               ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                             >  2:532638/1‑91 (MQ=255)
                                                               ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                             >  1:1507060/1‑91 (MQ=255)
                                                                ataaaaaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                             <  1:1579901/87‑1 (MQ=255)
                                                                   tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAGATGAACATTGAAAACTGAATGACACTATgcg                          >  2:50191/1‑89 (MQ=255)
                                                                     aaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                        >  2:2453917/1‑91 (MQ=255)
                                                                      aaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                       >  2:186052/1‑91 (MQ=255)
                                                                        aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                     >  2:3953402/1‑91 (MQ=255)
                                                                        aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                     >  1:2689497/1‑91 (MQ=255)
                                                                        aGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                     >  1:2391710/1‑91 (MQ=255)
                                                                         gatATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    >  2:2047473/3‑91 (MQ=255)
                                                                         gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    >  2:1878533/1‑91 (MQ=255)
                                                                         gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    >  2:22864/1‑91 (MQ=255)
                                                                         gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    >  1:694810/1‑91 (MQ=255)
                                                                         gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    <  1:2501147/91‑1 (MQ=255)
                                                                         gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    >  2:1352848/1‑91 (MQ=255)
                                                                           tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                  >  2:3602802/1‑91 (MQ=255)
                                                                           tATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                  >  2:3028784/1‑91 (MQ=255)
                                                                           tATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                  <  1:2850231/91‑1 (MQ=255)
                                                                            aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 <  1:77255/91‑1 (MQ=255)
                                                                            aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 >  1:2083518/1‑91 (MQ=255)
                                                                             ttGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                >  1:57113/1‑91 (MQ=255)
                                                                              tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc               >  2:3079826/1‑91 (MQ=255)
                                                                               gACTTAAATCTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc              <  1:3780703/91‑1 (MQ=255)
                                                                                 catAAATGTTGATAAAATGTATAATTAATTCTTGACGATAAAAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACaa            <  2:3355841/89‑1 (MQ=21)
                                                                                 cTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGACAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa            <  1:3953925/91‑1 (MQ=255)
                                                                                 cTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa            >  1:2333702/1‑91 (MQ=255)
                                                                                  ttAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa           <  2:2598411/91‑1 (MQ=255)
                                                                                  ttAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa           <  2:2918451/91‑1 (MQ=255)
                                                                                   tAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa          >  2:1802332/1‑91 (MQ=255)
                                                                                     aaTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa        >  2:1213168/1‑91 (MQ=255)
                                                                                     aaTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa        >  1:3096652/1‑91 (MQ=255)
                                                                                     aaTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa        >  1:3713063/1‑91 (MQ=255)
                                                                                          taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa   <  2:1535555/91‑1 (MQ=255)
                                                                                           aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc  >  1:1652676/1‑91 (MQ=255)
                                                                                           aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc  <  1:3583179/91‑1 (MQ=255)
                                                                                           aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc  <  1:3932348/91‑1 (MQ=255)
                                                                                           aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc  >  1:88613/1‑91 (MQ=255)
                                                                                           |                                                                                          
ATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAC  >  USA300TCH1516_ALE/512553‑512733

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: