Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,133,490 A→T 5.6% intergenic (+334/+38) dapE → / ← USA300HOU_RS10885 putative succinyl‑diaminopimelate desuccinylase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,133,4900AT 5.6% 1047.7 / 37.3 307intergenic (+334/+38)dapE/USA300HOU_RS10885putative succinyl‑diaminopimelate desuccinylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (160/129);  new base T (6/11);  total (166/141)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 322 total aligned reads displayed.

AAAGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGT  >  USA300TCH1516_ALE/2133400‑2133580
                                                                                          |                                                                                           
aaaGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAaata                                                                                             >  1:3414365/1‑91 (MQ=255)
aaaGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAaata                                                                                             >  1:1327507/1‑91 (MQ=255)
  aGGAACGGGACATAAATCCCTACAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAaataat                                                                                           <  2:3196641/91‑1 (MQ=255)
   ggAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAAtt                                                                                          >  1:2112150/1‑91 (MQ=255)
   ggAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAAtt                                                                                          <  2:2388160/91‑1 (MQ=255)
    gAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAAttt                                                                                         >  1:3133174/1‑91 (MQ=255)
    gAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAAttt                                                                                         >  1:2127797/1‑91 (MQ=255)
      aTGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtat                                                                                      >  1:1519022/1‑91 (MQ=255)
      aTGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtat                                                                                      >  2:3625381/1‑91 (MQ=255)
      aTGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTAtttgat                                                                                      <  1:3880461/91‑4 (MQ=25)
      aTGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTAtttgat                                                                                      <  1:1790172/91‑4 (MQ=25)
        gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                    >  2:3456360/1‑91 (MQ=255)
         ggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGa                                                                                   >  1:551711/1‑91 (MQ=255)
         ggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGa                                                                                   <  1:3274533/91‑1 (MQ=25)
          gACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGaa                                                                                  <  1:689935/91‑1 (MQ=255)
            cATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAAt                                                                                >  1:1556409/1‑91 (MQ=255)
            cATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAAt                                                                                <  1:2971966/91‑1 (MQ=255)
            cATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAAt                                                                                <  2:2252406/91‑1 (MQ=255)
             aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATc                                                                               <  2:1834359/91‑1 (MQ=21)
                 aTCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCt                                                                           >  1:3222158/1‑91 (MQ=255)
                 aTCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCt                                                                           >  1:2198722/1‑91 (MQ=255)
                  tCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCtt                                                                          >  1:1719401/1‑91 (MQ=255)
                  tCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCtt                                                                          >  2:2495690/1‑91 (MQ=255)
                  tCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATCAGCtt                                                                          <  1:1947056/91‑1 (MQ=21)
                  acccTAGAAAACAGCAGCTCAACGTGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCtt                                                                          <  2:962192/90‑1 (MQ=255)
                        gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                    >  1:2931472/1‑91 (MQ=255)
                        gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                    >  1:2233744/1‑91 (MQ=255)
                        gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                    >  2:3294760/1‑91 (MQ=255)
                        gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                    <  2:134270/91‑1 (MQ=255)
                         aaaaCAGCAGCTCAACGCGCTGTAAATTATCAAAGACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACa                                                                   <  2:2087210/91‑1 (MQ=255)
                         aaaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACa                                                                   >  2:3341545/1‑91 (MQ=255)
                           aaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACAtt                                                                 <  1:3700426/91‑1 (MQ=255)
                            aCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTg                                                                <  2:1962326/91‑1 (MQ=255)
                            aCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTg                                                                <  2:1999922/91‑1 (MQ=255)
                             cagcagCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGt                                                               >  1:1596567/1‑91 (MQ=255)
                             cagcagCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGt                                                               >  1:1476033/1‑91 (MQ=255)
                              agcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTc                                                              >  1:3801317/1‑91 (MQ=255)
                               gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                             >  2:1733588/1‑91 (MQ=255)
                               gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                             >  2:1611096/1‑91 (MQ=255)
                               gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                             <  2:1060911/91‑1 (MQ=255)
                               gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                             >  1:3635431/1‑91 (MQ=255)
                                cagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTa                                                            >  2:2791298/1‑91 (MQ=255)
                                cagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTa                                                            <  2:2910817/91‑1 (MQ=255)
                                cagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTa                                                            >  1:2486554/1‑91 (MQ=255)
                                     cgttgAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                       >  1:1839700/5‑91 (MQ=255)
                                     cAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                       <  2:844244/91‑1 (MQ=255)
                                       ccGAGCTGTAAATTATCAAAAACATTTGAAAAAAGCTGTTTTCTATCAAATAATTTATAGAAATCAGCTATTTTACATTGTCTAAGAACtt                                                     <  2:2711492/90‑1 (MQ=17)
                                       aCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACtt                                                     <  1:3084406/91‑1 (MQ=255)
                                       aCGAGCTGTAAATTATCAAATACATTGGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCGTTTTTACATTGTCTAAGAACtt                                                     <  2:306429/91‑1 (MQ=255)
                                        tgAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTa                                                    >  2:1064753/2‑91 (MQ=255)
                                         gAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTaa                                                   <  1:1120310/91‑1 (MQ=255)
                                          aGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAt                                                  <  1:226863/91‑1 (MQ=255)
                                           gCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtg                                                 <  2:605935/91‑1 (MQ=255)
                                           gCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtg                                                 >  1:1049093/1‑91 (MQ=255)
                                            cTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTGAtgt                                                >  1:3299041/1‑91 (MQ=255)
                                            cTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgt                                                >  2:971393/1‑91 (MQ=255)
                                               tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                             <  1:1373893/91‑1 (MQ=255)
                                                 aaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctg                                           >  2:1513894/1‑91 (MQ=255)
                                                  aTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgt                                          >  2:3507103/1‑91 (MQ=255)
                                                   ttATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtc                                         <  1:2013882/91‑1 (MQ=255)
                                                     aTCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtctc                                       >  2:1256588/1‑91 (MQ=255)
                                                     aTCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtctc                                       <  2:3767843/91‑1 (MQ=255)
                                                     aTCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtctc                                       <  1:3764487/91‑1 (MQ=255)
                                                      tCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATgtctctgtctct                                      >  2:2144220/1‑91 (MQ=255)
                                                        aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                    >  1:61150/1‑91 (MQ=255)
                                                        aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                    >  2:2969114/1‑91 (MQ=255)
                                                        aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                    >  2:1488310/1‑91 (MQ=255)
                                                         aaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTg                                   >  1:2323483/1‑91 (MQ=255)
                                                         aaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTCTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTg                                   >  2:3904821/1‑91 (MQ=255)
                                                           tACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGa                                 <  1:2622998/91‑1 (MQ=255)
                                                              aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                              >  1:2754370/1‑91 (MQ=255)
                                                               tttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACtt                             >  1:3471762/1‑91 (MQ=255)
                                                               tttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATCGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACtt                             >  2:3573269/1‑91 (MQ=255)
                                                                ttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACttt                            >  2:348777/1‑91 (MQ=255)
                                                                 tGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTta                           >  1:334619/1‑91 (MQ=255)
                                                                  gAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtat                          <  1:753110/91‑1 (MQ=255)
                                                                  gAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtat                          >  1:3528345/1‑91 (MQ=255)
                                                                  gAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtat                          >  2:714766/1‑91 (MQ=255)
                                                                   aaaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtata                         >  1:3634684/1‑91 (MQ=255)
                                                                    aaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTatat                        >  2:1890422/1‑91 (MQ=255)
                                                                     aaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATAtt                       <  1:1991404/91‑1 (MQ=255)
                                                                       aaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTc                     >  2:3937115/1‑91 (MQ=255)
                                                                       aaGCTGATTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTc                     >  1:1264383/1‑91 (MQ=21)
                                                                            gATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTACGAACTTAATGTCTCTGTCTCTTTGGAACTATATATTTCATAAg                >  1:3954696/1‑91 (MQ=255)
                                                                             aTTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt               <  1:477969/91‑1 (MQ=255)
                                                                              tttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTa              >  1:389632/1‑91 (MQ=255)
                                                                               ttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTAc             <  2:529879/91‑1 (MQ=255)
                                                                               ttCGATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGATCTTAATGTCTCTGTCTCTTTGGATCTTTATATTTCATAAGTAc             <  1:1429200/91‑1 (MQ=255)
                                                                                tCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACt            <  1:969087/91‑1 (MQ=255)
                                                                                 cTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTa           <  1:2448606/91‑1 (MQ=255)
                                                                                  tATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAt          <  2:3160512/91‑1 (MQ=255)
                                                                                   tataaaTTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATCTTTCATAAGTACTAt          >  2:417710/4‑91 (MQ=25)
                                                                                    tCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAttt        >  2:3253905/1‑91 (MQ=255)
                                                                                    tCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAttt        <  2:2565390/91‑1 (MQ=255)
                                                                                     cAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAtttt       >  2:3408277/1‑91 (MQ=255)
                                                                                      aaaTAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTa      >  2:3838575/1‑91 (MQ=255)
                                                                                      aaaTAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTa      >  2:819239/1‑91 (MQ=255)
                                                                                       aataatTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTAt     >  1:101144/1‑91 (MQ=255)
                                                                                        ataatTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGCCTCTTTGGAACTTTATATTTCATAAGTACTATTTTAtt    >  2:457520/1‑91 (MQ=255)
                                                                                         taatTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTg   >  1:1168860/1‑91 (MQ=255)
                                                                                          |                                                                                           
AAAGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGT  >  USA300TCH1516_ALE/2133400‑2133580

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: