Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,354,972 G→A intergenic (‑168/‑154) USA300HOU_RS06735 ← / → USA300HOU_RS06740 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,354,9720GA99.4% 1372.7 / NA 358intergenic (‑168/‑154)USA300HOU_RS06735/USA300HOU_RS06740hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (170/185);  minor base C (2/0);  total (173/185)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 358 total aligned reads displayed.

ACATTTTTATTAGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAGATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTGAAGATT  >  USA300TCH1516_ALE/1354884‑1355058
                                                                                        |                                                                                      
  aTTTTTATTAGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTc                                                                                    <  1:3306825/91‑1 (MQ=255)
  aTTTTTATTAGTTGATTTTAGTTCATCATTATTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTc                                                                                    <  1:5069119/91‑1 (MQ=255)
  aTTTTTATTAGCTGATTTTAGTACATCATTGTTATAATCAAATTATAAACTGACAGATATTAATGTTCATTGAATATGACGTGAAAAATTc                                                                                    >  2:4630483/1‑91 (MQ=255)
      ttattaGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGa                                                                                >  1:2227217/1‑91 (MQ=255)
       tattaGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGaa                                                                               <  2:2212161/91‑1 (MQ=255)
       tattaGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGaa                                                                               >  1:2097905/1‑91 (MQ=255)
        attaGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAAt                                                                              >  2:4314591/1‑91 (MQ=255)
         ttaGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAAtt                                                                             >  2:3028117/1‑91 (MQ=255)
          taGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTc                                                                            <  2:6185898/91‑1 (MQ=255)
           aGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCa                                                                           <  1:1351131/91‑1 (MQ=255)
            gTTGATTTTAGTTCATCATTGTTATAATCGAATTATAAACTGATAGATATTGAAGTTCATTGAATATGACGTGAAACAATCGTTAATTCaa                                                                          >  2:3769141/1‑91 (MQ=255)
            gTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCaa                                                                          >  2:546731/1‑91 (MQ=255)
             ttGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAg                                                                         >  1:5628873/1‑91 (MQ=255)
              tGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGt                                                                        <  1:3281305/91‑1 (MQ=255)
               gATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGtt                                                                       <  1:6423007/91‑1 (MQ=255)
               gATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGtt                                                                       >  1:2960807/1‑91 (MQ=255)
               gATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGtt                                                                       <  1:3003098/91‑1 (MQ=255)
               gATTTTAGTTCATCAGTGTTATAATCAAATTATAAACTGACAGATCTTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTGAAGtt                                                                       <  2:6454231/91‑1 (MQ=255)
                aTTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGttt                                                                      <  1:761210/91‑1 (MQ=255)
                  tttAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTAt                                                                    <  2:6303613/91‑1 (MQ=255)
                   ttAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCACTGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATg                                                                   <  2:3885400/91‑1 (MQ=255)
                   ttAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATg                                                                   <  2:1240306/91‑1 (MQ=255)
                   ttAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATg                                                                   <  2:2453455/91‑1 (MQ=255)
                    tAGTTCATCATTGTTCTAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGt                                                                  <  2:5073795/91‑1 (MQ=255)
                    tAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGt                                                                  <  1:346356/91‑1 (MQ=255)
                    tAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGt                                                                  <  1:2386741/91‑1 (MQ=255)
                     aGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTc                                                                 >  2:1699061/1‑91 (MQ=255)
                      gTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCg                                                                >  2:2004327/1‑91 (MQ=255)
                      gTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCg                                                                >  2:568648/1‑91 (MQ=255)
                      gTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCg                                                                <  1:2048083/91‑1 (MQ=255)
                       ttCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGa                                                               <  1:5703184/91‑1 (MQ=255)
                         catcatTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAAt                                                             >  2:4025893/1‑91 (MQ=255)
                          atcatTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAAtt                                                            >  2:1045452/1‑91 (MQ=255)
                           tcatTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAAttt                                                           <  2:1001873/91‑1 (MQ=255)
                               tGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGt                                                       >  1:4744860/1‑91 (MQ=255)
                                  tataATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTtat                                                    >  1:1863491/1‑91 (MQ=255)
                                  tataATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTtat                                                    <  2:4482712/91‑1 (MQ=255)
                                   ataATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTtata                                                   <  2:843345/91‑1 (MQ=255)
                                   ataATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTtata                                                   <  2:2365770/91‑1 (MQ=255)
                                   ataATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTtata                                                   >  1:5841275/1‑91 (MQ=255)
                                    taATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATaa                                                  <  2:5136213/91‑1 (MQ=255)
                                     aaTCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGt                                                                  <  1:174959/74‑1 (MQ=255)
                                     aaTCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAAc                                                 <  2:3671525/91‑1 (MQ=255)
                                     aaTCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAAc                                                 <  2:5855956/91‑1 (MQ=255)
                                      aTCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACg                                                <  1:3514308/91‑1 (MQ=255)
                                       tCAAATTATAAACTGACAGATATTGATGTTCAATGAATGTGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACgg                                               <  1:5442770/91‑1 (MQ=255)
                                        cAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGt                                              <  1:3087628/91‑1 (MQ=255)
                                        cAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGt                                              >  1:4373507/1‑91 (MQ=255)
                                          aaTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCa                                            >  1:4468247/1‑91 (MQ=255)
                                          aaTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCa                                            <  2:5209018/91‑1 (MQ=255)
                                          aaTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCa                                            <  1:2004327/91‑1 (MQ=255)
                                          aaTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCa                                            >  2:1080522/1‑91 (MQ=255)
                                          aaTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCa                                            <  2:1174241/91‑1 (MQ=255)
                                          aaTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCa                                            >  2:2116156/1‑91 (MQ=255)
                                          aaTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCa                                            >  1:4072639/1‑91 (MQ=255)
                                           aTTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCAt                                           >  1:4388174/1‑91 (MQ=255)
                                            ttATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAAAGTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCAtt                                          <  1:4122476/91‑1 (MQ=255)
                                                aaaCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTaaa                                      >  1:1896526/1‑91 (MQ=255)
                                                aaaCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTaaa                                      >  1:5226568/1‑91 (MQ=255)
                                                aaaCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTaaa                                      >  1:4812517/1‑91 (MQ=255)
                                                aaaCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTaaa                                      >  1:4513693/1‑91 (MQ=255)
                                                 aaCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAAt                                     <  1:1220183/91‑1 (MQ=255)
                                                  aCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATg                                    <  1:611883/91‑1 (MQ=255)
                                                  aCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTGATGTTATAACGGTCATTTAAATg                                    <  2:5634379/91‑1 (MQ=255)
                                                  aCTGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTAATGTCAAATTTATGTTATAACGGTCATTTAAAtt                                    >  2:4265206/1‑90 (MQ=255)
                                                    tGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGAc                                  <  1:867559/91‑1 (MQ=255)
                                                    tGACAGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGAc                                  <  1:6015601/91‑1 (MQ=255)
                                                        aGATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGaa                              <  1:5571316/91‑1 (MQ=255)
                                                         gATATTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAAt                             >  1:4548332/1‑91 (MQ=255)
                                                           tatTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTa                           <  1:2761159/91‑1 (MQ=255)
                                                            atTGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAg                          <  2:3327786/91‑1 (MQ=255)
                                                             ttGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAgg                         <  2:3487533/91‑1 (MQ=255)
                                                             ttGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAgg                         <  1:2607197/91‑1 (MQ=255)
                                                              tGATGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGt                        >  2:2617449/1‑91 (MQ=255)
                                                                 tGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCAc                     <  2:5436905/91‑1 (MQ=255)
                                                                 tGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCAc                     >  2:5913486/1‑91 (MQ=255)
                                                                 tGTTCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCAc                     >  2:2891856/1‑91 (MQ=255)
                                                                   ttCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTc                   <  2:6065858/91‑1 (MQ=255)
                                                                   ttCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTc                   <  2:1948901/91‑1 (MQ=255)
                                                                    tCAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATTACGGTCATTTAAATGACAGAATTAGGTCACTCa                  <  1:6531836/91‑1 (MQ=255)
                                                                     cAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCAt                 <  2:5293609/91‑1 (MQ=255)
                                                                     cAATGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCAt                 <  1:3826828/91‑1 (MQ=255)
                                                                      aaTGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATa                <  2:2794085/91‑1 (MQ=255)
                                                                       aTGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAg               <  2:3101423/91‑1 (MQ=255)
                                                                       aTGAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAg               <  1:3606024/91‑1 (MQ=255)
                                                                         gAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTa             >  1:904241/1‑91 (MQ=255)
                                                                         gAATATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTa             >  2:4874033/1‑91 (MQ=255)
                                                                          aaTATGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTAt            <  1:3498688/91‑1 (MQ=255)
                                                                           atatGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTAtt           >  1:6334584/1‑91 (MQ=255)
                                                                           atatGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACCCATAGTAtt           >  1:1676954/1‑91 (MQ=255)
                                                                            tatGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTAttt          >  2:6513034/1‑91 (MQ=255)
                                                                             atGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTAtttt         >  2:3377171/1‑91 (MQ=255)
                                                                              tGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTg        >  2:5550487/1‑91 (MQ=255)
                                                                              tGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTg        <  2:4545730/91‑1 (MQ=255)
                                                                              tGACGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTg        <  2:3175707/91‑1 (MQ=255)
                                                                                aCGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTGaa      >  2:1549322/1‑91 (MQ=255)
                                                                                 cGTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTGAAg     <  2:4766274/91‑1 (MQ=255)
                                                                                  gTGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTGAAGa    <  1:2897545/91‑1 (MQ=255)
                                                                                   tGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTGAAGAt   >  1:3653184/1‑91 (MQ=255)
                                                                                   tGAAAAATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTGAAGAt   >  1:2457693/1‑91 (MQ=255)
                                                                                        |                                                                                      
ACATTTTTATTAGTTGATTTTAGTTCATCATTGTTATAATCAAATTATAAACTGACAGATATTGATGTTCAATGAATATGACGTGAAAGATTCGTGAATTCAAGTTTATGTCGAATTTATGTTATAACGGTCATTTAAATGACAGAATTAGGTCACTCATAGTATTTTGAAGATT  >  USA300TCH1516_ALE/1354884‑1355058

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: