Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,825,599 A→G 8.4% intergenic (‑414/+434) cap8A_2 ← / ← icaR Capsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,825,5990AG 8.4% 769.6 / 49.8 263intergenic (‑414/+434)cap8A_2/icaRCapsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR
Reads supporting (aligned to +/- strand):  ref base A (134/105);  new base G (20/2);  total (156/107)
Fisher's exact test for biased strand distribution p-value = 1.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

Only 100 of 269 total aligned reads displayed.

GAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGGCGCTTATCA  >  USA300TCH1516_ALE/2825509‑2825688
                                                                                          |                                                                                         
gAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGa                                                                                           >  2:3254929/1‑91 (MQ=255)
gAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGa                                                                                           <  1:725336/91‑1 (MQ=255)
 aGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAt                                                                                          <  2:245570/91‑1 (MQ=255)
   cGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgt                                                                                        <  1:2613687/91‑1 (MQ=255)
   cGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgt                                                                                        >  1:1490422/1‑91 (MQ=255)
     tGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTgg                                                                                      >  2:2893892/1‑91 (MQ=255)
      gCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTggg                                                                                     >  2:2147351/1‑91 (MQ=255)
       cccGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGt                                                                                    >  1:164655/1‑91 (MQ=255)
         cGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTaa                                                                                  >  2:2651155/1‑91 (MQ=255)
           tACTCGGGAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAAtt                                                                                <  2:3527228/91‑1 (MQ=255)
             cTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTta                                                                              >  1:3082301/1‑91 (MQ=255)
             cTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTta                                                                              >  1:2044239/1‑91 (MQ=255)
              tCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtat                                                                             <  1:1600685/91‑1 (MQ=255)
              tCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtat                                                                             <  2:833245/91‑1 (MQ=255)
               cGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtata                                                                            >  1:2714983/1‑91 (MQ=255)
                ggTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatat                                                                           <  1:1990638/91‑1 (MQ=255)
                    aTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAtt                                                                       >  1:2773726/1‑91 (MQ=255)
                    aTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAtt                                                                       <  2:2171898/91‑1 (MQ=255)
                      tCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAtttt                                                                     <  1:2113176/91‑1 (MQ=255)
                       cAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAttttt                                                                    >  1:63642/1‑91 (MQ=255)
                        aaTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTttta                                                                   <  1:3839/91‑1 (MQ=255)
                        aaTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTttta                                                                   <  1:2109281/91‑1 (MQ=255)
                         aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                                  >  2:1796817/1‑91 (MQ=255)
                         aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                                  >  1:457432/1‑91 (MQ=255)
                         aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                                  <  1:3177580/91‑1 (MQ=255)
                          tACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttatt                                                                 >  1:258259/1‑91 (MQ=255)
                          tACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttatt                                                                 <  1:3520280/91‑1 (MQ=255)
                           aCCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttattt                                                                <  2:919237/91‑1 (MQ=255)
                           aCCTTGCATAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttattt                                                                <  2:2590490/91‑1 (MQ=255)
                                 cGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAttt                                                          <  1:3299849/91‑1 (MQ=255)
                                   tAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAttttt                                                        <  2:1619908/91‑1 (MQ=255)
                                    aaCCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTTGGTAATTTATATAATTTTTATTTACATTTTTa                                                       >  1:3353900/1‑91 (MQ=255)
                                     aCCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTaa                                                      >  2:3374950/1‑91 (MQ=255)
                                      cccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTCTTTACATTTTTAAc                                                     >  1:1191242/1‑91 (MQ=255)
                                      cccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAAc                                                     >  2:2008543/1‑91 (MQ=255)
                                        cGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTACTTTATATAATTTTTATTTACATTTTTAACtt                                                   >  2:2230159/1‑91 (MQ=255)
                                         gTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTta                                                  >  1:1467124/1‑91 (MQ=255)
                                         gTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTta                                                  >  1:2097523/1‑91 (MQ=255)
                                          tCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtat                                                 >  1:1056965/1‑91 (MQ=255)
                                          tCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtat                                                 <  2:1164577/91‑1 (MQ=255)
                                          tCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtat                                                 <  1:174102/91‑1 (MQ=255)
                                           cACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtata                                                >  2:2110120/1‑91 (MQ=255)
                                           cACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtata                                                >  1:602191/1‑91 (MQ=255)
                                            aCTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaa                                               <  1:2434551/91‑1 (MQ=255)
                                               gtgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaaa                                            >  1:2278236/1‑91 (MQ=255)
                                                tgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAACATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaaaa                                           >  1:3428624/1‑91 (MQ=255)
                                                   gTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgc                                        >  2:1011476/1‑91 (MQ=255)
                                                   gTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgc                                        <  2:1335324/91‑1 (MQ=255)
                                                   gTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgc                                        >  1:3130902/1‑91 (MQ=255)
                                                    ttGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgcg                                       <  1:277914/91‑1 (MQ=255)
                                                    ttGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGAGTAATTTATATAATTTTTATTTACATTTTTAAATTATAAAAAAAgcg                                       <  2:3914255/91‑1 (MQ=255)
                                                     tGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgcgc                                      >  2:666959/1‑91 (MQ=255)
                                                      gggTTAATGATAATAAAGCCCACACCTTTTGAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGcc                                     <  1:1276663/91‑1 (MQ=255)
                                                       ggTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCt                                    >  2:3412691/1‑91 (MQ=255)
                                                         ttAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTAt                                  >  1:160534/1‑91 (MQ=255)
                                                         ttAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTAt                                  >  2:853822/1‑91 (MQ=255)
                                                           aaTGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGt                                <  2:1037948/91‑1 (MQ=255)
                                                            aggATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAAATAATAAAAAAAGCGCCTATGTc                               <  2:3102977/89‑1 (MQ=255)
                                                            aTGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTTTTTACATTTTTAACTTATAAAAAAAGCGCCTATGTc                               >  2:1209918/1‑91 (MQ=255)
                                                            aTGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTc                               <  1:396253/91‑1 (MQ=255)
                                                             tGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCa                              <  1:2752431/91‑1 (MQ=255)
                                                             tGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCa                              >  1:3748791/1‑91 (MQ=255)
                                                               ataataAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATg                            >  2:2080940/1‑91 (MQ=255)
                                                                tcataaTGCCTACACCTTTTATAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGa                           <  2:2775896/89‑1 (MQ=255)
                                                                taataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGa                           <  1:1169004/91‑1 (MQ=255)
                                                                 aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATCTAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt                          >  2:3249411/1‑91 (MQ=255)
                                                                 aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt                          >  2:647699/1‑91 (MQ=255)
                                                                 aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt                          >  1:820154/1‑91 (MQ=255)
                                                                 aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt                          <  1:3931420/91‑1 (MQ=255)
                                                                   taaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAttt                        <  1:2515963/91‑1 (MQ=255)
                                                                   taaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAttt                        <  1:2892420/91‑1 (MQ=255)
                                                                    aaaGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTa                       <  2:2633674/91‑1 (MQ=255)
                                                                       gCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCa                    >  2:83591/1‑91 (MQ=255)
                                                                       gCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCa                    >  2:2308834/1‑91 (MQ=255)
                                                                        cccACATCTTTTTAAAAGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCAt                   >  2:1308896/1‑91 (MQ=39)
                                                                        cccACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCAt                   <  2:3706402/91‑1 (MQ=255)
                                                                        cccACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCAt                   <  1:1716107/91‑1 (MQ=255)
                                                                          cacaCCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATca                 >  1:235409/1‑91 (MQ=255)
                                                                           acacCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATcac                <  2:2864701/91‑1 (MQ=255)
                                                                             acCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACAt              <  1:956541/91‑1 (MQ=255)
                                                                             acCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACAt              >  1:2533061/1‑91 (MQ=255)
                                                                              ccTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATa             >  1:811257/1‑91 (MQ=255)
                                                                               cTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAg            >  1:2185538/1‑91 (MQ=255)
                                                                               cTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAg            >  1:1419131/1‑91 (MQ=255)
                                                                               cTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAg            >  1:1129621/1‑91 (MQ=255)
                                                                                ttttTAAAAGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCGATGTCATGATTTACCATCACATAgg           >  1:3670410/1‑91 (MQ=255)
                                                                                ttttAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAgg           <  1:3160908/91‑1 (MQ=255)
                                                                                 ttttaaaaGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGgc          >  2:3802263/5‑91 (MQ=255)
                                                                                 ttttaaaaGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGgc          >  1:1094871/5‑91 (MQ=255)
                                                                                 ttttaaaaGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGgc          >  1:1883184/5‑91 (MQ=255)
                                                                                 ttttAAAAGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCCTGATTTACCATCACATAGgc          >  2:186409/1‑91 (MQ=255)
                                                                                 tttAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGgc          >  1:2069971/1‑91 (MQ=255)
                                                                                  tttaaaaGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGgcg         >  2:541936/4‑91 (MQ=255)
                                                                                  tttaaaaGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGgcg         >  1:3634994/4‑91 (MQ=255)
                                                                                   tAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGgcgc        <  1:931001/91‑1 (MQ=255)
                                                                                    aaaaaGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGGCGCt       >  2:387796/1‑91 (MQ=255)
                                                                                      aaggATGTGAGTAATTTAGATAATTTTTATTTACATGTTTAACTTGTAAGAAGAGCGCCGATGTGATGGTAAATCATGACATAGGCGCttt     <  2:966011/88‑2 (MQ=9)
                                                                                      aaaGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTGATGGTAAATCATGACATAGGCGCttt     <  1:2980974/91‑2 (MQ=21)
                                                                                       aaGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGGCGCTTAt    >  1:386831/1‑91 (MQ=255)
                                                                                         taTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGGCGCTTATCa  <  2:1643627/90‑1 (MQ=255)
                                                                                          |                                                                                         
GAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATAGGCGCTTATCA  >  USA300TCH1516_ALE/2825509‑2825688

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: