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,6980CA94.5% 563.2 / 22.2 167coding (1216/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base C (1/2);  major base A (64/91);  minor base . (5/4);  total (70/97)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
*USA300TCH1516_ALE2,455,6990AT98.2% 583.7 / 0.1 168H406L (CAT→CTT) lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (69/96);  total (70/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 169 total aligned reads displayed.

GAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGC  >  USA300TCH1516_ALE/2455614‑2455785
                                                                                    ||                                                                                      
aaaGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGa                                                                                  <  2:1685693/90‑1 (MQ=255)
 aaGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaa                                                                                 <  2:194366/91‑1 (MQ=255)
  aGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTGATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaa                                                                                <  2:830912/91‑1 (MQ=255)
    tcgAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                              >  2:769074/3‑91 (MQ=255)
    aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                              >  1:2513979/1‑91 (MQ=255)
    aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                              >  2:854713/1‑91 (MQ=255)
    aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                              <  1:1090179/91‑1 (MQ=255)
    aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                              <  1:95541/91‑1 (MQ=255)
     tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaaa                                                                             >  2:1518612/1‑91 (MQ=255)
      gAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAAc                                                                            >  2:1648810/1‑91 (MQ=255)
         aTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTga                                                                         >  1:804536/1‑91 (MQ=255)
          tAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaa                                                                        >  1:192429/1‑91 (MQ=255)
          tAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaa                                                                        >  1:1805374/1‑91 (MQ=255)
            gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                                      <  1:1525224/91‑1 (MQ=255)
            gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                                      <  1:2538806/91‑1 (MQ=255)
             agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                     <  1:2625023/91‑1 (MQ=255)
              gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                    >  1:2712398/1‑91 (MQ=255)
              gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                    <  2:1753841/91‑1 (MQ=255)
               cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACa                                                                   >  1:3083662/1‑91 (MQ=255)
                aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACaa                                                                  <  2:2177777/91‑1 (MQ=255)
                  tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                                >  2:1182170/1‑91 (MQ=255)
                   cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                               <  1:2317800/91‑1 (MQ=255)
                     ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                             <  1:1013066/91‑1 (MQ=255)
                      tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                            >  2:2362274/1‑91 (MQ=255)
                        aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                          >  1:435314/1‑91 (MQ=255)
                        aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                          <  1:2231384/91‑1 (MQ=255)
                         aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                         <  2:2410892/91‑1 (MQ=255)
                         aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                         <  2:2412676/91‑1 (MQ=255)
                            aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaa                                                      <  2:2139238/91‑1 (MQ=255)
                            aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaa                                                      <  2:1863388/91‑1 (MQ=255)
                             aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                     <  2:78602/91‑1 (MQ=255)
                             aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                     >  1:763098/1‑91 (MQ=255)
                             aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                     >  1:2252126/1‑91 (MQ=255)
                              ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaaa                                                    >  1:2699596/1‑91 (MQ=255)
                                tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                                  <  1:99439/91‑1 (MQ=255)
                                tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                                  <  1:2293525/91‑1 (MQ=255)
                                tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                                  <  2:237588/91‑1 (MQ=255)
                                 gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaa                                                 >  2:1023328/1‑91 (MQ=255)
                                 gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaa                                                 >  1:1911316/1‑91 (MQ=255)
                                 gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaa                                                 >  2:1659507/1‑91 (MQ=255)
                                  cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaaa                                                >  2:2198286/1‑91 (MQ=255)
                                   gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                               <  1:811262/91‑1 (MQ=255)
                                     cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataa                                             >  1:157444/1‑91 (MQ=255)
                                     cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataa                                             >  2:1476169/1‑91 (MQ=255)
                                      ttATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                            <  2:2802653/91‑1 (MQ=255)
                                       tataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                           <  2:1596086/91‑1 (MQ=255)
                                          aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACa                                        >  2:3058509/1‑91 (MQ=255)
                                          aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACa                                        >  2:23498/1‑91 (MQ=255)
                                          aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACa                                        >  2:1862967/1‑91 (MQ=255)
                                          aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACa                                        <  1:2845347/91‑1 (MQ=255)
                                           aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                       <  1:1452891/91‑1 (MQ=255)
                                           aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                       >  1:1480052/1‑91 (MQ=255)
                                           aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                       <  1:1131888/91‑1 (MQ=255)
                                           aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                       <  1:1778508/91‑1 (MQ=255)
                                           aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                       <  2:1482918/91‑1 (MQ=255)
                                           aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                       <  2:213345/91‑1 (MQ=255)
                                             gaagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAg                                     <  2:456456/91‑1 (MQ=255)
                                              aagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                                    <  2:2532809/91‑1 (MQ=255)
                                              aagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                                    >  2:2811259/1‑91 (MQ=255)
                                               agacagacGTTCATCGGTTGTAATTAGTGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcg                                   <  2:2688888/91‑1 (MQ=255)
                                                gacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgc                                  <  2:1790664/91‑1 (MQ=255)
                                                gacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgc                                  <  1:1594754/91‑1 (MQ=255)
                                                 acagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAACAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCt                                 >  1:111225/1‑91 (MQ=255)
                                                   agacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgcgag                               >  1:66018/1‑88 (MQ=255)
                                                   agacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                               <  1:133456/91‑1 (MQ=255)
                                                    gacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                              >  1:1347717/1‑91 (MQ=255)
                                                    gacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                              <  1:1409995/91‑1 (MQ=255)
                                                    gacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                              <  2:3099159/91‑1 (MQ=255)
                                                    gacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                              >  2:1234104/1‑91 (MQ=255)
                                                      cGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTa                            <  2:306932/91‑1 (MQ=255)
                                                       gTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAg                           >  2:36339/1‑91 (MQ=255)
                                                        ttCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                          <  1:3058762/91‑1 (MQ=255)
                                                          catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                        <  2:2652571/91‑1 (MQ=255)
                                                          catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                        >  1:454578/1‑91 (MQ=255)
                                                          catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                        <  1:1280203/91‑1 (MQ=255)
                                                          catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                        <  1:497858/91‑1 (MQ=255)
                                                           atcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAcc                       <  1:1697656/91‑1 (MQ=255)
                                                            tcaGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAATTAATAACAAAGCGCTTGCTAGTAcc                       <  1:1972715/91‑1 (MQ=255)
                                                            tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                      >  1:1944787/1‑91 (MQ=255)
                                                            tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                      <  1:1034080/91‑1 (MQ=255)
                                                            tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                      <  1:1639628/91‑1 (MQ=255)
                                                             caGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctc                     <  1:3048019/91‑1 (MQ=255)
                                                              aGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                    <  1:1462006/91‑1 (MQ=255)
                                                              aGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                    <  2:1189348/91‑1 (MQ=255)
                                                              aGTTGTAATGGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                    <  1:2977064/91‑1 (MQ=255)
                                                                ttGTAATTGATGCAAATAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                  <  2:2971555/91‑1 (MQ=255)
                                                                ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                  >  1:87873/1‑91 (MQ=255)
                                                                ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                  <  2:2972840/91‑1 (MQ=255)
                                                                ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                  <  2:663508/91‑1 (MQ=255)
                                                                 tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                 <  1:384881/91‑1 (MQ=255)
                                                                 tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                 <  1:384869/91‑1 (MQ=255)
                                                                  gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa                <  1:1863512/91‑1 (MQ=255)
                                                                   tAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaa               <  1:2016198/91‑1 (MQ=255)
                                                                     aTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg             >  2:538480/1‑91 (MQ=255)
                                                                     aTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAGAGCGCTTGCTAGTACCTCTTAAAAAg             <  2:1758657/91‑1 (MQ=255)
                                                                     aTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg             >  2:726482/1‑91 (MQ=255)
                                                                      ttGATGCAAAAAAATTATCGAAAAAACTGAAGAACAACCTTCAGATAAAAATAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAga            >  1:2885057/1‑91 (MQ=255)
                                                                       tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat           <  2:2788760/91‑1 (MQ=255)
                                                                       tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat           >  1:319193/1‑91 (MQ=255)
                                                                                 agaaatatCGATAAAACTGAAGAACAATCTTCAGATAAAAACAGATAATCGCATAGCGGTTGCTAGTAACACATAAAAAGATGATGCTAGc  <  2:3045798/87‑1 (MQ=255)
                                                                                    ||                                                                                      
GAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGC  >  USA300TCH1516_ALE/2455614‑2455785

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