Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 41,171 T→C 6.7% intergenic (‑38/‑62) pbp ← / → mecR1 Beta‑lactam‑inducible penicillin‑binding protein/Methicillin resistance mecR1 protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE41,1710TC 6.7% 478.1 / 16.6 148intergenic (‑38/‑62)pbp/mecR1Beta‑lactam‑inducible penicillin‑binding protein/Methicillin resistance mecR1 protein
Reads supporting (aligned to +/- strand):  ref base T (69/69);  new base C (7/3);  total (76/72)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 156 total aligned reads displayed.

CCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGTA  >  USA300TCH1516_ALE/41083‑41260
                                                                                        |                                                                                         
cccGACAACTACAACTATTAAAATAAGTGTATCAATTTTTATCTCTTTCATCATTATTCTCCTTATATAAGACTACATTTGTAATATACTa                                                                                         <  1:3317151/91‑1 (MQ=255)
cccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTa                                                                                         <  1:979200/91‑1 (MQ=255)
 ccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATAtactac                                                                                        <  1:2346249/91‑1 (MQ=255)
  cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa                                                                                       >  1:2878784/1‑91 (MQ=255)
  cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa                                                                                       >  1:1356896/1‑91 (MQ=255)
  cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa                                                                                       >  1:1356811/1‑91 (MQ=255)
   gACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATATACTACaa                                                                                      <  1:3906517/91‑1 (MQ=255)
    acaactacaactATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACaaa                                                                                     >  1:3891640/1‑91 (MQ=255)
     caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt                                                                                    >  1:2433281/1‑91 (MQ=255)
     caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt                                                                                    >  1:634998/1‑91 (MQ=255)
      aactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATg                                                                                   <  1:3719715/91‑1 (MQ=255)
      aactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATATACTACAAATg                                                                                   <  1:1624225/91‑1 (MQ=255)
       actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgt                                                                                  <  1:3125769/91‑1 (MQ=255)
       actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgt                                                                                  <  1:3129705/91‑1 (MQ=255)
       actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGCATATTACAAATgt                                                                                  <  1:3124710/91‑1 (MQ=255)
              ctaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTAttt                                                                           <  1:810549/91‑1 (MQ=255)
                 ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg                                                                        >  1:12880/1‑91 (MQ=255)
                 ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg                                                                        >  1:2628741/1‑91 (MQ=255)
                    aaaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCa                                                                     >  1:93012/1‑91 (MQ=255)
                     aaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCaa                                                                    <  1:2671176/91‑1 (MQ=255)
                     aaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCaa                                                                    >  1:1053564/1‑91 (MQ=255)
                        aaGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaat                                                                 >  1:649159/1‑91 (MQ=255)
                          gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa                                                               >  1:1842909/1‑91 (MQ=255)
                          gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa                                                               >  1:2236317/1‑91 (MQ=255)
                          gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa                                                               >  1:1842807/1‑91 (MQ=255)
                            ggAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg                                                             >  1:1211714/1‑91 (MQ=255)
                            ggAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg                                                             >  1:2285655/1‑91 (MQ=255)
                            ggAACAATNTTTNTCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg                                                             <  1:422027/91‑1 (MQ=255)
                              aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCATAATAATGtt                                                           <  1:2091198/91‑1 (MQ=255)
                              aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                           <  1:1088068/91‑1 (MQ=255)
                              aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                           <  1:1157088/91‑1 (MQ=255)
                              aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                           <  1:869004/91‑1 (MQ=255)
                              aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                           <  1:1541637/91‑1 (MQ=255)
                              aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                           <  1:3543157/91‑1 (MQ=255)
                                cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat                                                         <  1:1514244/91‑1 (MQ=255)
                                cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat                                                         <  1:304695/91‑1 (MQ=255)
                                cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat                                                         <  1:1784841/91‑1 (MQ=255)
                                 aaTTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtata                                                        <  1:3139337/91‑1 (MQ=255)
                                 aaTTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtata                                                        >  1:3800018/1‑91 (MQ=255)
                                     tttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttt                                                    >  1:1026460/1‑91 (MQ=255)
                                      ttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAtttt                                                   >  1:1878572/1‑91 (MQ=255)
                                       tATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttttt                                                  <  1:1094713/91‑1 (MQ=255)
                                       tATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttttt                                                  >  1:2860091/1‑91 (MQ=255)
                                         tctcTTTCATCAATATACTCCTTAAATAAGACTACATTTGTAGTATATTACAAATGTAGGATCCATGTCAAAATAATGTTATAATTTTtgt                                                <  1:2976150/91‑1 (MQ=255)
                                         tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                                <  1:2806776/91‑1 (MQ=255)
                                         tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                                >  1:2750313/1‑91 (MQ=255)
                                         tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                                <  1:986049/91‑1 (MQ=255)
                                         tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                                >  1:2709413/1‑91 (MQ=255)
                                             tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGata                                            >  1:2646537/1‑91 (MQ=255)
                                             tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGata                                            <  1:2021960/91‑1 (MQ=255)
                                              ttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat                                           <  1:117452/91‑1 (MQ=255)
                                              ttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat                                           <  1:770793/91‑1 (MQ=255)
                                              ttCATCAATATACTCCTCATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat                                           <  1:3556564/91‑1 (MQ=255)
                                                  tcaATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAg                                       <  1:2694215/91‑1 (MQ=255)
                                                  tcaATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAg                                       <  1:595764/91‑1 (MQ=255)
                                                    aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                     >  1:1414036/1‑91 (MQ=255)
                                                    aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                     >  1:85107/1‑91 (MQ=255)
                                                    aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                     <  1:2634061/91‑1 (MQ=255)
                                                     atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                    >  1:3781354/1‑91 (MQ=255)
                                                     atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                    >  1:2034120/1‑91 (MQ=255)
                                                     atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                    >  1:3303407/1‑91 (MQ=255)
                                                     atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                    <  1:3238268/91‑1 (MQ=255)
                                                           tCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAg                              >  1:2311903/1‑91 (MQ=255)
                                                           tCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAg                              <  1:1498097/91‑1 (MQ=255)
                                                               tatataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGgtgt                          >  1:907967/1‑91 (MQ=255)
                                                               tatataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGgtgt                          >  1:2158678/1‑91 (MQ=255)
                                                                atataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGtt                         >  1:172600/1‑91 (MQ=255)
                                                                atataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGtt                         <  1:1591117/91‑1 (MQ=255)
                                                                 tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                        <  1:562934/91‑1 (MQ=255)
                                                                 tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                        >  1:2406334/1‑91 (MQ=255)
                                                                 tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                        >  1:2461675/1‑91 (MQ=255)
                                                                 tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                        <  1:1449401/91‑1 (MQ=255)
                                                                 tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                        >  1:2405431/1‑91 (MQ=255)
                                                                  ataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTat                       <  1:2633670/91‑1 (MQ=255)
                                                                   taAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatc                      >  1:2088175/1‑91 (MQ=255)
                                                                   taAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatc                      <  1:4082468/91‑1 (MQ=255)
                                                                      gACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATGATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc                   <  1:3506623/91‑1 (MQ=255)
                                                                      gACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc                   <  1:2054786/91‑1 (MQ=255)
                                                                      gACTACATTTGTAGTATATTACAAATGTAGTATGTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc                   <  1:3681032/91‑1 (MQ=255)
                                                                         tACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttt                <  1:3291291/91‑1 (MQ=255)
                                                                         tACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttt                <  1:3291175/91‑1 (MQ=255)
                                                                          aCATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt               >  1:974680/1‑91 (MQ=255)
                                                                          aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt               >  1:2776684/1‑91 (MQ=255)
                                                                          aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt               >  1:406530/1‑91 (MQ=255)
                                                                            aTTTGTAGTATATTACAGATGTAGTATTTATGGCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttttt             <  1:3546833/91‑1 (MQ=255)
                                                                            aTTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttttt             >  1:3567834/1‑91 (MQ=255)
                                                                            aTTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttttt             <  1:554697/91‑1 (MQ=255)
                                                                             tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTa            >  1:2563806/1‑91 (MQ=255)
                                                                              ttGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTaa           >  1:4101680/1‑91 (MQ=255)
                                                                               tGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAAt          >  1:1233240/1‑91 (MQ=255)
                                                                                gtagtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATg         <  1:4008565/91‑1 (MQ=255)
                                                                                 tagtaTATTACGAATGTGGGATTTATGTCGAAATAATGTTAGAATTTTTGTGGTATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt        <  1:3333042/91‑1 (MQ=255)
                                                                                 taatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt        >  1:2034848/4‑91 (MQ=255)
                                                                                  agtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt       <  1:345111/91‑1 (MQ=255)
                                                                                  aatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt       >  1:2430913/3‑91 (MQ=255)
                                                                                   gtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTa      >  1:2174465/1‑91 (MQ=255)
                                                                                     atatTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAg    <  1:3957117/91‑1 (MQ=255)
                                                                                     atatTACAAATGGAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAg    >  1:3200796/1‑91 (MQ=255)
                                                                                      tatTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGt   <  1:237308/91‑1 (MQ=255)
                                                                                       atTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGta  <  1:3114916/91‑1 (MQ=255)
                                                                                        |                                                                                         
CCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGTA  >  USA300TCH1516_ALE/41083‑41260

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