Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,455,698 2 bp→AT 100% coding (1216‑1217/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,6980CA78.7% 189.9 / 53.2 75coding (1216/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (29/29);  minor base . (13/3);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*USA300TCH1516_ALE2,455,6990AT98.2% 217.1 / ‑3.4 77H406L (CAT→CTT) lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (40/34);  minor base G (2/0);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAG  >  USA300TCH1516_ALE/2455566‑2455839
                                                                                                                                    ||                                                                                                                                             
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)
                                                                                                                             caaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACAcc          <  1:714541/139‑1 (MQ=255)
                                                                                                                               aaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCct        <  1:46826/141‑1 (MQ=255)
                                                                                                                                     ataaCGAAAAAACTGAAGAACAATCTGCAAATAAAAACAAAGAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAg  <  1:342897/141‑1 (MQ=255)
                                                                                                                                    ||                                                                                                                                             
TCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAG  >  USA300TCH1516_ALE/2455566‑2455839

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: