Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,455,698 2 bp→AT 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,6980CA85.6% 245.0 / 39.6 92coding (1216/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (35/42);  minor base . (9/4);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01
*USA300TCH1516_ALE2,455,6990AT94.5% 263.5 / ‑6.7 93H406L (CAT→CTT) lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base T (43/43);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTG  >  USA300TCH1516_ALE/2455558‑2455841
                                                                                                                                             ||                                                                                                                                              
cagcagAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCACCAGTGGTAATTGATGCAAAAAAATTATCGaaa                                                                                                                                        >  2:771416/1‑149 (MQ=255)
   cagAATCAGTTAAATAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAAGAGAGCAATCATTAAAAGACCCTGCGACTTATAAAGAAGACAGACGTTCATTAGTTGTAATTGATGCAAAAAAATTATCTaaaaaa                                                                                                                                     <  2:111541/149‑1 (MQ=255)
   cagAATCAGTTAAACAAGCACCACAAGGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGGTGTAATTGATGCAAAAAAATTATCGaaaaaa                                                                                                                                     >  2:516410/1‑149 (MQ=255)
   cagAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTAAAATGAAGATGAAATAGAGCAATGATTAAAAGAACCTGAGACTTATAAAGAAGACAGACGTTCTGGAGTGGGAATTGATGCAAAAAAATTATCGaaaaaa                                                                                                                                     >  2:539652/1‑149 (MQ=255)
    agAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAAc                                                                                                                                    <  1:345446/149‑1 (MQ=255)
            ttAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGAAGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGACGAAAGACCTTCATCAGTCGAAAGTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                                                                            >  2:193254/1‑149 (MQ=255)
             tAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGTCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACa                                                                                                                           >  2:146631/1‑149 (MQ=255)
               aaCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAAt                                                                                                                         >  1:330864/1‑149 (MQ=255)
                aCAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        <  1:479495/149‑1 (MQ=255)
                 cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                                                                       >  2:597630/1‑149 (MQ=255)
                  aaGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCtt                                                                                                                      <  2:581124/149‑1 (MQ=255)
                  aaGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAAGAGAGCAATCATTAAAAGAACCTGCGACGTATAAAGAAGACAGACGTTCATAAGTTGTAATTGATGCAAAAAAAGTATCGAAAAAACTGTAGAACAAGCtt                                                                                                                      >  2:674080/1‑149 (MQ=255)
                   aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                                                                                     >  2:370117/1‑149 (MQ=255)
                     caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                                   <  2:236240/149‑1 (MQ=255)
                     caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                                   <  2:668034/149‑1 (MQ=255)
                     caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                                   <  2:286424/149‑1 (MQ=255)
                     caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                                   <  1:550577/149‑1 (MQ=255)
                        caaaaaGAGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAAAAAAAGAAAAAACTGAAGAACAATCTTCAGATa                                                                                                                <  1:281179/146‑1 (MQ=255)
                         acaAAGTGACGATTTAACAAACGATTCAGATGAAGAAGAAGTAGAGCAATCATTAAAAGAACCTGCGACTTATACAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGCTaa                                                                                                               >  2:385821/1‑149 (MQ=255)
                         acaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaa                                                                                                               >  2:293707/1‑149 (MQ=255)
                          caAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaa                                                                                                              >  1:256943/1‑149 (MQ=255)
                             aGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                                                                           <  2:163150/149‑1 (MQ=255)
                             aGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                                                                           <  1:367962/149‑1 (MQ=255)
                             aGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                                                                           <  1:374124/149‑1 (MQ=255)
                                gACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTAGAAAGAAGACAGACGTGCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACGGAAGAACAATCTTCAGATAAAAACaaa                                                                                                        <  2:114875/149‑1 (MQ=255)
                                gACGATTAAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaaa                                                                                                        <  2:753176/149‑1 (MQ=255)
                                 aCGATTTAACAAACGATTCAGATGAAGATAAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                                       >  1:481691/1‑149 (MQ=255)
                                 aCGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCGGATAAAAACAaat                                                                                                       >  2:288968/1‑149 (MQ=255)
                                 aCGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAACACTGCAGAACAATCTTCAGATAAAACCaccg                                                                                                       >  1:909717/1‑146 (MQ=255)
                                 aCGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                                       <  2:330864/149‑1 (MQ=255)
                                  cGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTACGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaata                                                                                                      >  2:381595/1‑149 (MQ=255)
                                     tttAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                                                                                   >  1:745216/1‑149 (MQ=255)
                                         aCAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                                                                               >  2:349355/1‑149 (MQ=255)
                                            aaagATTCAAATGAAGATGAAACATATCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                                                                                            >  2:905448/4‑149 (MQ=255)
                                             aCGATTCAAATGAAGATGAAATAGAGCAATCACTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcg                                                                                           >  2:453191/1‑149 (MQ=255)
                                               gATTCAAAGCAAGATGAACTCGAGCAATCATTAAAAGACCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCt                                                                                         <  2:917392/149‑1 (MQ=255)
                                                 ttCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                                                                                       <  1:642028/149‑1 (MQ=255)
                                                 ttCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                                                                                       <  1:828163/149‑1 (MQ=255)
                                                      aTGAAGATGAAATAGAGCAATCATTAAAAGAGCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAACTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  <  2:392156/149‑1 (MQ=255)
                                                      aTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  >  1:871934/1‑149 (MQ=255)
                                                      aTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  <  1:515107/149‑1 (MQ=255)
                                                        gAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                                                                                >  1:112771/1‑149 (MQ=255)
                                                          aGATGAAATCGAGCAACCATGAAAAGAACCTGCGACTTACAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAAAAAAAGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                                                                              <  2:597067/149‑1 (MQ=255)
                                                          aGATGAAATATCACAACAATTAAAAGAACGCCACCATTAAAAAGAACACAGACGTTCATCAGTTGAAATTGATGCAAAAAAATTACCGAAAACACTGAACAACAATTTGCAGATAAAAACAAATAATAACAAATCGCCTGCTAGTACct                                                                              <  2:154291/149‑1 (MQ=255)
                                                          aGATGAAATAGAGAAATCAATAGAAGAACCTGCGACTGATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAAATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                                                                              <  2:541715/149‑1 (MQ=255)
                                                           gATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGCAAGTACctc                                                                             >  2:277156/1‑149 (MQ=255)
                                                            gtGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAGAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                                            <  1:44041/148‑1 (MQ=255)
                                                            aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                                            <  1:339056/149‑1 (MQ=255)
                                                                 aTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaa                                                                       <  1:666322/149‑1 (MQ=255)
                                                                  tAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                                      <  1:535003/149‑1 (MQ=255)
                                                                  tAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAATACCGCTTaaaaa                                                                      >  1:791591/1‑149 (MQ=255)
                                                                  tAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                                      >  1:847589/1‑149 (MQ=255)
                                                                  aagagGAATAATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCGGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                                      <  2:677893/148‑1 (MQ=255)
                                                                    gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACATCTTAAAAAga                                                                    >  2:728371/1‑149 (MQ=255)
                                                                    gagCAATCATTAAAAGAACCTGCGACTTATAAAAAAGACAGACATTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACATCTTAAAAAga                                                                    >  2:132330/1‑149 (MQ=255)
                                                                     agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat                                                                   >  2:652803/1‑149 (MQ=255)
                                                                      gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCCTGCTAGTACCTCTTAACAAgatt                                                                  >  2:139534/1‑148 (MQ=255)
                                                                      gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAAAAAAGCACGAGCTAGTACCTCTTAAAAAgaaa                                                                  >  2:830902/1‑147 (MQ=255)
                                                                      gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAATTGTAATTGATGCAAAAAAATTATCGAAAAACCTGAAGAACAATCTTCAGATAAAAGCAAATAATAACAAAGCGCGTGCTATTACCACTTacaagcatg                                                                  >  1:261735/1‑144 (MQ=255)
                                                                          tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACACAGAATAACAAAGCCCTTGCTAGTACCTCTTAAAAAGATGATGc                                                              >  1:171401/1‑149 (MQ=255)
                                                                               aaaaGAACCGGTGACTTATAAACAAGACAGACGTTCATCAGTTTTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                                         <  2:463448/149‑1 (MQ=255)
                                                                                aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCaa                                                        >  1:924430/1‑149 (MQ=255)
                                                                                   gAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcg                                                     <  2:783726/149‑1 (MQ=255)
                                                                                   gAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcg                                                     <  1:96742/149‑1 (MQ=255)
                                                                                   gAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAAGTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAACAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcg                                                     <  2:59753/149‑1 (MQ=255)
                                                                                    aaaacTGCTACTCATAAAGAACACAGACGTTCATCAGTTGTAATTGATGCAAAAACATTATCTTAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcg                                                     <  1:215121/145‑1 (MQ=255)
                                                                                     aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                                   >  1:748012/1‑149 (MQ=255)
                                                                                         cgCGACATATAAAGGAGACAGACGTTCATCAGTGGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtttt                                                <  2:784139/148‑1 (MQ=255)
                                                                                          gtgACTTATAGCGACGCCAGATGTTCATGAGTTGTATTTTCTGCAAAAAAAATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  2:527903/147‑1 (MQ=255)
                                                                                          gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  2:283359/149‑1 (MQ=255)
                                                                                          gCGACTTATAAAGAAGACAGACGTGCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAACAGATGATGCTAGCAAGCGCTTTTc                                               <  1:552056/149‑1 (MQ=255)
                                                                                            gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATAATCAAAAAAACAGAATAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTcta                                             >  2:196616/1‑149 (MQ=255)
                                                                                            gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTcta                                             >  1:655741/1‑149 (MQ=255)
                                                                                            gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCATGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctt                                             >  2:722463/1‑148 (MQ=255)
                                                                                             aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctat                                            <  1:14494/149‑1 (MQ=255)
                                                                                             aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctat                                            <  1:294732/149‑1 (MQ=255)
                                                                                              cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACACAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctata                                           <  1:772209/149‑1 (MQ=255)
                                                                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatat                                     <  2:918832/149‑1 (MQ=255)
                                                                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatat                                     <  2:810744/149‑1 (MQ=255)
                                                                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatat                                     >  2:750798/1‑149 (MQ=255)
                                                                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatat                                     <  2:67436/149‑1 (MQ=255)
                                                                                                    aagaagACAGACGTTCATCAGTTGTAATTGATCCCAAAAAGGTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatat                                     >  2:199204/1‑149 (MQ=255)
                                                                                                              acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTc                           >  2:431717/1‑149 (MQ=255)
                                                                                                              acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTc                           >  2:918841/1‑149 (MQ=255)
                                                                                                               cGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCGTTTCTATACTATATATTATTTTTCt                          >  2:450888/1‑149 (MQ=255)
                                                                                                               cGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCt                          >  2:125660/1‑149 (MQ=255)
                                                                                                                     tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCAAATAACAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAtt                    <  1:727059/149‑1 (MQ=255)
                                                                                                                       aGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTGTAAAAtttt                  >  2:272679/1‑149 (MQ=255)
                                                                                                                           gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTGGCTAGTACCTCTTAAAAAGATGATGCTAGCACGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAaca              >  1:636903/1‑149 (MQ=255)
                                                                                                                            tAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAacac             <  2:460793/149‑1 (MQ=255)
                                                                                                                            tAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAacac             >  1:611460/1‑149 (MQ=255)
                                                                                                                                 gATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCTAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCctct        <  2:397388/149‑1 (MQ=255)
                                                                                                                                        aaaaaaTTATCGAAAAAACGGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCTAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTg  <  2:72267/149‑1 (MQ=255)
                                                                                                                                             ||                                                                                                                                              
CAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTG  >  USA300TCH1516_ALE/2455558‑2455841

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: