Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,133,485 Δ1 bp 6.3% intergenic (+329/+43) 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,4850C. 6.3% 871.3 / 48.3 252intergenic (+329/+43)dapE/USA300HOU_RS10885putative succinyl‑diaminopimelate desuccinylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (108/124);  new base . (6/10);  total (114/138)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 255 total aligned reads displayed.

TCGTAAAAGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTT  >  USA300TCH1516_ALE/2133395‑2133575
                                                                                          |                                                                                           
 cGTAAAAGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCa                                                                                            >  1:2250466/1‑91 (MQ=255)
 cGTAAAAGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCa                                                                                            >  2:2144863/1‑91 (MQ=255)
 cGTAAAAGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCa                                                                                            >  1:296843/1‑91 (MQ=255)
  ggaaaaGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACAATTGAAAAAAGCTGATTTCTATCaa                                                                                           <  2:3029892/89‑1 (MQ=255)
    aaaaGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAaat                                                                                         <  1:1187830/91‑1 (MQ=255)
    aaaaGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAaat                                                                                         <  1:1290353/91‑1 (MQ=255)
         gAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTAtttg                                                                                   <  2:2609583/91‑2 (MQ=25)
            tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtata                                                                                >  1:2082256/1‑91 (MQ=255)
            tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTAtttgata                                                                                <  2:1681731/91‑5 (MQ=25)
             gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                               >  2:1335475/1‑91 (MQ=255)
             gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                               >  2:2945489/1‑91 (MQ=255)
             gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                               >  1:3666968/1‑91 (MQ=255)
             gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                               >  2:2836411/1‑91 (MQ=255)
              ggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGa                                                                              <  1:2830459/91‑1 (MQ=25)
               gACATAAATCCCTAGAAAACAGCAGCTCAACGGGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTGCTATCAAATAATTTATAGaa                                                                             <  1:173126/91‑1 (MQ=255)
               gACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATAAAAATAAAAAAAGCTGATTTCTATAAATTATTTGATAGaa                                                                             <  2:3423724/91‑1 (MQ=18)
                 cATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAAt                                                                           >  1:2443881/1‑91 (MQ=255)
                  aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATc                                                                          <  1:2704389/91‑1 (MQ=255)
                  aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATc                                                                          <  2:2909121/91‑1 (MQ=255)
                  aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATc                                                                          <  2:3733886/91‑1 (MQ=255)
                   tAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAACTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCa                                                                         <  1:3084453/91‑1 (MQ=255)
                   tAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATCa                                                                         <  1:1072447/91‑1 (MQ=21)
                        cccTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttt                                                                    <  1:1057015/91‑1 (MQ=255)
                        cccTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttt                                                                    >  2:1767039/1‑91 (MQ=255)
                          cTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttttt                                                                  <  2:3737003/91‑1 (MQ=255)
                             gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                               >  1:2582118/1‑91 (MQ=255)
                             gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                               >  1:541878/1‑91 (MQ=255)
                             gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                               >  1:520368/1‑91 (MQ=255)
                             gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                               <  1:439338/91‑1 (MQ=255)
                             gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATCAGCttttttca                                                               <  1:573989/91‑3 (MQ=18)
                              aaaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACa                                                              >  2:2463937/1‑91 (MQ=255)
                              aaaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACa                                                              <  2:69777/91‑1 (MQ=255)
                              aaaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACa                                                              >  2:2465479/1‑91 (MQ=255)
                               aaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACAt                                                             <  1:951531/91‑1 (MQ=255)
                                aaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACAtt                                                            >  2:1265663/1‑91 (MQ=255)
                                   agcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTc                                                         <  2:149563/91‑1 (MQ=255)
                                    gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                        <  1:406795/91‑1 (MQ=255)
                                    gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                        <  2:1051521/91‑1 (MQ=255)
                                      agcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTaa                                                      <  2:170995/91‑1 (MQ=255)
                                         tCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGaa                                                   >  2:766463/1‑91 (MQ=255)
                                         tCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGaa                                                   >  1:715557/1‑91 (MQ=255)
                                          cgttgAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                  >  1:3219352/5‑91 (MQ=255)
                                           aaCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACt                                                 <  2:2413574/91‑1 (MQ=255)
                                           aaCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACt                                                 <  2:2292708/91‑1 (MQ=255)
                                            aCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACtt                                                <  1:3835214/91‑1 (MQ=255)
                                            aCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACtt                                                <  2:3688403/91‑1 (MQ=255)
                                            aCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACtt                                                >  2:2385820/1‑91 (MQ=255)
                                              gAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTaa                                              <  2:323502/91‑1 (MQ=255)
                                              gAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTaa                                              >  2:3131883/1‑91 (MQ=255)
                                                gCTGTGAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtg                                            >  2:900418/1‑91 (MQ=255)
                                                 cTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgt                                           >  2:264954/1‑91 (MQ=255)
                                                  tGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtc                                          >  2:2278207/1‑91 (MQ=255)
                                                    tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                        >  1:423334/1‑91 (MQ=255)
                                                    tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                        >  2:3495456/1‑91 (MQ=255)
                                                    tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                        <  1:1891840/91‑1 (MQ=255)
                                                    tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                        >  1:2107530/1‑91 (MQ=255)
                                                    tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                        >  1:3495161/1‑91 (MQ=255)
                                                    tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                        <  1:629179/91‑1 (MQ=255)
                                                    tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                        >  1:3192811/1‑91 (MQ=255)
                                                     aaaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctct                                       >  2:2633642/1‑91 (MQ=255)
                                                      aaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctg                                      <  2:3618522/91‑1 (MQ=255)
                                                       aTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgt                                     <  2:1547691/91‑1 (MQ=255)
                                                        ttATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtc                                    <  2:520368/91‑1 (MQ=255)
                                                          aTCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtctc                                  <  2:18238/91‑1 (MQ=255)
                                                           tctaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATgtctctgtctct                                 <  2:3202742/88‑1 (MQ=255)
                                                           tCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATgtctctgtctct                                 <  2:3528000/91‑1 (MQ=255)
                                                             aaaTACATTTGAAAAAAGCTGCTTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                               >  2:572847/1‑91 (MQ=255)
                                                             aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                               >  1:1009539/1‑91 (MQ=255)
                                                             aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                               >  1:2114545/1‑91 (MQ=255)
                                                              aaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTg                              <  2:1195822/91‑1 (MQ=255)
                                                               agaCATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTgg                             <  1:882012/89‑1 (MQ=255)
                                                                  cATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAAc                          <  1:37484/91‑1 (MQ=255)
                                                                  cATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAAc                          <  1:2886030/91‑1 (MQ=255)
                                                                   aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                         >  1:240474/1‑91 (MQ=255)
                                                                   aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                         >  1:2766650/1‑91 (MQ=255)
                                                                   aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                         <  2:165138/91‑1 (MQ=255)
                                                                   aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                         <  2:1942712/91‑1 (MQ=255)
                                                                    tttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACtt                        <  1:1569996/91‑1 (MQ=255)
                                                                     ttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACttt                       <  1:185453/91‑1 (MQ=255)
                                                                     ttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACttt                       <  1:2048975/91‑1 (MQ=255)
                                                                       gAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtat                     <  2:75187/91‑1 (MQ=255)
                                                                        aaaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtata                    <  1:254664/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)
                                                                               cTGATTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTACGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATa             >  2:3408098/1‑91 (MQ=21)
                                                                                 gATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAg           >  2:1573353/1‑91 (MQ=255)
                                                                                  aTTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt          <  1:36560/91‑1 (MQ=255)
                                                                                  aTTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTTGAACTTTATATTTCATAAGt          >  2:2770791/1‑91 (MQ=21)
                                                                                  aTTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt          >  2:205267/1‑91 (MQ=25)
                                                                                   tttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTa         >  1:3234242/1‑91 (MQ=255)
                                                                                    ttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTAc        >  2:1986099/1‑91 (MQ=255)
                                                                                     tCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACt       <  1:959795/91‑1 (MQ=255)
                                                                                     tCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACt       <  1:3319521/91‑1 (MQ=255)
                                                                                     tCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCACAAGTACt       >  2:3459578/1‑91 (MQ=255)
                                                                                      cTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTa      <  2:602709/91‑1 (MQ=255)
                                                                                       tATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAt     <  1:238551/91‑1 (MQ=255)
                                                                                       tATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAt     <  2:3793273/91‑1 (MQ=255)
                                                                                        aTCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAtt    >  1:3039062/1‑91 (MQ=255)
                                                                                          cAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATCAGTACTCtttt  >  1:1286645/1‑91 (MQ=255)
                                                                                          |                                                                                           
TCGTAAAAGGAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTT  >  USA300TCH1516_ALE/2133395‑2133575

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: