Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 779,970 G→A 5.4% intergenic (+26/+46) csbB → / ← saeS Putative glycosyltransferase CsbB/Histidine protein kinase SaeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE779,9700GA 5.4% 781.0 / 18.5 238intergenic (+26/+46)csbB/saeSPutative glycosyltransferase CsbB/Histidine protein kinase SaeS
Reads supporting (aligned to +/- strand):  ref base G (118/105);  new base A (7/7);  total (126/112)
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 260 total aligned reads displayed.

AAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAGT  >  USA300TCH1516_ALE/779880‑780060
                                                                                          |                                                                                          
 aGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATggg                                                                                           >  1:2933273/1‑91 (MQ=255)
 aGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATggg                                                                                           >  1:1499086/1‑91 (MQ=255)
  gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGTCCTCATTAATGGGa                                                                                          <  1:2550226/91‑1 (MQ=255)
  gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                          <  1:3140035/91‑1 (MQ=255)
  gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                          >  1:1393527/1‑91 (MQ=255)
    tACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGc                                                                                        <  2:909772/91‑1 (MQ=255)
      cAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGctt                                                                                      >  2:552092/1‑91 (MQ=255)
      cAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGctt                                                                                      >  2:3513117/1‑91 (MQ=255)
       aaaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttc                                                                                     >  1:1548197/1‑91 (MQ=255)
        aaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                                    >  2:43709/1‑91 (MQ=255)
        aaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                                    >  2:2530372/1‑91 (MQ=255)
          tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCttt                                                                                  <  2:1499086/91‑1 (MQ=255)
          tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCttt                                                                                  <  1:2464767/91‑1 (MQ=255)
          tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCttt                                                                                  <  1:20526/91‑1 (MQ=255)
          tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCttt                                                                                  <  1:2096834/91‑1 (MQ=255)
           cTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCtttt                                                                                 >  2:1922549/1‑91 (MQ=255)
            tttCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGCTTCttttt                                                                                <  1:3823267/91‑1 (MQ=37)
             ttCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTa                                                                               >  2:846954/1‑91 (MQ=255)
              tCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAg                                                                              >  1:3375898/1‑91 (MQ=255)
              tCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAg                                                                              >  2:1637152/1‑91 (MQ=255)
               cAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATGCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGt                                                                             >  2:2951490/1‑91 (MQ=255)
                aagacTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTc                                                                            >  2:1833123/5‑91 (MQ=255)
                  gTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCtt                                                                          >  1:686556/1‑91 (MQ=255)
                  gTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCtt                                                                          >  1:2937289/1‑91 (MQ=255)
                   tCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCttt                                                                         <  2:1751153/91‑1 (MQ=255)
                    cTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTg                                                                        <  2:1199145/91‑1 (MQ=255)
                     ttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGc                                                                       <  1:1889280/91‑1 (MQ=255)
                     ttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGc                                                                       <  2:1029227/91‑1 (MQ=255)
                     ttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGc                                                                       <  2:1875118/91‑1 (MQ=255)
                      tGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                                      <  2:1548197/91‑1 (MQ=255)
                      tGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                                      >  1:3421153/1‑91 (MQ=255)
                      tGAAAATGATGAGAAGGATACCCATAAAGTTTATTATAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                                      <  2:106423/91‑1 (MQ=255)
                       gAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                                     <  1:1319776/91‑1 (MQ=255)
                       gAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                                     <  1:1726858/91‑1 (MQ=255)
                       gAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                                     >  1:1878225/1‑91 (MQ=255)
                         aaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattt                                                                   >  1:873402/1‑91 (MQ=255)
                           atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                 >  2:3064143/1‑91 (MQ=255)
                           atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                 >  2:965797/1‑91 (MQ=255)
                           atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                 >  1:2724983/1‑91 (MQ=255)
                           atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                 >  1:2902694/1‑91 (MQ=255)
                           atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                 <  1:3007132/91‑1 (MQ=255)
                           atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                 >  1:3323943/1‑91 (MQ=255)
                            tgatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCttttt                                                                                <  2:2831335/75‑1 (MQ=255)
                              atgaGAAGGATACCCATAAAGTTTATACTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttattt                                                              >  2:2218478/1‑91 (MQ=255)
                                gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                                            <  1:122990/91‑1 (MQ=255)
                                gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                                            <  1:2178652/91‑1 (MQ=255)
                                  gaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTtata                                                          >  1:27670/1‑91 (MQ=255)
                                     ggATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaat                                                       <  2:1984095/91‑1 (MQ=255)
                                      gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTataaaata                                                     <  1:1864879/91‑4 (MQ=255)
                                        tACCCATAAAGTTTATTCGAAATAAGCAAAAAAAGAAGCCCTCATTAATGGGAGTTTCTGTTTAGTCTTTGCATTTTATTGTAtaaagaaa                                                    <  2:2418145/91‑1 (MQ=255)
                                         aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                                   <  1:3064066/91‑1 (MQ=255)
                                          cccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATc                                                  <  2:2683286/91‑1 (MQ=255)
                                                aaGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTa                                            <  2:923621/91‑1 (MQ=255)
                                                 aGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTAt                                           <  1:835697/91‑1 (MQ=255)
                                                  gTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATg                                          >  2:3819541/1‑91 (MQ=255)
                                                  gTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATg                                          <  2:2166521/91‑1 (MQ=255)
                                                  gTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATg                                          >  2:2108866/1‑91 (MQ=255)
                                                      aTTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGt                                      <  2:3154204/91‑1 (MQ=255)
                                                         cTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAAt                                   <  1:2808079/91‑1 (MQ=255)
                                                            aataaaCGAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTc                                <  1:1943599/91‑1 (MQ=255)
                                                               aaaCAAAAAAAGAAGCTCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTaa                             >  2:1312132/1‑91 (MQ=255)
                                                               aaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTaa                             <  2:2326263/91‑1 (MQ=255)
                                                                 aCAAAAAAAGAAGCCCTCATTAATGGGAACTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAtt                           <  2:2254518/91‑1 (MQ=255)
                                                                  cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt                          <  2:264752/91‑1 (MQ=255)
                                                                    aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt                        <  2:3511659/91‑1 (MQ=255)
                                                                    aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt                        >  1:3567349/1‑91 (MQ=255)
                                                                    aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt                        >  1:413642/1‑91 (MQ=255)
                                                                     aaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgtg                       >  2:2475591/1‑91 (MQ=255)
                                                                      aaaaGAAGCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTgtgt                      >  2:3298841/1‑91 (MQ=37)
                                                                        aagaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAAt                                   >  2:2400889/1‑76 (MQ=255)
                                                                         agaagCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAAt                   >  1:3764452/1‑91 (MQ=37)
                                                                          gaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATg                  >  1:2344717/1‑91 (MQ=255)
                                                                          gaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATg                  >  2:3125483/1‑91 (MQ=255)
                                                                          gaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATg                  <  2:2994325/91‑1 (MQ=255)
                                                                           aagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGt                 <  1:10050/91‑1 (MQ=255)
                                                                           aagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGt                 >  1:305601/1‑91 (MQ=255)
                                                                           aagCCCTCATTAATGGGAGCTTCTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                <  2:243861/91‑1 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                >  2:111943/1‑91 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                >  2:720222/1‑91 (MQ=255)
                                                                              cccTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTAc              <  1:945748/91‑1 (MQ=255)
                                                                               ccTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACa             >  1:640825/1‑91 (MQ=255)
                                                                                cccATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAg            >  1:1163132/3‑91 (MQ=255)
                                                                                  cATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTc          <  2:1407907/91‑1 (MQ=255)
                                                                                  cATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTc          <  2:1776127/91‑1 (MQ=255)
                                                                                   aTTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCa         <  2:3039870/91‑1 (MQ=255)
                                                                                   aTTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCa         <  1:3681769/91‑1 (MQ=255)
                                                                                    ttAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCAt        >  1:3515351/1‑91 (MQ=255)
                                                                                    ttAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCAt        >  1:3514846/1‑91 (MQ=255)
                                                                                    ttAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCAt        <  1:2366981/91‑1 (MQ=255)
                                                                                     tAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATc       <  2:3682982/91‑1 (MQ=255)
                                                                                      aaTGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCg      >  1:1600611/1‑91 (MQ=255)
                                                                                      aaTGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCg      <  1:1082853/91‑1 (MQ=255)
                                                                                      aaTGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTATAATGTTACAGTCATCg      >  1:3326905/1‑91 (MQ=255)
                                                                                        tGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTa    <  1:1272010/91‑1 (MQ=255)
                                                                                         gggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAg   <  2:338465/91‑1 (MQ=255)
                                                                                         gggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAg   >  2:3183669/1‑91 (MQ=255)
                                                                                         gggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAg   <  1:3465202/91‑1 (MQ=255)
                                                                                         gggAGCTTCTTTTTAGTCTTTGCATGTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAg   <  1:3629010/91‑1 (MQ=255)
                                                                                          ggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAGt  <  2:371959/91‑1 (MQ=255)
                                                                                          ggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAGt  >  2:818654/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAGT  >  USA300TCH1516_ALE/779880‑780060

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