Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 779,970 G→A 6.6% 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 6.6% 854.4 / 30.3 265intergenic (+26/+46)csbB/saeSPutative glycosyltransferase CsbB/Histidine protein kinase SaeS
Reads supporting (aligned to +/- strand):  ref base G (113/134);  new base A (10/8);  total (123/142)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

AAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAGT  >  USA300TCH1516_ALE/779880‑780060
                                                                                          |                                                                                          
aaGCTACAAATGTTTCAAGTCTTGAAAATAATGAGAAATATAGCAATAAAGTTTATTTTAAATATATAAAAAAAGAAGCCCTCATTAATgg                                                                                            <  2:2860536/91‑1 (MQ=255)
  gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                          <  1:120228/91‑1 (MQ=255)
  gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                                          <  1:1502373/91‑1 (MQ=255)
      cAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGctt                                                                                      >  2:577957/1‑91 (MQ=255)
       aaaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGcttc                                                                                     <  1:2956082/91‑1 (MQ=37)
        aaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                                    >  2:3345977/1‑91 (MQ=255)
        aaTCTTTCAAGTCTTGAAAATGATGAGAAGAATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                                    <  2:1204860/91‑1 (MQ=255)
         aTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttctt                                                                                   >  1:1333465/1‑91 (MQ=255)
          tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCttt                                                                                  <  1:438897/91‑1 (MQ=255)
              tCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAg                                                                              <  2:979108/91‑1 (MQ=255)
                aaGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTc                                                                            <  1:3535689/91‑1 (MQ=255)
                 aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                                           <  2:904138/91‑1 (MQ=255)
                 aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                                           >  2:3643360/1‑91 (MQ=255)
                  gTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCtt                                                                          <  1:567482/91‑1 (MQ=255)
                   tCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCttt                                                                         <  2:3166192/91‑1 (MQ=255)
                   tCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCttt                                                                         >  1:664263/1‑91 (MQ=255)
                    gttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTg                                                                        <  2:3921471/90‑1 (MQ=255)
                    cTTCAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTg                                                                        <  1:3061746/91‑1 (MQ=255)
                     ttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGc                                                                       <  1:1710587/91‑1 (MQ=255)
                     ttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGc                                                                       >  1:1289620/1‑91 (MQ=255)
                       gAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                                     <  2:2818436/91‑1 (MQ=255)
                       gAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                                     <  1:2335105/91‑1 (MQ=255)
                       gAAAATGATGAGAAGGATACCCATAAAGTTTATTATAAATAAAAAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                                     <  2:3261791/91‑1 (MQ=255)
                        aaaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatt                                                                    <  2:2573264/91‑1 (MQ=255)
                         aaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattt                                                                   <  2:2677740/91‑1 (MQ=255)
                           atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                                 >  2:1876710/1‑91 (MQ=255)
                            tgatgaGAAGGATACCCATAAAGTTTATTCTAGATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttat                                                                <  2:2539525/91‑1 (MQ=255)
                            tgatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttat                                                                >  1:3320121/1‑91 (MQ=255)
                            tgatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttat                                                                <  2:2539037/91‑1 (MQ=255)
                              atgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttattt                                                              <  2:3054198/91‑1 (MQ=255)
                              atgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttattt                                                              >  2:3278683/1‑91 (MQ=255)
                               tgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttatttt                                                             >  1:3391899/1‑91 (MQ=255)
                                gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                                            >  2:1263300/1‑91 (MQ=255)
                                gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                                            >  1:2046004/1‑91 (MQ=255)
                                 agaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATtttattttat                                                           <  2:1291922/91‑1 (MQ=255)
                                 agaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATtttattttat                                                           <  2:1292244/91‑1 (MQ=255)
                                  gaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttatttataa                                                          <  2:2194705/91‑5 (MQ=255)
                                  gaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTtata                                                          >  1:179280/1‑91 (MQ=255)
                                  gaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTtata                                                          <  1:852400/91‑1 (MQ=255)
                                   aaGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataa                                                         <  1:1962066/91‑1 (MQ=255)
                                    aGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaa                                                        <  2:1115247/91‑1 (MQ=255)
                                     ggATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaat                                                       <  1:48011/91‑1 (MQ=255)
                                      gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTattaatg                                                      >  2:1854017/1‑90 (MQ=255)
                                      gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaata                                                      <  1:3417237/91‑1 (MQ=255)
                                      gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaata                                                      >  1:1325796/1‑91 (MQ=255)
                                      gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTGAGTCTTTGCATTTTATTTTataaata                                                      <  2:2234821/91‑1 (MQ=255)
                                         aCCCATAGAGTTGATTGTAAATGAACAAAAAAAGAAGCCCTCATTAATGGGTGCTTCTTTTTAGTCTTTGTATTTTATTTTATaaataaat                                                   <  2:1174581/91‑1 (MQ=255)
                                          cccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATc                                                  >  2:77447/1‑91 (MQ=255)
                                          cccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATc                                                  <  2:356583/91‑1 (MQ=255)
                                          cccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATc                                                  <  1:1506017/91‑1 (MQ=255)
                                           ccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCg                                                 <  2:3729416/91‑1 (MQ=255)
                                            cATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCgg                                                <  1:289596/91‑1 (MQ=255)
                                            cATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCgg                                                <  1:1872315/91‑1 (MQ=255)
                                                 aGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTAt                                           <  2:182733/91‑1 (MQ=255)
                                                 aGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTACAAATAAATCGGATTAt                                           >  2:3789413/1‑91 (MQ=255)
                                                  gTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATg                                          <  1:214432/91‑1 (MQ=255)
                                                    ttATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGAc                                        <  2:85917/91‑1 (MQ=255)
                                                    ttATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGAc                                        >  1:1597962/1‑91 (MQ=255)
                                                    ttATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGAc                                        <  1:1081643/91‑1 (MQ=255)
                                                      aTTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGt                                      >  1:1739243/1‑91 (MQ=255)
                                                      aTTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGt                                      >  2:108062/1‑91 (MQ=255)
                                                         ttaaataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAAt                                   <  2:3926741/90‑1 (MQ=255)
                                                            aataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTc                                >  1:246403/1‑91 (MQ=255)
                                                             ataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCt                               <  1:933332/91‑1 (MQ=255)
                                                             ataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCt                               <  2:2525999/91‑1 (MQ=255)
                                                               aaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTaa                             >  1:2774213/1‑91 (MQ=255)
                                                                aaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAt                            <  1:1114257/91‑1 (MQ=255)
                                                                aaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAt                            <  1:3324652/91‑1 (MQ=255)
                                                                aaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAt                            >  2:2344216/1‑91 (MQ=255)
                                                                 aCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAtt                           >  2:1507440/1‑91 (MQ=255)
                                                                  cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt                          <  1:1744194/91‑1 (MQ=255)
                                                                  cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt                          >  2:2021976/1‑91 (MQ=255)
                                                                   aaaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtg                         <  1:296679/91‑1 (MQ=255)
                                                                    aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt                        <  1:2661255/91‑1 (MQ=255)
                                                                     aaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgtg                       >  2:1545327/1‑91 (MQ=255)
                                                                     aaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgtg                       >  2:225078/1‑91 (MQ=255)
                                                                        aagaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTaa                    >  1:2482566/1‑91 (MQ=255)
                                                                         agaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAAt                   >  1:1394927/1‑91 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                >  1:732833/1‑91 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                >  1:1787449/1‑91 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                >  2:2636068/1‑91 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                >  2:68885/1‑91 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                <  2:206503/91‑1 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt                >  2:1377524/1‑91 (MQ=255)
                                                                            agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGCAATGtt                <  1:3078364/91‑1 (MQ=255)
                                                                              cccTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTAc              <  2:1417443/91‑1 (MQ=255)
                                                                               tcccATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACa             >  2:3002687/4‑91 (MQ=39)
                                                                               tcccATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACa             >  2:1888424/4‑91 (MQ=39)
                                                                                cccATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAg            >  2:3210435/3‑91 (MQ=255)
                                                                                 tCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATTAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGt           <  2:3046620/91‑1 (MQ=255)
                                                                                 ccATTAATGAGGGCTTCTCTCTAGTCTTTGCACTTTACTGTCTAAATCAATCGGATTATGACGTAATGTCTAGTTTGTGTAATGTCACAGt           >  1:1432138/2‑91 (MQ=255)
                                                                                  cATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTc          >  1:3827932/1‑91 (MQ=255)
                                                                                   aTTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCa         <  1:2632026/91‑1 (MQ=255)
                                                                                   aTTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCa         >  1:3622088/1‑91 (MQ=255)
                                                                                    ttAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCAt        <  2:389486/91‑1 (MQ=255)
                                                                                    ttAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCAt        <  2:474920/91‑1 (MQ=255)
                                                                                        tGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTa    <  1:3896063/91‑1 (MQ=255)
                                                                                         gggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAg   <  2:1929848/91‑1 (MQ=255)
                                                                                         gggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAg   <  2:25174/91‑1 (MQ=255)
                                                                                          ggAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAGt  >  1:2823465/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTACAGTCATCGTAGT  >  USA300TCH1516_ALE/779880‑780060

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