Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,133,490 A→T 6.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 6.6% 806.2 / 33.4 242intergenic (+334/+38)dapE/USA300HOU_RS10885putative succinyl‑diaminopimelate desuccinylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (103/123);  new base T (6/10);  total (109/133)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 244 total aligned reads displayed.

GGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGT  >  USA300TCH1516_ALE/2133403‑2133580
                                                                                       |                                                                                           
ggAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAAtt                                                                                          <  2:52598/91‑1 (MQ=255)
 gAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAAttt                                                                                         >  2:2599214/1‑91 (MQ=255)
   aTGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtat                                                                                      >  1:3178463/1‑91 (MQ=255)
    tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtata                                                                                     >  1:2255772/1‑91 (MQ=255)
    tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtata                                                                                     >  2:2315710/1‑91 (MQ=255)
    tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtata                                                                                     >  1:2082256/1‑91 (MQ=255)
    tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtata                                                                                     <  1:470024/91‑1 (MQ=255)
    tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTAtttgata                                                                                     <  2:1681813/91‑5 (MQ=25)
    tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGACTTCTATCAAATAATTtata                                                                                     >  1:1515520/1‑91 (MQ=255)
     gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                    <  2:3429259/91‑1 (MQ=255)
     gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                    >  1:3666968/1‑91 (MQ=255)
     gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                    <  1:1583288/91‑1 (MQ=255)
      ggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGa                                                                                   <  1:2830459/91‑1 (MQ=25)
       gACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGaa                                                                                  >  1:2416382/1‑91 (MQ=255)
        aCATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGaaa                                                                                 >  2:3473543/1‑91 (MQ=255)
         cATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAAt                                                                                <  1:3238982/91‑1 (MQ=21)
          aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATc                                                                               <  1:2704389/91‑1 (MQ=255)
          aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATc                                                                               <  2:2909179/91‑1 (MQ=255)
          aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATc                                                                               <  2:2398404/91‑1 (MQ=21)
           tAAATCCCTAGAAAACAGCAGCTCAACGCGCTGTAAATTATCAAATACATGTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATCa                                                                              <  2:759044/91‑1 (MQ=18)
           tAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCa                                                                              <  2:1909206/91‑1 (MQ=255)
           tAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCa                                                                              <  1:304816/91‑1 (MQ=255)
           tAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATCa                                                                              <  1:1072447/91‑1 (MQ=21)
             aaTCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGc                                                                            <  1:684145/91‑1 (MQ=255)
             aaTCCCTAGAAAACAGCAGCTCAACGAGCTGTAAAGTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATCAGc                                                                            <  2:2672034/91‑1 (MQ=18)
              aTCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCGGCt                                                                           <  2:899888/91‑1 (MQ=255)
                cccTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttt                                                                         >  2:1767435/1‑91 (MQ=255)
                     gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                    >  1:2582118/1‑91 (MQ=255)
                     gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                    >  1:541648/1‑91 (MQ=255)
                      aaaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACa                                                                   >  2:2464972/1‑91 (MQ=255)
                       aaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACAt                                                                  <  1:951636/91‑1 (MQ=255)
                        aaCAGCAGCTCGACGAGCTGTAACTTATCAAATACATTTGAAAAAAGCTGATTTCTATCGAATAATTTATAGAAATCAGCTTTTTTACAtt                                                                 <  2:2967026/91‑1 (MQ=255)
                        aaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACAtt                                                                 >  2:1265663/1‑91 (MQ=255)
                          cagcagCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGt                                                               <  2:3516526/91‑1 (MQ=255)
                           agcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTc                                                              >  1:911465/1‑91 (MQ=255)
                            gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                             >  1:186584/1‑91 (MQ=255)
                              agcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTaa                                                           <  2:3734502/91‑1 (MQ=255)
                              agcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTaa                                                           <  2:3771095/91‑1 (MQ=255)
                                cTCAACGGGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGa                                                         <  2:2840897/91‑1 (MQ=255)
                                 tCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGaa                                                        <  1:3746410/91‑1 (MQ=255)
                                  cgttgAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                       >  1:3219352/5‑91 (MQ=255)
                                  cgttgAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                       >  2:3087802/5‑91 (MQ=255)
                                  cAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                       <  1:434910/91‑1 (MQ=255)
                                  cAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                       <  1:96868/91‑1 (MQ=255)
                                    aCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACtt                                                     <  1:756344/91‑1 (MQ=255)
                                    aCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACtt                                                     <  1:3835214/91‑1 (MQ=255)
                                     cGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                       <  2:1749743/88‑1 (MQ=255)
                                      gAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTaa                                                   <  2:323524/91‑1 (MQ=255)
                                         cTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgt                                                >  2:686339/1‑91 (MQ=255)
                                          tGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtc                                               >  2:2278207/1‑91 (MQ=255)
                                          tGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtc                                               <  2:2635442/91‑1 (MQ=255)
                                          tGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtc                                               <  1:1057984/91‑1 (MQ=255)
                                            tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                             >  2:2802889/1‑91 (MQ=255)
                                            tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                             >  1:3495161/1‑91 (MQ=255)
                                            tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                             <  2:444732/91‑1 (MQ=255)
                                            tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                             >  1:2107530/1‑91 (MQ=255)
                                            tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                             >  1:423334/1‑91 (MQ=255)
                                             aaaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCtttttt                                                                      >  1:2456615/1‑65 (MQ=255)
                                             aaaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCtttttt                                                                      <  2:2456615/65‑1 (MQ=255)
                                             aaaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctct                                            <  2:1969553/91‑1 (MQ=255)
                                               aTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgt                                          <  2:585511/91‑1 (MQ=255)
                                                 tATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtct                                        <  2:2110324/91‑1 (MQ=255)
                                                  aTCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtctc                                       <  2:18238/91‑1 (MQ=255)
                                                  aTCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtctc                                       <  1:614836/91‑1 (MQ=255)
                                                   tCAAATACATTTGAAAAAAGCTGACTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATgtctctgtctct                                      <  2:1515520/91‑1 (MQ=255)
                                                    cAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCtt                                     >  2:3628603/1‑91 (MQ=255)
                                                    cAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCtt                                     <  1:3207087/91‑1 (MQ=255)
                                                     aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                    <  1:573536/91‑1 (MQ=255)
                                                     aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                    >  1:2114545/1‑91 (MQ=255)
                                                     aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                    >  1:1737483/1‑91 (MQ=255)
                                                     aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                    >  1:1009539/1‑91 (MQ=255)
                                                      aaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTg                                   <  2:1195822/91‑1 (MQ=255)
                                                        tACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGa                                 >  2:2613008/1‑91 (MQ=255)
                                                          cATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAAc                               <  2:342243/91‑1 (MQ=255)
                                                          cATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAAc                               <  1:2886030/91‑1 (MQ=255)
                                                           aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                              >  1:2766650/1‑91 (MQ=255)
                                                           aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                              <  2:165138/91‑1 (MQ=255)
                                                           aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAACCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                              >  1:240396/1‑91 (MQ=255)
                                                            tttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTGCATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACtt                             <  2:3192811/91‑1 (MQ=255)
                                                             ttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACttt                            <  1:2048975/91‑1 (MQ=255)
                                                              tGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTta                           >  2:371427/1‑91 (MQ=255)
                                                                aaaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtata                         <  1:254713/91‑1 (MQ=255)
                                                                aaaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtata                         <  1:3588707/91‑1 (MQ=255)
                                                                    aaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAAATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTc                     >  2:2989116/1‑91 (MQ=255)
                                                                     aGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTACATTTCa                    >  2:3876806/1‑91 (MQ=255)
                                                                      gCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCAt                   <  2:726581/91‑1 (MQ=255)
                                                                      gCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCAt                   >  1:390810/1‑91 (MQ=255)
                                                                       cTGATTTCTATCAAATAATTTATATAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATa                  <  2:3618406/91‑1 (MQ=255)
                                                                         gATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAg                >  2:1573353/1‑91 (MQ=255)
                                                                         gATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCAGAAg                >  1:3009726/1‑91 (MQ=255)
                                                                          aTTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt               <  1:36560/91‑1 (MQ=255)
                                                                          aTTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt               >  2:1408046/1‑91 (MQ=25)
                                                                          aTTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt               >  1:1908907/1‑91 (MQ=25)
                                                                          aTTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt               >  2:205231/1‑91 (MQ=25)
                                                                           tttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTa              <  2:3219352/91‑1 (MQ=255)
                                                                            ttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTAc             >  2:1291844/1‑91 (MQ=255)
                                                                             tCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACt            <  1:3319521/91‑1 (MQ=255)
                                                                              cTATGAAGTAAGTTATAGGGATCAGCTTTTTTACATTGTCTAAGAACTTAATGGCTCTGTCTCTTTGGAACTTTATATTTCATGAGTACTa           <  2:1051932/91‑1 (MQ=25)
                                                                               tATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAt          <  2:911465/91‑1 (MQ=255)
                                                                                  cAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATCAGTACTCtttt       >  1:1286645/1‑91 (MQ=255)
                                                                                       |                                                                                           
GGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGT  >  USA300TCH1516_ALE/2133403‑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: