Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,430,685 Δ1 bp 11.6% intergenic (‑109/‑276) suhB_2 ← / → birA_2 Inositol‑1‑monophosphatase/Bifunctional ligase/repressor BirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,430,6850T.11.6% 758.4 / 87.9 242intergenic (‑109/‑276)suhB_2/birA_2Inositol‑1‑monophosphatase/Bifunctional ligase/repressor BirA
Reads supporting (aligned to +/- strand):  ref base T (92/121);  new base . (2/26);  total (94/148)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

TATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTAT  >  USA300TCH1516_ALE/2430595‑2430775
                                                                                          |                                                                                           
tatcattatcatTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGt                                                                                             >  1:193048/1‑91 (MQ=255)
 atcattatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGta                                                                                            >  2:3865401/1‑90 (MQ=255)
 atcattatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTc                                                                                            >  2:752049/1‑91 (MQ=255)
  tcattatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCa                                                                                           <  1:3309204/91‑1 (MQ=255)
  tcattatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCa                                                                                           <  1:2890740/91‑1 (MQ=255)
   cattatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCAc                                                                                          >  2:1940329/1‑91 (MQ=255)
   cattatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCAATATATGCATAAAAAAATCGGCTAGTCAc                                                                                          <  2:590770/91‑1 (MQ=255)
     ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGt                                                                                        >  1:3396877/1‑91 (MQ=255)
     ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGt                                                                                        >  2:2887862/1‑91 (MQ=255)
     ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTg                                                                                       <  2:641295/91‑1 (MQ=255)
     ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTg                                                                                       <  1:920912/91‑1 (MQ=255)
     ttagcattTTATGATCGCCCGATATTAAAGGAATCAACAATACAGCTGACGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGt                                                                                        <  1:2276706/91‑1 (MQ=255)
      tatcattTTATTATTACCCTATATTAAAAGATTCAACAATATAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTga                                                                                      <  2:3480229/91‑2 (MQ=255)
      tatcattTTATTATTACCCTATATTAAAAGAATCAACAATACCACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTg                                                                                       <  2:1858428/91‑1 (MQ=255)
      tatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTg                                                                                       <  1:3506381/91‑1 (MQ=255)
      tatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTg                                                                                       <  2:2343757/91‑1 (MQ=255)
      tatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTg                                                                                       <  1:2135208/91‑1 (MQ=255)
        tcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTgact                                                                                    <  1:882797/91‑4 (MQ=255)
        tcattTAATTATTACCCTATAGTAAAAGAATCAACAATACAACTGAAGCCTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTgact                                                                                    <  2:1707085/91‑4 (MQ=255)
          attTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAg                                                                                  >  2:2510066/1‑91 (MQ=255)
          attTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAATCGGCTAGTCACGTGCTAGc                                                                                 >  1:2793939/1‑91 (MQ=255)
           ttttATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGc                                                                                 >  1:935582/1‑91 (MQ=255)
            tttATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGcc                                                                                >  2:2796561/1‑91 (MQ=255)
            tttATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGcc                                                                                <  2:2943725/91‑1 (MQ=255)
             ttattattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTGACTAGCCg                                                                               <  1:2993700/91‑1 (MQ=255)
              tattattaCCCTATATTGAAAGAATCAACAATACAACGGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTGACTAGCCGa                                                                              <  2:413033/91‑1 (MQ=255)
              tattattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGa                                                                              >  1:3851684/1‑91 (MQ=255)
              tattattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTGACTAGCCGa                                                                              <  1:427403/91‑1 (MQ=255)
               attattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGAc                                                                             <  2:2844353/91‑1 (MQ=255)
                ttattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACa                                                                            <  1:2393644/91‑1 (MQ=255)
                  attaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACaaa                                                                          <  1:517919/91‑1 (MQ=255)
                   ttaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAAt                                                                         >  1:1537799/1‑91 (MQ=255)
                   ttaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAAt                                                                         >  1:2615108/1‑91 (MQ=255)
                    taCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGACTAGTCACGTGCTAGCCGACAAATa                                                                        >  2:3458292/1‑91 (MQ=255)
                     aCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAg                                                                       >  1:1352114/1‑91 (MQ=255)
                      cccTATATTAAAAGAATCAACAATACAGCGGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAga                                                                      <  2:1708306/91‑1 (MQ=255)
                      cccTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAga                                                                      <  1:1010817/91‑1 (MQ=255)
                      cccTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCAGTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAga                                                                      <  2:3271763/91‑1 (MQ=255)
                           tatTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaagga                                                                 <  2:3012771/91‑1 (MQ=255)
                           tatTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaagga                                                                 <  2:1023756/91‑1 (MQ=255)
                             ttAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaaggaaa                                                               >  2:314356/1‑91 (MQ=255)
                             ttAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaaggaaa                                                               >  1:2461897/1‑91 (MQ=255)
                             ttAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaaggaaa                                                               <  1:3188849/91‑1 (MQ=255)
                                 aaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAaagtaa                                                           >  2:3655138/1‑91 (MQ=255)
                                 aaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAaagtaa                                                           >  2:2475667/1‑91 (MQ=255)
                                 aaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCGAGTCACGTGCTAGCCGACAAATAGAAAGGAaagtac                                                           >  2:3915600/1‑90 (MQ=255)
                                    aaTCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAagtaagta                                                        <  1:1416841/91‑1 (MQ=255)
                                    aaTCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAagtaagta                                                        <  2:3779706/91‑1 (MQ=255)
                                      tCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaat                                                      >  1:3656824/1‑91 (MQ=255)
                                      tCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaat                                                      <  1:1238239/91‑1 (MQ=255)
                                      tCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaat                                                      >  1:2546948/1‑91 (MQ=255)
                                      tCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaat                                                      <  1:797387/91‑1 (MQ=255)
                                       caacaaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaata                                                     <  1:1131828/91‑1 (MQ=255)
                                       caacaaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaata                                                     <  2:2835710/91‑1 (MQ=255)
                                       caacaaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaata                                                     <  2:2998681/91‑1 (MQ=255)
                                         acaaTACAACTGAAGACTTCCTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTCGCCGACAAAGAGAAAGGAAAGTAAGTaataaa                                                   >  2:211857/1‑91 (MQ=255)
                                            aTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAatat                                                >  2:1460287/1‑91 (MQ=255)
                                             tACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATAtt                                               <  2:1177284/91‑1 (MQ=255)
                                             tACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATAtt                                               >  1:3088303/1‑91 (MQ=255)
                                             tACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATAtt                                               >  1:2966998/1‑91 (MQ=255)
                                               cAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGa                                             <  1:2742239/91‑1 (MQ=255)
                                                aaCTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGaa                                            <  2:3142920/91‑1 (MQ=255)
                                                 aCTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAg                                           >  2:2205898/1‑91 (MQ=255)
                                                 aCTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAg                                           <  1:2044593/91‑1 (MQ=255)
                                                  cTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGCAAGTAAGTAATAAATATTGAAGa                                          >  2:721570/1‑91 (MQ=255)
                                                  cTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGa                                          >  1:353606/1‑91 (MQ=255)
                                                   tGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAt                                         >  2:2250687/1‑91 (MQ=255)
                                                     aaGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgt                                       <  1:3348169/91‑1 (MQ=255)
                                                      aGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgtt                                      <  2:1615471/91‑1 (MQ=255)
                                                      aGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgtt                                      <  1:1068196/91‑1 (MQ=255)
                                                           tcttcATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGAt                                 <  1:687713/91‑1 (MQ=255)
                                                            cttcATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATg                                <  1:2904139/91‑1 (MQ=255)
                                                            cttcATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATg                                >  1:1810837/1‑91 (MQ=255)
                                                              tcATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTa                              >  2:449600/1‑91 (MQ=255)
                                                                 ttttATGCATAAAAAAATCGGCTAGTCCCGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGCAGATGTTGTGATGTAACt                           >  2:3845676/1‑91 (MQ=255)
                                                                 ttttATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACt                           <  1:45620/91‑1 (MQ=255)
                                                                     aTGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGaa                       <  2:295596/91‑1 (MQ=255)
                                                                      tGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                      <  2:2717785/91‑1 (MQ=255)
                                                                      tGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                      <  1:3000857/91‑1 (MQ=255)
                                                                      tGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                      <  1:1010073/91‑1 (MQ=255)
                                                                      tGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                      <  2:838419/91‑1 (MQ=255)
                                                                      tGCATAAAAAAATCGACTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAACAAATATTGAAGATGTTGTGATGTAACTTGAAc                      >  1:3365647/1‑91 (MQ=255)
                                                                        cATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGa                    >  2:3325318/1‑91 (MQ=255)
                                                                         aTAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGAt                   >  1:32285/1‑91 (MQ=255)
                                                                          tAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGAtt                  <  1:1896163/91‑1 (MQ=255)
                                                                          tAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGAtt                  <  2:2013987/91‑1 (MQ=255)
                                                                              aaaaTCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTaaaa              >  1:3910940/1‑91 (MQ=255)
                                                                               aaaTCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAg             <  1:3933286/91‑1 (MQ=255)
                                                                               aaaTCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAg             <  1:3727706/91‑1 (MQ=255)
                                                                                 aTCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCt           >  2:1409898/1‑91 (MQ=255)
                                                                                   cGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTAt         <  2:3060960/91‑1 (MQ=255)
                                                                                   cGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTAt         >  2:3960612/1‑91 (MQ=255)
                                                                                   cGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTAt         <  2:868892/91‑1 (MQ=255)
                                                                                   cGGCTAGCACGTGACTAGCCGATAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTAt         >  1:2251183/1‑91 (MQ=255)
                                                                                      cTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTg      >  1:148767/1‑91 (MQ=255)
                                                                                      cTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTg      >  2:2880211/1‑91 (MQ=255)
                                                                                      cTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTg      >  2:1541867/1‑91 (MQ=255)
                                                                                       tAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGt     <  1:3456762/91‑1 (MQ=255)
                                                                                       tAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGt     >  1:2147802/1‑91 (MQ=255)
                                                                                        aGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGtt    >  2:2288446/1‑91 (MQ=255)
                                                                                          |                                                                                           
TATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTAT  >  USA300TCH1516_ALE/2430595‑2430775

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