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.93.6% 757.9 / 33.4 185coding (1210/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base G (0/3);  major base . (84/89);  minor base A (5/4);  total (89/96)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
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.

Only 100 of 205 total aligned reads displayed.

GATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGAT  >  USA300TCH1516_ALE/2455605‑2455776
                                                                                       |                                                                                    
gATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaa                                                                                   <  1:210107/91‑1 (MQ=255)
 aTTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCaaaaaa                                                                                  >  1:2739483/1‑91 (MQ=255)
  ttCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaa                                                                                 >  1:563495/1‑91 (MQ=255)
  ttCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaa                                                                                 >  2:1669289/1‑91 (MQ=255)
   tCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaat                                                                                >  1:1795167/1‑90 (MQ=255)
    cAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaatt                                                                               >  1:1531092/1‑89 (MQ=255)
     aaaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaatta                                                                              >  2:2119459/1‑88 (MQ=255)
      aaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattat                                                                             <  1:2294311/91‑5 (MQ=255)
      aaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAaaaaaattat                                                                             >  2:3334108/1‑87 (MQ=255)
          aaGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaa                                                                        >  2:893261/1‑91 (MQ=255)
            gATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaa                                                                      <  1:1041646/91‑1 (MQ=255)
             aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                     <  1:3392750/91‑1 (MQ=255)
             aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaa                                                                     <  1:3655708/91‑1 (MQ=255)
              tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaaa                                                                    >  1:852969/1‑91 (MQ=255)
                 aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTg                                                                 >  1:1433285/1‑91 (MQ=255)
                 aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTg                                                                 <  1:3843882/91‑1 (MQ=255)
                    agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                              <  1:3584397/91‑1 (MQ=255)
                    agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                              >  1:2720567/1‑91 (MQ=255)
                    agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                              <  1:3854504/91‑1 (MQ=255)
                    agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                              <  2:2006676/91‑1 (MQ=255)
                     gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                             >  2:178431/1‑91 (MQ=255)
                     gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                             <  1:693152/91‑1 (MQ=255)
                     gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                             <  2:1843783/91‑1 (MQ=255)
                      agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                            <  2:2934412/91‑1 (MQ=255)
                      agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                            >  1:2945014/1‑91 (MQ=255)
                      agCAATCATTAAAAGAACCTGCGACGTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                            <  2:3804783/91‑1 (MQ=255)
                       gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                           <  2:1054529/91‑1 (MQ=255)
                        cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATAATCGAAAAAACTGAAGAACa                                                          <  2:330975/91‑1 (MQ=255)
                        cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACa                                                          >  1:2082039/1‑91 (MQ=255)
                         aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACaa                                                         >  1:2226040/1‑91 (MQ=255)
                         aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACaa                                                         >  2:1501954/1‑91 (MQ=255)
                         aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACaa                                                         <  1:2031173/91‑1 (MQ=255)
                          atcatcAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAAt                                                        <  2:3579311/91‑1 (MQ=255)
                           tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                       >  1:3292946/1‑91 (MQ=255)
                           tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                       <  2:1544051/91‑1 (MQ=255)
                           acATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATAGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                       <  2:103877/90‑1 (MQ=255)
                            cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                      >  1:3718850/1‑91 (MQ=255)
                             aTTAAAAGAACCTGCGACTTATACAGAAGACAGACGTTCATCAGAAGTAATTGCTGCAAAAAAACTATCGAAAAAACTGAAGAACAATCtt                                                     >  2:2605055/1‑91 (MQ=255)
                             aTTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCt                                                      <  2:3457364/91‑1 (MQ=255)
                             aTTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCtt                                                     >  1:1294002/1‑91 (MQ=255)
                              ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                    >  2:1405897/1‑91 (MQ=255)
                              ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                    <  2:2465467/91‑1 (MQ=255)
                              ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                    <  2:3004755/91‑1 (MQ=255)
                              ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                    <  1:2611704/91‑1 (MQ=255)
                               tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                   >  2:2494663/1‑91 (MQ=255)
                               tAAAAGAACCTGCGACTTATAAAGAAGACAGACGGTCATCAGTTGTAAGAGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                   <  2:474368/91‑1 (MQ=255)
                                aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                  >  1:523122/1‑91 (MQ=255)
                                 aaaGAACCTGCGAGTTATAAAGAAGACAGACGTTCATCAGTAGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                 <  2:3420855/91‑1 (MQ=255)
                                 aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                 >  1:3745088/1‑91 (MQ=255)
                                 aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                 <  1:2111340/91‑1 (MQ=255)
                                 aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                 >  1:3865390/1‑91 (MQ=255)
                                  aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                <  1:405293/91‑1 (MQ=255)
                                  aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                <  1:2280536/91‑1 (MQ=255)
                                   aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATAGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATa                                               <  2:816771/91‑1 (MQ=255)
                                      aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                            >  1:3443053/1‑91 (MQ=255)
                                      aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                            >  2:368155/1‑91 (MQ=255)
                                      aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAAAAGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                            <  2:705883/91‑1 (MQ=255)
                                       ccTGCGACTAATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCACATaaaaa                                           >  1:1946098/1‑91 (MQ=255)
                                        cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                          <  2:2226040/91‑1 (MQ=255)
                                         tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                         <  1:3124441/91‑1 (MQ=255)
                                         tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                         <  2:2378399/91‑1 (MQ=255)
                                         tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                         >  2:599580/1‑91 (MQ=255)
                                          gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaa                                        >  2:3111083/1‑91 (MQ=255)
                                          gCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaa                                        >  1:2782686/1‑91 (MQ=255)
                                            ggcTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAATAAAATTATTGAACAAACTGAAGAACAATCTTCAGATAAAAACAaat                                      <  1:3093168/89‑1 (MQ=255)
                                                 ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAAtaataa                                 <  1:1787372/91‑1 (MQ=255)
                                                  taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAAc                                >  2:1349292/1‑91 (MQ=255)
                                                   aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACa                               <  1:1384163/91‑1 (MQ=255)
                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                              <  1:2120930/91‑1 (MQ=255)
                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                              <  1:3322671/91‑1 (MQ=255)
                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                              <  1:3796202/91‑1 (MQ=255)
                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                              <  2:2808104/91‑1 (MQ=255)
                                                    aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                              <  1:2846453/91‑1 (MQ=255)
                                                      gaagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAg                            >  1:2249130/1‑91 (MQ=255)
                                                       aagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc                           >  2:3151481/1‑91 (MQ=255)
                                                         gacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgc                         <  2:2825682/91‑1 (MQ=255)
                                                         aacagacGTTCATCAGTTGTAATTGATGCAAAAAAATTAAAGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgc                         <  2:1610158/90‑1 (MQ=255)
                                                           cagacGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCt                        <  1:3143785/91‑1 (MQ=255)
                                                             gacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAGAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                     >  1:3365396/1‑91 (MQ=255)
                                                              acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCt                    <  1:1161532/91‑1 (MQ=255)
                                                               cGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTa                   >  1:3248909/1‑91 (MQ=255)
                                                                   catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTa                   <  2:1744405/87‑1 (MQ=255)
                                                                   catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc               >  2:2895820/1‑91 (MQ=255)
                                                                   catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc               >  2:1779774/1‑91 (MQ=255)
                                                                     tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct             >  1:2715215/1‑91 (MQ=255)
                                                                     tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct             >  2:3378609/1‑91 (MQ=255)
                                                                      caGTTGTAATTGATGCAAAAAAATTATCGAAAACACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTCGCTAGTACctc            >  2:3790925/1‑91 (MQ=255)
                                                                      caGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctc            <  1:2150368/91‑1 (MQ=255)
                                                                       aGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct           <  2:2087660/91‑1 (MQ=255)
                                                                         ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa         >  2:879775/1‑91 (MQ=255)
                                                                         ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa         >  1:696014/1‑91 (MQ=255)
                                                                          tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa        >  1:174213/1‑91 (MQ=255)
                                                                          tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa        >  2:371075/1‑91 (MQ=255)
                                                                           gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa       <  1:812505/91‑1 (MQ=255)
                                                                           gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa       <  2:3069309/91‑1 (MQ=255)
                                                                              aTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg    <  2:2179060/91‑1 (MQ=255)
                                                                                tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat  <  1:440205/91‑1 (MQ=255)
                                                                                tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat  >  2:1241076/1‑91 (MQ=255)
                                                                                tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat  >  2:3190691/1‑91 (MQ=255)
                                                                                tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat  <  2:2420357/91‑1 (MQ=255)
                                                                                       |                                                                                    
GATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGAT  >  USA300TCH1516_ALE/2455605‑2455776

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: