Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,840 T→G 13.5% intergenic (+34/+89) bcr_2 → / ← USA300HOU_RS13190 Bicyclomycin resistance protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,559,8400TG13.5% 511.2 / 64.8 171intergenic (+34/+89)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (73/74);  new base G (9/14);  total (83/88)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

Only 100 of 198 total aligned reads displayed.

TTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAT  >  USA300TCH1516_ALE/2559750‑2559929
                                                                                          |                                                                                         
ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                           <  2:1179450/91‑1 (MQ=255)
ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                           >  2:823324/1‑91 (MQ=255)
ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                           >  2:449481/1‑91 (MQ=255)
 tAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATGACACTTCAACATAATTAGAACCCTAGCAAAGAtatc                                                                                          >  2:3828972/1‑91 (MQ=255)
 tAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatc                                                                                          >  2:3767551/1‑91 (MQ=255)
 tAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatc                                                                                          >  2:3792173/1‑91 (MQ=255)
  agtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatct                                                                                         >  2:1745143/1‑91 (MQ=255)
  agtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatct                                                                                         <  2:2300637/91‑1 (MQ=255)
   gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                        <  1:1532573/91‑1 (MQ=255)
   gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                        <  1:2740320/91‑1 (MQ=255)
   gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                        <  2:1213478/91‑1 (MQ=255)
   gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                        <  2:1707610/91‑1 (MQ=255)
    tcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAAAtatctat                                                                                       >  2:3363574/1‑91 (MQ=255)
       gtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCtt                                                                                    <  1:2317674/91‑1 (MQ=255)
          tACAAATCATTTACTTTAAAATGATTAAAAAGCGACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                                 <  1:3181304/91‑2 (MQ=37)
           aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                <  2:2799738/91‑1 (MQ=255)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                               <  1:1898783/91‑4 (MQ=37)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                               <  2:2445609/91‑4 (MQ=37)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                               <  1:131816/91‑4 (MQ=37)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                               <  2:1988244/91‑4 (MQ=37)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                               >  2:389742/1‑91 (MQ=255)
             aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                              >  2:1601247/1‑91 (MQ=255)
             aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                              <  1:2452438/91‑1 (MQ=255)
               aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                            >  2:2258036/1‑91 (MQ=255)
               aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                            <  2:3413379/91‑1 (MQ=255)
                tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                <  1:2301718/86‑1 (MQ=255)
                tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                           <  1:1302342/91‑1 (MQ=255)
                 cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACCTAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGtt                                                                          <  1:3370186/91‑1 (MQ=21)
                 cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGtt                                                                          >  1:1794935/1‑91 (MQ=255)
                   tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttct                                                                        <  2:2654795/91‑1 (MQ=21)
                   tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttct                                                                        >  1:846819/1‑91 (MQ=255)
                    ttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttctt                                                                       >  1:830057/1‑91 (MQ=255)
                     tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                      >  1:2089875/1‑91 (MQ=255)
                        tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                   <  1:2041042/91‑1 (MQ=255)
                        tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                   >  1:1200937/1‑91 (MQ=255)
                        tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAGTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                   <  2:3170456/91‑1 (MQ=255)
                          tAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTAt                                                                 >  2:1415794/1‑91 (MQ=255)
                             aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                              >  2:1880978/1‑91 (MQ=255)
                             aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                              <  2:2083403/91‑1 (MQ=255)
                               tGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAAtt                                                            >  2:2199187/1‑91 (MQ=255)
                                gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                           >  1:1616999/1‑91 (MQ=255)
                                gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                           <  1:218980/91‑1 (MQ=255)
                                gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                           <  1:2843663/91‑1 (MQ=255)
                                 aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                          >  1:1877137/1‑91 (MQ=255)
                                 aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                          >  2:211613/1‑91 (MQ=255)
                                  ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                         <  1:1569733/91‑1 (MQ=255)
                                  ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                         <  1:721136/91‑1 (MQ=255)
                                   tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                                        >  1:2894949/1‑91 (MQ=255)
                                     aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                      <  1:1031882/91‑1 (MQ=255)
                                         gCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                                  <  1:1914124/91‑1 (MQ=255)
                                            aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                               >  2:3205472/1‑91 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                             <  2:2608569/91‑1 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                             >  1:323327/1‑91 (MQ=255)
                                                  cGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                         <  1:3771826/91‑1 (MQ=255)
                                                    cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                                       <  2:3734548/91‑1 (MQ=255)
                                                    cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                                       <  1:3458098/91‑1 (MQ=255)
                                                     aTAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaa                                      >  2:2429809/1‑91 (MQ=255)
                                                     aTAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaa                                      >  1:2965607/1‑91 (MQ=255)
                                                      tAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaa                                     <  2:1576757/91‑1 (MQ=255)
                                                          acTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                                 <  1:2073301/91‑1 (MQ=255)
                                                            ttCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACg                               >  2:1875919/1‑91 (MQ=255)
                                                             tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                              <  1:2137577/91‑1 (MQ=255)
                                                             tCAACATAATTAGAACCCTAGCAAAAATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                              >  2:273199/1‑91 (MQ=255)
                                                              cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                             >  2:3567056/1‑91 (MQ=255)
                                                                aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                           <  2:2341279/91‑1 (MQ=255)
                                                                 cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTc                          >  2:3507781/1‑91 (MQ=255)
                                                                 cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTc                          <  2:1214778/91‑1 (MQ=255)
                                                                  aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCa                         >  1:28132/1‑91 (MQ=255)
                                                                  aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCa                         <  1:3693693/91‑1 (MQ=255)
                                                                     aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                      >  1:2369220/1‑91 (MQ=255)
                                                                      ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                     >  1:2926884/1‑91 (MQ=255)
                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                    <  1:543150/91‑1 (MQ=255)
                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                    <  1:3793868/91‑1 (MQ=255)
                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                    <  1:543556/91‑1 (MQ=255)
                                                                        aGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                   >  1:3205959/1‑91 (MQ=21)
                                                                        aGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCTAATCTTCAACTAAAATTACGCCTTCATAGCaa                   >  1:3300768/1‑91 (MQ=255)
                                                                        aGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                   >  2:845037/1‑91 (MQ=255)
                                                                         gAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                  >  1:188455/1‑91 (MQ=21)
                                                                         gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                  <  1:595875/91‑1 (MQ=255)
                                                                         gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                  <  1:3736695/91‑1 (MQ=255)
                                                                         gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                  >  2:3473813/1‑91 (MQ=255)
                                                                          aaCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                 >  1:31005/1‑91 (MQ=21)
                                                                           aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                <  1:2014648/91‑1 (MQ=255)
                                                                           aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                <  1:2551231/91‑1 (MQ=255)
                                                                           aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                >  1:2714702/1‑91 (MQ=255)
                                                                             ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa              >  2:2177838/1‑91 (MQ=255)
                                                                             ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa              >  1:444454/1‑91 (MQ=255)
                                                                              ctgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt             >  1:3821609/5‑91 (MQ=25)
                                                                              ctgacAAAGATACATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt             >  1:2991829/5‑91 (MQ=37)
                                                                              cTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt             <  1:1031586/91‑1 (MQ=255)
                                                                               tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt            >  1:988217/1‑91 (MQ=255)
                                                                                aGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt           <  1:3146948/91‑1 (MQ=255)
                                                                                  cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt         <  2:692838/91‑1 (MQ=255)
                                                                                  cAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt         >  2:300636/1‑91 (MQ=38)
                                                                                    aaGATATCTATCTTTGTCAGGGATCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAt       <  2:3125617/91‑1 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa   <  1:3558778/91‑1 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa   >  1:1762199/1‑91 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa   >  1:1418947/1‑91 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa   >  1:1139384/1‑91 (MQ=255)
                                                                                         atctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAt  >  1:753920/1‑91 (MQ=255)
                                                                                          |                                                                                         
TTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAT  >  USA300TCH1516_ALE/2559750‑2559929

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: