Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 779,970 G→A 10.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,9700GA10.4% 402.6 / 28.8 135intergenic (+26/+46)csbB/saeSPutative glycosyltransferase CsbB/Histidine protein kinase SaeS
Reads supporting (aligned to +/- strand):  ref base G (63/58);  new base A (8/6);  total (71/64)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 137 total aligned reads displayed.

AGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAG  >  USA300TCH1516_ALE/779881‑780059
                                                                                         |                                                                                         
 gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCTATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                         <  1:1538125/91‑1 (MQ=255)
 gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                         >  1:3871866/1‑91 (MQ=255)
 gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                         <  1:3640773/91‑1 (MQ=255)
 gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                         <  1:2788662/91‑1 (MQ=255)
 gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                         <  1:2015244/91‑1 (MQ=255)
  cTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAg                                                                                        <  1:2939254/91‑1 (MQ=255)
   tACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGc                                                                                       <  1:366500/91‑1 (MQ=255)
   tACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGc                                                                                       <  1:2489437/91‑1 (MQ=255)
   tACAAATCTCTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGc                                                                                       <  1:3967846/91‑1 (MQ=255)
    aCAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGct                                                                                      >  1:526692/1‑91 (MQ=255)
     cAAATCTTTCAAGTCTTGAACATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGctt                                                                                     <  1:2522786/91‑1 (MQ=255)
      aaaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGcttc                                                                                    <  1:1111664/91‑1 (MQ=37)
      aaaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGcttc                                                                                    <  1:2917143/91‑1 (MQ=37)
      aaaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGcttc                                                                                    <  1:587875/91‑1 (MQ=38)
       aaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                                   >  1:2363487/1‑91 (MQ=255)
       aaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                                   <  1:3832468/91‑1 (MQ=255)
        aTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttctt                                                                                  >  1:3625904/1‑91 (MQ=255)
        aTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttctt                                                                                  >  1:2910690/1‑91 (MQ=255)
         tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGCTTCttt                                                                                 <  1:1291955/91‑1 (MQ=37)
         tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCttt                                                                                 <  1:2973559/91‑1 (MQ=255)
          cTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGCTTCtttt                                                                                <  1:3983295/91‑1 (MQ=37)
               aaGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTc                                                                           <  1:2454060/91‑1 (MQ=255)
                aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                                          >  1:1601521/1‑91 (MQ=255)
                aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                                          >  1:1399328/1‑91 (MQ=255)
                aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                                          >  1:1952273/1‑91 (MQ=255)
                 gTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAAGAAACAAAAAAAGAAGCCCTCATTAAGGGGAGCTTCTTTTTAGTCtt                                                                         <  1:3296817/91‑1 (MQ=255)
                    ttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGc                                                                      >  1:373344/1‑91 (MQ=255)
                     tGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                                     <  1:1982861/91‑1 (MQ=255)
                     tGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                                     <  1:3566517/91‑1 (MQ=255)
                      gAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                                    <  1:1853026/91‑1 (MQ=255)
                       aaaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGCAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatt                                                                   >  1:75437/1‑91 (MQ=255)
                       aaaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatt                                                                   <  1:3432025/91‑1 (MQ=255)
                        aaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattt                                                                  >  1:837283/1‑91 (MQ=255)
                         aaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttt                                                                 <  1:1670276/91‑1 (MQ=255)
                         aaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttt                                                                 >  1:1990662/1‑91 (MQ=255)
                          atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                >  1:3549546/1‑91 (MQ=255)
                          atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                >  1:3549496/1‑91 (MQ=255)
                           tgatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttat                                                               >  1:1552388/1‑91 (MQ=255)
                            gatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttatt                                                              >  1:362975/1‑91 (MQ=255)
                             atgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttattt                                                             <  1:2820658/91‑1 (MQ=255)
                             atgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttattt                                                             >  1:2841204/1‑91 (MQ=255)
                              tgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttatttt                                                            >  1:2005729/1‑91 (MQ=255)
                               gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                                           >  1:2040908/1‑91 (MQ=255)
                               gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                                           >  1:1137533/1‑91 (MQ=255)
                                agaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATtttattttat                                                          >  1:658925/1‑91 (MQ=255)
                                agaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATtttattttat                                                          >  1:659030/1‑91 (MQ=255)
                                 gaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTtata                                                         <  1:2858168/91‑1 (MQ=255)
                                  aaGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataa                                                        >  1:3714605/1‑91 (MQ=255)
                                     gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaata                                                     >  1:914648/1‑91 (MQ=255)
                                     gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaata                                                     >  1:527508/1‑91 (MQ=255)
                                     gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaata                                                     >  1:527554/1‑91 (MQ=255)
                                      aTACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaataa                                                    <  1:843677/91‑1 (MQ=255)
                                       tACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTAtaaataaa                                                   <  1:2480988/91‑1 (MQ=255)
                                        aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                                  <  1:584665/91‑1 (MQ=255)
                                        aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                                  <  1:1986578/91‑1 (MQ=255)
                                        aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                                  <  1:3932510/91‑1 (MQ=255)
                                        aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                                  <  1:583243/91‑1 (MQ=255)
                                          ccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCg                                                <  1:1581721/91‑1 (MQ=255)
                                          ccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCg                                                <  1:1492247/91‑1 (MQ=255)
                                           cATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCgg                                               >  1:585906/1‑91 (MQ=255)
                                                 gTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATg                                         <  1:856498/91‑1 (MQ=255)
                                                    tATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACg                                      >  1:1641550/1‑91 (MQ=255)
                                                      ttCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTa                                    >  1:364147/1‑91 (MQ=255)
                                                         taaataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATg                                 >  1:1929003/1‑91 (MQ=255)
                                                          aaataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGt                                <  1:2912304/91‑1 (MQ=255)
                                                               aaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAt                           >  1:2175011/1‑91 (MQ=255)
                                                                 cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt                         >  1:1963111/1‑91 (MQ=255)
                                                                 cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGCCTCTCCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt                         >  1:475733/1‑91 (MQ=255)
                                                                  aaaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtg                        >  1:2927139/1‑91 (MQ=255)
                                                                   aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt                       <  1:2527830/91‑1 (MQ=255)
                                                                   aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt                       <  1:2526866/91‑1 (MQ=255)
                                                                     aaaaGAAGCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTgtgt                     >  1:45642/1‑91 (MQ=37)
                                                                     aaaaGAAGCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTgtgt                     >  1:2667524/1‑91 (MQ=37)
                                                                      aaaGAAGCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTa                    >  1:2494021/1‑91 (MQ=37)
                                                                        agaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAAt                  >  1:929841/1‑91 (MQ=255)
                                                                         gaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTACGACGTAATGTCTAATTTGTGTAATg                 >  1:3847150/1‑91 (MQ=255)
                                                                          aagCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGt                >  1:2320111/1‑91 (MQ=38)
                                                                          aagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGt                <  1:2331382/91‑1 (MQ=255)
                                                                           agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt               >  1:1478149/1‑91 (MQ=255)
                                                                           agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt               >  1:702688/1‑91 (MQ=255)
                                                                           agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt               >  1:685368/1‑91 (MQ=255)
                                                                           agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt               >  1:2870800/1‑91 (MQ=255)
                                                                           agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt               <  1:2602243/91‑1 (MQ=255)
                                                                           agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt               >  1:2595843/1‑91 (MQ=255)
                                                                            gCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTa              <  1:1335369/91‑1 (MQ=255)
                                                                            gCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTa              >  1:3154343/1‑91 (MQ=255)
                                                                            gCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTa              <  1:1334562/91‑1 (MQ=255)
                                                                              tcccATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACa            >  1:3818696/4‑91 (MQ=39)
                                                                               cccATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAg           >  1:3252107/3‑91 (MQ=255)
                                                                                  aTTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCa        <  1:526130/91‑1 (MQ=255)
                                                                                  aTTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCa        <  1:458752/91‑1 (MQ=255)
                                                                                   ttAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCAt       <  1:3325323/91‑1 (MQ=255)
                                                                                   ttAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCAt       >  1:411564/1‑91 (MQ=255)
                                                                                    tAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATc      >  1:1439262/1‑91 (MQ=255)
                                                                                    tAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATc      >  1:1439506/1‑91 (MQ=255)
                                                                                     aaTGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCg     <  1:1125506/91‑1 (MQ=255)
                                                                                      aTGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGt    >  1:2205011/1‑91 (MQ=255)
                                                                                      aTGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGt    >  1:2509998/1‑91 (MQ=255)
                                                                                      aTGGGAGCTGCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGt    <  1:3720563/91‑1 (MQ=255)
                                                                                        gggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAg  <  1:1858017/91‑1 (MQ=255)
                                                                                         |                                                                                         
AGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAG  >  USA300TCH1516_ALE/779881‑780059

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: