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,6980CA96.2% 621.4 / 15.0 185coding (1216/1260 nt)lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base C (1/3);  major base A (83/89);  minor base . (4/3);  total (89/96)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*USA300TCH1516_ALE2,455,6990AT97.2% 626.6 / 3.5 189H406L (CAT→CTT) lyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base T (87/94);  total (90/99)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 189 total aligned reads displayed.

AAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAG  >  USA300TCH1516_ALE/2455615‑2455788
                                                                                   ||                                                                                         
  gATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaa                                                                                  <  1:1041646/91‑1 (MQ=255)
  gATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaa                                                                                  >  2:1290321/1‑91 (MQ=255)
    tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaaa                                                                                >  1:852969/1‑91 (MQ=255)
    tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGaaaaaa                                                                                >  1:3694215/1‑91 (MQ=255)
      aaaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAATTATCGAAAAAACTg                                                                             >  1:1822645/1‑91 (MQ=255)
       aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTg                                                                             >  1:1433285/1‑91 (MQ=255)
       aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTg                                                                             >  1:2754916/1‑91 (MQ=255)
       aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATAGATACAAAAAAATTATCGAAAAAACTg                                                                             <  2:2499731/91‑1 (MQ=255)
          agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                                          <  2:1904259/91‑1 (MQ=255)
          agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                                          >  1:2720567/1‑91 (MQ=255)
          agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                                          <  2:2006820/91‑1 (MQ=255)
          agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaag                                                                          <  1:3854504/91‑1 (MQ=255)
           gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                                         >  2:3618850/1‑91 (MQ=255)
           gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaaga                                                                         <  1:693152/91‑1 (MQ=255)
            agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                        >  1:2945014/1‑91 (MQ=255)
            agCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                        >  2:3524308/1‑91 (MQ=255)
            agCAATCATTAAAAGAACCTGCGACGTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTgaagaa                                                                        <  2:3804783/91‑1 (MQ=255)
             gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTAATGCAAAAAAATTATCGAAAAAACTGAAGAAc                                                                       >  1:3311230/1‑91 (MQ=255)
              cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACa                                                                      <  1:1204295/91‑1 (MQ=255)
              cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACa                                                                      >  1:2082039/1‑91 (MQ=255)
               aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACaa                                                                     >  1:3573615/1‑91 (MQ=255)
               aaTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACaa                                                                     >  2:1501954/1‑91 (MQ=255)
                 tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                                   <  1:3783300/91‑1 (MQ=255)
                 tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                                   >  1:3627582/1‑91 (MQ=255)
                 tCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATc                                                                   <  2:1544051/91‑1 (MQ=255)
                  cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                  <  1:2238316/91‑1 (MQ=255)
                  cATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCt                                                                  >  1:3718850/1‑91 (MQ=255)
                   aTTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCtt                                                                 >  2:339794/1‑91 (MQ=255)
                   aTTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCtt                                                                 >  2:166893/1‑91 (MQ=255)
                    ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                                <  2:2465467/91‑1 (MQ=255)
                    ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                                >  2:1405897/1‑91 (MQ=255)
                    ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                                <  1:2611704/91‑1 (MQ=255)
                    ttAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGCTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                                >  2:3002102/1‑91 (MQ=255)
                    ttAAAAGAACCTGCGACTTATAAAAAAAACAGACGTCCATCAGTTCTAATTGATGCAAAAAAATAATCGAAAAAACGGAAAAACAACCTTc                                                                >  2:104646/1‑91 (MQ=255)
                     tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                               >  2:2494663/1‑91 (MQ=255)
                     tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                               >  1:2804685/1‑91 (MQ=255)
                     tAAAAGAACCTGCGACTTATAAAGAAGACAGACGGTCATCAGTTGTAAGAGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                               <  2:474368/91‑1 (MQ=255)
                      aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                              >  1:523122/1‑91 (MQ=255)
                       aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                             <  2:3132962/91‑1 (MQ=255)
                       aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                             <  1:2110817/91‑1 (MQ=255)
                       aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                             >  1:3865390/1‑91 (MQ=255)
                       aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                             >  1:3745088/1‑91 (MQ=255)
                        aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                            <  1:2280494/91‑1 (MQ=255)
                        aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                            <  1:405293/91‑1 (MQ=255)
                        aaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                            <  1:2281634/91‑1 (MQ=255)
                         aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATa                                                           <  1:1625856/91‑1 (MQ=255)
                           aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaa                                                         <  1:3253078/91‑1 (MQ=255)
                             ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaaa                                                       >  1:540541/1‑91 (MQ=255)
                             ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaaa                                                       >  2:1809508/1‑91 (MQ=255)
                              cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                      <  2:2226040/91‑1 (MQ=255)
                               tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                                     <  2:2378399/91‑1 (MQ=255)
                               tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                                     <  1:3124441/91‑1 (MQ=255)
                               tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACa                                                     >  2:846860/1‑91 (MQ=255)
                                 cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaaa                                                   >  2:782731/1‑91 (MQ=255)
                                 cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaaa                                                   >  2:3345594/1‑91 (MQ=255)
                                  ggcTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAATAAAATTATTGAACAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                  <  1:3093168/89‑1 (MQ=255)
                                  gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAAATATAGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                  <  2:1230464/91‑1 (MQ=255)
                                   aCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaata                                                 >  1:2049676/1‑91 (MQ=255)
                                       ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAAtaataa                                             <  1:211703/91‑1 (MQ=255)
                                        taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAAc                                            >  1:3841949/1‑91 (MQ=255)
                                         aaaGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACa                                           >  1:337193/1‑91 (MQ=255)
                                          aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                          <  2:2808104/91‑1 (MQ=255)
                                          aagaagACAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaa                                          <  2:1364265/91‑1 (MQ=255)
                                           agaagaCAGACGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACaaa                                         >  1:334809/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)
                                                acagacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCt                                    >  1:2009985/1‑91 (MQ=255)
                                                 cagacGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCt                                    <  1:3143785/91‑1 (MQ=255)
                                                  agacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                                  >  1:1071973/1‑91 (MQ=255)
                                                  agacGTTCATCAGTTGTAAATAAAACAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                                  <  2:3429380/91‑1 (MQ=255)
                                                   gacGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                                 <  2:2910645/91‑1 (MQ=255)
                                                    acGTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCt                                >  1:3404425/1‑91 (MQ=255)
                                                      gTTCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAg                              <  1:3160506/91‑1 (MQ=255)
                                                       ttCATCAGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                             <  1:2147755/91‑1 (MQ=255)
                                                        tcatcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTa                            <  2:3303858/91‑1 (MQ=255)
                                                         catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTa                               >  1:1744405/1‑87 (MQ=255)
                                                         catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTa                               <  2:1744405/87‑1 (MQ=255)
                                                         catcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                           >  2:1779774/1‑91 (MQ=255)
                                                           tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                         <  2:881915/91‑1 (MQ=255)
                                                           tcaGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                         >  1:2715215/1‑91 (MQ=255)
                                                            caGTTGTAATTGATGCAAAAAAATTATCGAAAACACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTCGCTAGTACctc                        >  2:3790925/1‑91 (MQ=255)
                                                            caGTTGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctc                        <  1:2151039/91‑1 (MQ=255)
                                                               ttGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                     >  1:696014/1‑91 (MQ=255)
                                                                tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                    <  2:975923/91‑1 (MQ=255)
                                                                tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                    >  2:371074/1‑91 (MQ=255)
                                                                tGTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTCGCTAGTACCTCTTaa                    >  2:372643/1‑91 (MQ=255)
                                                                 gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa                   <  2:3069309/91‑1 (MQ=255)
                                                                 gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa                   <  2:3066919/91‑1 (MQ=255)
                                                                 gTAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaa                   <  1:812505/91‑1 (MQ=255)
                                                                  tAATTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaa                  <  2:696014/91‑1 (MQ=255)
                                                                    aTTGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg                <  2:2179060/91‑1 (MQ=255)
                                                                     ttGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAga               <  2:211758/91‑1 (MQ=255)
                                                                      tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat              <  1:440205/91‑1 (MQ=255)
                                                                      tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat              <  1:2553446/91‑1 (MQ=255)
                                                                      tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat              >  1:1364127/1‑91 (MQ=255)
                                                                      tGATGCAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat              <  2:2420357/91‑1 (MQ=255)
                                                                          tgcaaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatgat           <  2:660192/87‑1 (MQ=255)
                                                                           gcaaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGatgatg          <  2:191669/88‑1 (MQ=255)
                                                                                   cATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAg  >  2:2658623/1‑91 (MQ=255)
                                                                                   ||                                                                                         
AAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAG  >  USA300TCH1516_ALE/2455615‑2455788

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: