Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 41,171 T→C 14.1% 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,1710TC14.1% 891.5 / 119.0 349intergenic (‑38/‑62)pbp/mecR1Beta‑lactam‑inducible penicillin‑binding protein/Methicillin resistance mecR1 protein
Reads supporting (aligned to +/- strand):  ref base T (178/121);  new base C (35/14);  total (214/135)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-01

Only 100 of 381 total aligned reads displayed.

AACCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGTAT  >  USA300TCH1516_ALE/41081‑41261
                                                                                           |                                                                                          
aaCCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATGCTCCTTATATAAGACTACATTTGTAGTatat                                                                                            <  1:2647675/91‑1 (MQ=255)
 aCCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATAtt                                                                                           >  2:29372/1‑91 (MQ=255)
  cccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTa                                                                                          <  2:285117/91‑1 (MQ=255)
  cccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTa                                                                                          <  1:2245233/91‑1 (MQ=255)
   ccGACCGCTACAACTATTAAAATAAGTGGAACAATTTTTATCTTGTTCATCAATATACTCCTTATATAAGACTACATTTGTAATAtactac                                                                                         <  2:3642355/91‑1 (MQ=255)
   ccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTATCATCAATATACTCCTTATATAAGACTACATTTGTAATAtactac                                                                                         <  1:876080/91‑1 (MQ=255)
    cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATAtactac                                                                                         <  1:3043787/91‑1 (MQ=255)
    cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATATACTACa                                                                                        <  1:1329838/91‑1 (MQ=255)
      acaactacaactATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACaaa                                                                                      >  2:311399/1‑91 (MQ=255)
      acaactacaactATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATATACTACaaa                                                                                      <  1:2261291/91‑1 (MQ=255)
       caactacaactaTTAAAATAAGTGGAACCATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACaaaa                                                                                     >  2:1391945/1‑90 (MQ=255)
         actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACCCCTTATATAAGACTACATTTGTAGTATATTACAAATgt                                                                                   >  2:763508/1‑91 (MQ=255)
            acaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgtagt                                                                                <  2:1563789/91‑1 (MQ=255)
            acaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgtagt                                                                                >  1:468190/1‑91 (MQ=255)
                 taTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTa                                                                           <  1:3233105/91‑1 (MQ=255)
                  aTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTAt                                                                          >  1:1619673/1‑91 (MQ=255)
                   ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAATATTTATg                                                                         >  2:3935061/1‑91 (MQ=255)
                      aaaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTGCATTTGTAGTATATTACAAATGTAGTATTTATGTCa                                                                      >  2:1010888/1‑91 (MQ=255)
                      aaaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCa                                                                      >  2:1540465/1‑91 (MQ=255)
                      aaaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCa                                                                      <  1:2118424/91‑1 (MQ=255)
                       aaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCaa                                                                     >  2:1102251/1‑91 (MQ=255)
                       aaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCaa                                                                     >  2:1857123/1‑91 (MQ=255)
                          aaGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaat                                                                  >  1:1739175/1‑91 (MQ=255)
                          aaGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaat                                                                  >  2:1867693/1‑91 (MQ=255)
                          aaGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaat                                                                  >  2:3796835/1‑91 (MQ=255)
                            gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa                                                                >  2:2085621/1‑91 (MQ=255)
                             tGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataat                                                               >  1:2727425/1‑91 (MQ=255)
                              ggAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg                                                              <  2:2989414/91‑1 (MQ=255)
                               gAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGt                                                             >  1:2777103/1‑91 (MQ=255)
                               gAACAATTTTTATCTTTTTCACCAATATACTCCCTATATACGACTACATTTGTAGTATATTACCAATGTAGTATTTATGTCAAcataatcc                                                             >  2:3817441/1‑89 (MQ=255)
                                aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                            <  2:2798256/91‑1 (MQ=255)
                                aaCAAGTTTTATCTTTTTCATCAATATACTCCGTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                            <  2:3042173/91‑1 (MQ=255)
                                  cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat                                                          >  2:654540/1‑91 (MQ=255)
                                     tttttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAt                                                       >  2:338397/1‑91 (MQ=255)
                                      ttttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAtt                                                      >  2:520391/1‑91 (MQ=255)
                                        ttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAtttt                                                    >  2:1325163/1‑91 (MQ=255)
                                          aTCTTTTTTTTCTTTATACTCCTTATATAAGACTACATTTGTAGTATATTTCAATTGTATTATTTATGTCATACTAATTGTATAATATTtg                                                  >  2:1454559/1‑91 (MQ=255)
                                           tCTTTTTCATCAATATACTCCTTACATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                                 <  2:3687321/91‑1 (MQ=255)
                                              tttttCATCAATATACTCCTTATATAAGGCTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGa                                               <  2:2323584/91‑1 (MQ=255)
                                              tttttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGa                                               >  2:3226045/1‑91 (MQ=255)
                                                tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAATATAATGTTATAATTTTTGTGata                                             <  2:3574184/91‑1 (MQ=255)
                                                    atcaATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGa                                         >  1:1354535/1‑91 (MQ=255)
                                                    atcaATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGa                                         >  2:1583730/1‑91 (MQ=255)
                                                      caATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAgg                                       >  1:329698/1‑91 (MQ=255)
                                                       aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                      >  2:2720485/1‑91 (MQ=255)
                                                       aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                      >  2:1626156/1‑91 (MQ=255)
                                                        atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                     <  1:1213878/91‑1 (MQ=255)
                                                        atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                     >  1:1818216/1‑91 (MQ=255)
                                                          ataCTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTa                                   <  1:2690149/91‑1 (MQ=255)
                                                           taCTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTCATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAg                                  >  2:1812052/1‑91 (MQ=255)
                                                           taCTCCTTATATAAGACTACATTTGTAGTATATGACACATGTAGTACTTATGTCAAAATGATGTTATAATTTTTGTGATATGGAGGTGTAg                                  >  1:3158459/1‑91 (MQ=255)
                                                               ccTTATATAAGACTACATTTGTAGTATATTACAGATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAgcatg                              >  2:2519821/1‑87 (MQ=255)
                                                                 ttATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGgtg                            >  1:2814225/1‑91 (MQ=255)
                                                                  tatataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGgtgt                           <  1:791105/91‑1 (MQ=255)
                                                                   atataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGtt                          <  1:366210/91‑1 (MQ=255)
                                                                    tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                         <  2:1642747/91‑1 (MQ=255)
                                                                     ataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTat                        >  2:1844566/1‑91 (MQ=255)
                                                                     ataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTat                        >  1:1912281/1‑91 (MQ=255)
                                                                         taCTACATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc                    >  1:1429348/2‑91 (MQ=255)
                                                                         gACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc                    <  2:2483164/91‑1 (MQ=255)
                                                                           cTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtt                  >  2:212611/1‑91 (MQ=255)
                                                                           cTACATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtt                  >  2:1435358/1‑91 (MQ=255)
                                                                            tACATTTGTAGTATATTACAAATGTAGTATTTATGGCAAAGTAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCAGCttt                 <  2:579149/91‑1 (MQ=255)
                                                                            tACATTTGTAGTATATTACAAATGTAGAATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttt                 >  2:2027637/1‑91 (MQ=255)
                                                                             aCATTTGTAGTATATTGCAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt                <  2:3776301/91‑1 (MQ=255)
                                                                             aCATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt                <  1:1855412/91‑1 (MQ=255)
                                                                             aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt                >  1:3927201/1‑91 (MQ=255)
                                                                             aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt                >  1:28867/1‑91 (MQ=255)
                                                                             aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt                >  1:3435498/1‑91 (MQ=255)
                                                                              cATTTGTAGTATATTACAAATGTAGTCTTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttttt               >  2:2239689/1‑91 (MQ=255)
                                                                              cATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttttt               >  1:1042040/1‑91 (MQ=255)
                                                                              cATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttttt               <  1:820540/91‑1 (MQ=255)
                                                                              cATTTGTAGTATATTACAAATGTAGTATTTATATCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttttt               <  1:3371380/91‑1 (MQ=255)
                                                                               aTTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGCTATGGAGGTGTAGAAGGTGTTATCATCtttttt              >  2:932276/1‑91 (MQ=255)
                                                                                tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATACTGTTATAATTTTTGTGATAAGGAGGTGTAGAAGGTGTTATCATCtttctga             >  1:2008217/1‑87 (MQ=255)
                                                                                tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTa             >  1:642662/1‑91 (MQ=255)
                                                                                tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTa             <  1:1136271/91‑1 (MQ=255)
                                                                                tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTa             >  2:1087607/1‑91 (MQ=255)
                                                                                tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTa             >  2:1279035/1‑91 (MQ=255)
                                                                                 ttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGGTATGGAGGTGTAGAAGGTGTTATCATCTTTTTTaa            <  2:228753/91‑1 (MQ=255)
                                                                                 ttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTaa            >  2:3551043/1‑91 (MQ=255)
                                                                                  tGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAAt           >  2:3086715/1‑91 (MQ=255)
                                                                                    tagtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt         >  2:2761950/1‑91 (MQ=255)
                                                                                    taatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt         >  2:3460509/4‑91 (MQ=255)
                                                                                    taatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt         >  2:2422705/4‑91 (MQ=255)
                                                                                    taatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt         >  1:2496236/4‑91 (MQ=255)
                                                                                     agtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt        >  1:1000001/1‑91 (MQ=255)
                                                                                     agtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt        >  1:3107807/1‑91 (MQ=255)
                                                                                     aatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt        >  2:1158207/3‑91 (MQ=255)
                                                                                     aatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt        >  1:240148/3‑91 (MQ=255)
                                                                                     aatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt        >  2:118762/3‑91 (MQ=255)
                                                                                     aatataCTACAAACGTAGTATTTATGTCAAAATAATGTTATAATATTAGTGATATGGCGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt        >  2:412983/3‑91 (MQ=255)
                                                                                       tataTTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTaa      >  2:715648/1‑91 (MQ=255)
                                                                                       tataCTACAAATGTTGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTaa      >  2:589343/1‑91 (MQ=255)
                                                                                       tataCTACAAATGTTGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTaa      >  2:2622245/1‑91 (MQ=255)
                                                                                       tataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTaa      >  2:3204542/1‑91 (MQ=255)
                                                                                        atatTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAg     >  2:221507/1‑91 (MQ=255)
                                                                                        atatTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAg     <  1:3651551/91‑1 (MQ=255)
                                                                                        atatTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAg     <  1:602290/91‑1 (MQ=255)
                                                                                          atTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGta   >  1:564107/1‑91 (MQ=255)
                                                                                           |                                                                                          
AACCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGTAT  >  USA300TCH1516_ALE/41081‑41261

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