Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,455,692 Δ1 bp 100% 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.76.0% 183.1 / 65.1 75coding (1210/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base . (29/28);  minor base A (11/2);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 3.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

AGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTT  >  USA300TCH1516_ALE/2455559‑2455824
                                                                                                                                     |                                                                                                                                     
agcatAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAAGGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaatt                                                                                                                                <  1:263058/141‑3 (MQ=255)
agcagAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaatt                                                                                                                                <  1:792002/141‑3 (MQ=255)
    gAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcg                                                                                                                            <  1:490024/141‑7 (MQ=255)
    gAATCAGTTAAACAAGCACCACAAAGTGACGATTGAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcg                                                                                                                            >  2:592686/1‑135 (MQ=255)
    gAAACAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcg                                                                                                                            <  1:284247/141‑7 (MQ=255)
    aaaTCAGTTAAACAAGCACCACAAAGTTACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcg                                                                                                                            <  1:382362/140‑7 (MQ=255)
     aaTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattatcga                                                                                                                           >  1:692060/1‑134 (MQ=255)
       tCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaa                                                                                                                       <  1:761292/141‑1 (MQ=255)
         aGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAACAAAACTATCGaaaaa                                                                                                                     >  1:714826/1‑141 (MQ=255)
             aaaCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGACGACAGACGTTCATCAGTTGAAATTCTGGCGAATAAATTATCGAAAAAACTg                                                                                                                 <  2:593979/141‑1 (MQ=255)
             aaaCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTg                                                                                                                 >  1:807011/1‑141 (MQ=255)
             aaaCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAAACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTg                                                                                                                 >  2:70219/1‑141 (MQ=255)
              aaCAAGCACCACAAAGCGTCGTTGTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGATTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTga                                                                                                                <  1:491548/141‑1 (MQ=255)
                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                                                                              >  2:832887/1‑141 (MQ=255)
                 aaGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAACAGAGCAAGCATTAAAAGAACCTGCGACTTATAGAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                                                                             <  1:618632/141‑1 (MQ=255)
                  aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                                                            >  1:205754/1‑141 (MQ=255)
                   acaccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCACTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATTCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                                                           <  2:94948/140‑1 (MQ=255)
                        aaaaaGTTACGAATTACCTAACGAATAAAAATAAGATGAAATAGAGCAATTATTAAAAGAAATTGCGACTTAAAAAGAATACAGACGTTCATCAGTTGTAAAAAAAAAAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                                                      <  2:653521/139‑1 (MQ=255)
                          aaaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                                                                    <  1:656697/141‑1 (MQ=255)
                            atcgACGATTCAACAAACAAAACAAACAAAGATGAAAGAGCGCAATCATTAAAAGAACCTACGACTTAAAAAAAATACAGACGTTCATCAGTTGAAAAAAAAAAAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                  <  2:767694/138‑1 (MQ=255)
                                  gATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                            >  2:833793/1‑141 (MQ=255)
                                  gATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                            >  2:601247/1‑141 (MQ=255)
                                   aTTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGCAGAATGGCATCTTTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaaa                                                                                           >  1:823350/1‑141 (MQ=255)
                                   aTTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaaa                                                                                           >  1:182089/1‑141 (MQ=255)
                                        aCAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACgaaa                                                                                      >  2:589076/1‑140 (MQ=255)
                                           aaCGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                   <  1:674119/141‑1 (MQ=255)
                                              gATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAATTATCGAAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAAc                                                                                <  1:810622/141‑1 (MQ=255)
                                                   aaaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCACAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                                                                           >  2:648872/1‑141 (MQ=255)
                                                    aaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACCAAgcg                                                                          >  1:82026/1‑141 (MQ=255)
                                                    aaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcg                                                                          <  1:110022/141‑1 (MQ=255)
                                                     aTGAAGATGAAATAGAGCAATCATTACAAGAACCTGCGACTTATAAAGAAGACAGACGATCATCAGTTGTAATTGATGCAAGAAAATTATAGAAAAAACTGAGGAAAAAACTTCAGATAAAAACCAAAAATAACAAAgccc                                                                         >  2:295378/1‑139 (MQ=255)
                                                        aaGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAACCGCTTg                                                                      <  2:462410/141‑1 (MQ=255)
                                                         aGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCATATGAAAACAAATAATAACAAAGCGCTTgg                                                                     >  1:56223/1‑140 (MQ=255)
                                                           aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCGGATAAAAACAAATAATAACAAAGCGCTTGCTa                                                                   <  2:205754/141‑1 (MQ=255)
                                                           aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTa                                                                   <  2:140085/141‑1 (MQ=255)
                                                            tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGGTGTAATTGATGCAAAAAAATTATCGAAAAAACTGCAGAACAATCTTCAGATAAAAACAAATAATAACAAAGAGCATGCTAg                                                                  >  2:434183/1‑141 (MQ=255)
                                                            tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATAAGTTGTAATTGATGCAAAAAAAGTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAACGCATGCTaa                                                                  >  2:486811/1‑140 (MQ=255)
                                                                aTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGAAGCTCATCAGTTGTAATGGATGCAAAAAAAGGATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTGGTAcc                                                              >  1:166105/1‑141 (MQ=255)
                                                                   gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                           >  2:526276/1‑141 (MQ=255)
                                                                   gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                           <  1:597198/141‑1 (MQ=255)
                                                                    agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCtt                                                          <  1:234514/141‑1 (MQ=255)
                                                                     gCAATCATTAAAAGAACCTGTGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAATGCGCCTGCTAGTACCTCGTa                                                         >  1:312989/1‑141 (MQ=255)
                                                                     gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                                                         >  2:58024/1‑141 (MQ=255)
                                                                        aTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaa                                                      <  1:144533/141‑1 (MQ=255)
                                                                        aTCAGTAAAAGGACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATACAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaa                                                      <  2:602030/141‑1 (MQ=255)
                                                                         tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAATAACAGGCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                     >  1:821815/1‑141 (MQ=255)
                                                                           aTTAAAAGAACCTGCGACTTATAAAGACGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACCAAGCGCTTGCTGGTACCACATAAAAAga                                                   >  2:450013/1‑141 (MQ=255)
                                                                            ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACCAATAATAACAAAGCGATTGCTAGTACATCTTAAAAAgag                                                  >  2:499758/1‑140 (MQ=255)
                                                                            ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTCCCTCTTAAAAAgat                                                  >  2:271792/1‑141 (MQ=255)
                                                                            ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat                                                  >  1:643733/1‑141 (MQ=255)
                                                                            ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGGTGTAATCGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCCTCAGATAAAAACAAATAATACCAAAGCGCTCGCTAGTACCTCTTaacaaaaa                                                  >  2:187055/1‑138 (MQ=255)
                                                                                  gAACCGGCAGCTTATAAATAAGACAGACGTTCATCAGGTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCt                                            <  2:696823/141‑1 (MQ=255)
                                                                                    aCCTGCGACTTATAAAGAAGACATACGTTCAGCAGTTGTAATTGATGCACAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATTCTAg                                          <  2:508202/141‑1 (MQ=255)
                                                                                      cTGCGACTTGTAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                        <  2:135175/141‑1 (MQ=255)
                                                                                      cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                        <  1:601247/141‑1 (MQ=255)
                                                                                      cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                        <  1:833793/141‑1 (MQ=255)
                                                                                       tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCaa                                       <  1:70219/141‑1 (MQ=255)
                                                                                       tGCGACTCATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGTAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAACaa                                       >  2:589469/1‑141 (MQ=255)
                                                                                       tGCGAATTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAACGCGCTTGCTCAAGTCCGTTATTAAGATGATGCTAGCaa                                       <  2:385452/141‑1 (MQ=255)
                                                                                        gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTATAAAGATGATGCTAGCAAg                                      >  1:777329/1‑141 (MQ=255)
                                                                                         cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgc                                     <  1:681894/141‑1 (MQ=255)
                                                                                         cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgc                                     >  2:799925/1‑141 (MQ=255)
                                                                                           aCTTATAGAGAACACAGACGTTCATCAGTTCTCATTGAGTTAATTAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTATCTTTTAAATAGATGATGCTAGCACgcgc                                   <  2:174058/141‑1 (MQ=255)
                                                                                           aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                   >  1:217117/1‑141 (MQ=255)
                                                                                               ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCCTGCTAGTACCTCTTAAAAGGATGATGCTAGCAAGCGAtttt                               >  2:662821/1‑141 (MQ=255)
                                                                                                taAAGAAGACAGACGTCCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                              >  1:498766/1‑141 (MQ=255)
                                                                                                 aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                             <  2:802917/141‑1 (MQ=255)
                                                                                                     agtacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatac                         <  2:474622/138‑1 (MQ=255)
                                                                                                     aagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAATACTGAAGGACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCTAGCGCTTTTctatac                         <  1:618272/141‑1 (MQ=255)
                                                                                                        acagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactat                      >  2:157095/1‑141 (MQ=255)
                                                                                                          agacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATCCtatat                    >  1:181988/1‑141 (MQ=255)
                                                                                                          agacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatat                    >  1:268343/1‑141 (MQ=255)
                                                                                                            acGTTCATCAGTTGTGATTGATGCACAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTatatat                  >  1:421706/1‑141 (MQ=255)
                                                                                                            acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAGACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCATGCTAGTACCTCTTCAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatatct                  >  2:232502/1‑139 (MQ=255)
                                                                                                             cGTTCATCAGTTGTAATTGCTGCATAACAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGAAAGCGCTTTTCTATACtatatagt                 >  2:180447/1‑139 (MQ=255)
                                                                                                             cGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATAtatt                 >  2:250005/1‑141 (MQ=255)
                                                                                                                 catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAttt             >  1:324568/1‑141 (MQ=255)
                                                                                                                   tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTCttttt           >  1:139023/1‑141 (MQ=255)
                                                                                                                     aGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTCGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCt         >  1:526299/1‑141 (MQ=255)
                                                                                                                           aaTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAtt   <  1:387033/141‑1 (MQ=255)
                                                                                                                            aTTGATGCAAAAAAATTATCGAACAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAttt  <  1:635854/141‑1 (MQ=255)
                                                                                                                            aTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGAGGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAttt  >  1:123535/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                     
AGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTT  >  USA300TCH1516_ALE/2455559‑2455824

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