Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,455,692 Δ1 bp 59.1% coding (1210/1260 nt) lyrA → Lysostaphin resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,455,6920G.59.1% 114.6 / 213.3 122coding (1210/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base G (23/13);  new base . (35/37);  total (66/56)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 130 total aligned reads displayed.

GAATCAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCC  >  USA300TCH1516_ALE/2455554‑2455832
                                                                                                                                          |                                                                                                                                            
gAATCAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCaaaa                                                                                                                                            >  1:384229/1‑141 (MQ=255)
gAATCAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGAGTGAACAAATGAGTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCaaaa                                                                                                                                            >  2:879069/1‑141 (MQ=255)
gAATCAGCAGAATAAGTTAAACGAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCaaaa                                                                                                                                            >  2:258489/1‑141 (MQ=255)
      gcagAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaatta                                                                                                                                      <  1:1199384/141‑4 (MQ=255)
         gAATCAGTTAAACAAGCACCACAGAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcg                                                                                                                                   >  2:1017634/1‑135 (MQ=255)
         gAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcg                                                                                                                                   <  2:1006411/141‑7 (MQ=255)
          aaTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcga                                                                                                                                  >  2:458724/1‑134 (MQ=255)
          aaTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCg                                                                                                                                  >  2:1075839/1‑141 (MQ=255)
            tCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaa                                                                                                                               <  1:668845/141‑1 (MQ=255)
              aGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGaaaa                                                                                                                              >  2:382043/1‑141 (MQ=255)
              aGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                                                                             <  2:1051871/141‑1 (MQ=255)
               gTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATAGTCGaaaaaa                                                                                                                            >  2:1016628/1‑141 (MQ=255)
               aatAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACCTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaat                                                                                                                            <  2:320077/139‑2 (MQ=255)
                 aaaaCAAGCACCACAAAGTGCCGATATAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAAATGATGCAAAAAAATTATCGAAAAAACt                                                                                                                          <  2:1054325/140‑1 (MQ=255)
                  aaaCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTg                                                                                                                         <  2:976007/141‑1 (MQ=255)
                    aCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTga                                                                                                                        <  2:30077/141‑1 (MQ=255)
                       aGCACCAAAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAAAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                                                                    <  2:270473/141‑1 (MQ=255)
                             acaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                                                              >  1:210713/1‑141 (MQ=255)
                             acaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                                                              >  2:950414/1‑141 (MQ=255)
                                aaGTGACGATTTAACAAACGATTCAACTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGATGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                                                                           <  1:723515/141‑1 (MQ=255)
                                  gTGACGATGTAACGAACGATTCAAAGGAAGATGAAATAGAGGAATCATTAAAAGAACCTGCGAGTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                         <  2:809519/141‑1 (MQ=255)
                                   tGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                                                                        <  1:1202121/141‑1 (MQ=255)
                                    gACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAgca                                                                                                        >  2:626699/1‑139 (MQ=255)
                                     aagATATAACAAACGATTCAAATAAAGATGAAAAAAAGCAATCATTAAAAGAACCTGTGACTTATAAAGAAGACAGACTTTCATCAGTTGTAATTGATGCAAAAAAAAAAAAAAAAAAACTGAAGAACAATCTTCAGATaa                                                                                                      <  2:242265/139‑1 (MQ=255)
                                     aCGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTGCAGATaa                                                                                                      >  2:954802/1‑141 (MQ=255)
                                       gATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATaaa                                                                                                     <  1:100035/141‑1 (MQ=255)
                                       gATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                                    >  1:701066/1‑141 (MQ=255)
                                       gATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACCTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATaaa                                                                                                     >  2:298780/1‑141 (MQ=255)
                                                aaCGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataa                                                                                            >  1:47206/1‑141 (MQ=255)
                                                 aCGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGAAGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTGCAGATGAAAACAAataata                                                                                          >  1:382381/1‑141 (MQ=255)
                                                  agATTGAAATGAAGATAAAATAAAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAAAAAAAAAAAAAACTGAAGAACAATCTTCAGATAAAAACAAAtaataa                                                                                         <  2:925442/140‑1 (MQ=255)
                                                     atCAAATAAAGACAAAATAGAACAATAATTAAAAGAACCAGCGACTTATAAAGAAGACAGACGTACATCAGTTGAAATTGATACAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                                                                      <  2:1097263/140‑1 (MQ=255)
                                                      tCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                                                                      <  2:807197/141‑1 (MQ=255)
                                                        aaaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                                                                                   >  1:187264/1‑141 (MQ=255)
                                                         aaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcg                                                                                  >  2:919394/1‑141 (MQ=255)
                                                         aaTGAAGATGAAATAGAGCAACCATTAAAAGAACCTGCTCCTCATAAAGCAGACCGACAGACATCAGCACATTCGGATGTAAAAAAATCATCGAAAAAACTAAAGAACAATCTTCAGATAAAACCAAAAAAAAACaacacg                                                                                  >  2:400625/1‑137 (MQ=255)
                                                          aTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcg                                                                                  >  2:295897/1‑141 (MQ=255)
                                                             aaGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCtt                                                                               >  1:1029348/1‑141 (MQ=255)
                                                               gATTAAATAGAGCAATCATTAAAAGAACCTGCGACGTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAACAAACTGCAGAACAATCTTCAGATCAAACCAACTAATAACAAAGCGCTTGCt                                                                            <  2:699278/141‑1 (MQ=255)
                                                               gATCAAATAGAGCAATCATAAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCACCAGTTGAAATTGATGCAAAAAAAATAACGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCt                                                                            <  2:1095116/141‑1 (MQ=255)
                                                                 tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCATGCTac                                                                          >  2:969635/1‑140 (MQ=255)
                                                                  gAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAg                                                                          >  1:603892/1‑141 (MQ=255)
                                                                      aagagGAAACATTAAAAGAACCTGCGACTTACAAAGAAGAAAGAAGTTCATCAGTTGTAATTGATGCAGAAAAAAATAAAAAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAcc                                                                      <  2:27484/140‑1 (MQ=255)
                                                                       agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctc                                                                    >  2:649973/1‑141 (MQ=255)
                                                                       agaaCAAGCATTAAAAAAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGACAAAACTGAAGAACAAGTTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctc                                                                    <  2:492301/141‑1 (MQ=255)
                                                                        gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                                   <  2:178035/141‑1 (MQ=255)
                                                                        gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                                   <  1:898236/141‑1 (MQ=255)
                                                                          gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGCTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCtt                                                                  >  2:879240/1‑141 (MQ=255)
                                                                           cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                                                                >  1:686918/1‑141 (MQ=255)
                                                                            aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGCTAGTACCTCTTaa                                                                >  1:328332/1‑141 (MQ=255)
                                                                             agagatAAAAAAACCTGCGACTTAGAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAACGCGCTTGCTAGTACCTCTTaaaa                                                              <  2:434281/136‑1 (MQ=255)
                                                                             aTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa                                                               <  2:859502/141‑1 (MQ=255)
                                                                              tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTacaa                                                              >  1:407308/1‑141 (MQ=255)
                                                                               cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg                                                            <  1:1017634/141‑1 (MQ=255)
                                                                                 ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAga                                                           >  1:172532/1‑141 (MQ=255)
                                                                                   aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGGTGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAATCGCTTGCTAGTACCTCTTAAAAAgatga                                                        >  1:833231/1‑141 (MQ=255)
                                                                                   aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatga                                                        >  1:1161600/1‑141 (MQ=255)
                                                                                    aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatga                                                        >  1:217078/1‑141 (MQ=255)
                                                                                    aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatgat                                                       >  2:818868/1‑141 (MQ=255)
                                                                                     aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGatgatg                                                      >  1:355377/1‑141 (MQ=255)
                                                                                       gAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCt                                                    <  2:787664/141‑1 (MQ=255)
                                                                                        aaCCTGCGACTTAAAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCt                                                    <  1:298780/141‑1 (MQ=255)
                                                                                         aCCTGCGACTTATAAAGAAGACAGACGTTAACCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGAAGATGCTa                                                   >  1:284669/1‑141 (MQ=255)
                                                                                          ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGATAGc                                                 >  2:889020/1‑141 (MQ=255)
                                                                                          ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTCGCTAGTACATCTTAAAAAGATGATGCTaac                                                 >  2:342196/1‑139 (MQ=255)
                                                                                           cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                                <  1:551511/141‑1 (MQ=255)
                                                                                            tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCaa                                               <  2:213822/141‑1 (MQ=255)
                                                                                             gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAg                                              >  2:802257/1‑141 (MQ=255)
                                                                                             gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAg                                              <  2:904730/141‑1 (MQ=255)
                                                                                              cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgc                                             >  1:4143/1‑141 (MQ=255)
                                                                                                aCTTATAAAGACGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                           <  1:252142/141‑1 (MQ=255)
                                                                                                aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                           <  1:13458/141‑1 (MQ=255)
                                                                                                  ttATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                          <  2:407308/141‑1 (MQ=255)
                                                                                                    ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCatt                                        >  2:459468/1‑141 (MQ=255)
                                                                                                    ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtttt                                       >  1:370321/1‑141 (MQ=255)
                                                                                                    ataAAGAAGACAGACGTTCATCAGTTGAAAGTGATGCAGAAAAACATATCGAAAAAACTGAATAACAATCTTCAGATAAAAACAAATAAGAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCttt                                        <  2:1131154/141‑1 (MQ=255)
                                                                                                          aagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatac                                 <  1:474461/141‑1 (MQ=255)
                                                                                                           agacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                                >  1:696761/1‑141 (MQ=255)
                                                                                                              cagacGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactat                              <  2:845668/141‑1 (MQ=255)
                                                                                                              cagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactata                             >  2:106024/1‑141 (MQ=255)
                                                                                                              cagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatacgata                             >  1:954457/1‑141 (MQ=255)
                                                                                                                 acGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGGAGCTAGCAATCGATTTTctatactatttt                           >  1:61889/1‑138 (MQ=255)
                                                                                                                   gTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATAtatta                        >  1:967262/1‑141 (MQ=255)
                                                                                                                     tcatcaGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATAtattat                       >  1:509156/1‑141 (MQ=255)
                                                                                                                     tcatcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATattatt                      >  1:581586/1‑141 (MQ=255)
                                                                                                                     tcatcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACAATATattatt                      >  1:1163455/1‑141 (MQ=255)
                                                                                                                      catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATCCTATATATTAttt                     >  2:306880/1‑141 (MQ=255)
                                                                                                                      catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCGTGCTAGTACCTCTTAAAAAGATACTGCTCGCAAGCGCTTTTCTATACTATATATTAttt                     >  1:1069352/1‑141 (MQ=255)
                                                                                                                       atcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCCAGCGCTTTTCTATACTATATATTAtttt                    <  1:240033/141‑1 (MQ=255)
                                                                                                                       atcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAtttt                    <  1:919394/141‑1 (MQ=255)
                                                                                                                       atcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAAATGAAGAACAATCTTAAGATACAACCAAATAACCATCAACCGCCTGATCGAACCCCTTAAAAAGATGATGCTAACAAGCGCGTTTCTATACTATATATTAtttt                    >  2:694256/1‑141 (MQ=255)
                                                                                                                        tcaGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATattttatttt                    >  2:220004/1‑141 (MQ=255)
                                                                                                                        tcaGTTGTAATTGATGCAAAAAAATTATCGAAAGAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAttttt                   <  2:665533/141‑1 (MQ=255)
                                                                                                                          aGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTCTAATATATATTATTTTTc                  >  1:447077/1‑141 (MQ=255)
                                                                                                                          aGTTGTAATTGATGCAAAAAAATTATCGAAAAACCTGGAGAAGAATGTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCt                 <  1:1118386/141‑1 (MQ=255)
                                                                                                                            ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaa               >  2:619816/1‑141 (MQ=255)
                                                                                                                            ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaa               <  1:950414/141‑1 (MQ=255)
                                                                                                                                aaTTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTACAAt            >  1:344170/1‑141 (MQ=255)
                                                                                                                                aaTTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCGTTTCTATCCTATATATTATTGTTCTAAAAt            >  1:531429/1‑141 (MQ=255)
                                                                                                                                          gAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACAccc  >  1:487269/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
GAATCAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCC  >  USA300TCH1516_ALE/2455554‑2455832

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