Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,133,494:1 +G 5.6% intergenic (+338/+34) 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,4941.G 5.6% 1017.3 / 52.4 232intergenic (+338/+34)dapE/USA300HOU_RS10885putative succinyl‑diaminopimelate desuccinylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (95/124);  new base G (6/7);  total (101/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 251 total aligned reads displayed.

GAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGTTTA  >  USA300TCH1516_ALE/2133404‑2133583
                                                                                           |                                                                                         
gAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTAtttg                                                                                           <  2:2609583/91‑2 (MQ=25)
  aTGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtat                                                                                         >  1:3178463/1‑91 (MQ=255)
   tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtata                                                                                        <  1:470024/91‑1 (MQ=255)
   tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTtata                                                                                        >  2:2315818/1‑91 (MQ=255)
   tGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTAtttgata                                                                                        <  2:1681813/91‑5 (MQ=25)
    gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAGATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                       >  1:3357201/1‑91 (MQ=255)
    gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                       <  1:1583288/91‑1 (MQ=255)
    gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                       >  1:3666968/1‑91 (MQ=255)
    gggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAg                                                                                       <  2:3429259/91‑1 (MQ=255)
     ggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGa                                                                                      >  1:3078014/1‑91 (MQ=255)
     ggACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATAAATTATTTGATAGa                                                                                      <  1:2830459/91‑1 (MQ=25)
         aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATc                                                                                  <  2:3733886/91‑1 (MQ=255)
         aTAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATc                                                                                  <  2:2909179/91‑1 (MQ=255)
          tAAATCCCTAGAAAACAGCAGCTCAACGCGCTGTAAATTATCAAATACATGTGAAAAAAGCTGATTTCTATAAATTATTTGATAGAAATCa                                                                                 <  2:759044/91‑1 (MQ=18)
          tAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCa                                                                                 <  1:2075154/91‑1 (MQ=255)
             aTCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCt                                                                              <  2:3198352/91‑1 (MQ=255)
               cccTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttt                                                                            >  2:2245210/1‑91 (MQ=255)
               cccTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttt                                                                            >  2:1767435/1‑91 (MQ=255)
                ccTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCtttt                                                                           >  1:2330525/1‑91 (MQ=255)
                 cTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttttt                                                                          <  2:1033846/91‑1 (MQ=255)
                 cTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCttttt                                                                          <  1:943195/91‑1 (MQ=255)
                    gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                       >  1:541878/1‑91 (MQ=255)
                    gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                       <  1:439338/91‑1 (MQ=255)
                    gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                       >  1:3220583/1‑91 (MQ=255)
                    gAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTAc                                                                       >  1:2169807/1‑91 (MQ=255)
                     aaaaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACa                                                                      >  2:2463937/1‑91 (MQ=255)
                       aaCAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACAtt                                                                    >  2:1265663/1‑91 (MQ=255)
                         cagcagCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGt                                                                  <  2:3516526/91‑1 (MQ=255)
                           gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                                <  1:406795/91‑1 (MQ=255)
                           gcagcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCt                                                                >  1:186584/1‑91 (MQ=255)
                             agcTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTaa                                                              <  2:3734502/91‑1 (MQ=255)
                                 cgttgAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                          >  1:1969553/5‑91 (MQ=255)
                                  gttgAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACt                                                         >  1:2927498/4‑91 (MQ=255)
                                  aaCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACt                                                         <  1:2102329/91‑1 (MQ=255)
                                  aaCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACt                                                         <  2:2292708/91‑1 (MQ=255)
                                  aaCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACt                                                         <  2:2413574/91‑1 (MQ=255)
                                   aCGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACtt                                                        <  1:3835214/91‑1 (MQ=255)
                                    cGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                          <  2:1749743/88‑1 (MQ=255)
                                    cGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAAc                                                          >  1:1749743/1‑88 (MQ=255)
                                     gAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTaa                                                      <  2:323502/91‑1 (MQ=255)
                                     gAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTaa                                                      >  2:3131883/1‑91 (MQ=255)
                                        cTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgt                                                   >  2:264954/1‑91 (MQ=255)
                                        cTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgt                                                   >  2:686339/1‑91 (MQ=255)
                                         tGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtc                                                  <  1:1057984/91‑1 (MQ=255)
                                           tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                                >  2:573536/1‑91 (MQ=255)
                                           tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                                <  1:629179/91‑1 (MQ=255)
                                           tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                                >  2:3495456/1‑91 (MQ=255)
                                           tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                                >  1:2109988/1‑91 (MQ=255)
                                           tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                                >  2:2802889/1‑91 (MQ=255)
                                           tAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctc                                                >  1:423334/1‑91 (MQ=255)
                                            aaaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctct                                               <  2:1969553/91‑1 (MQ=255)
                                            aaaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctct                                               >  2:2633642/1‑91 (MQ=255)
                                             aaTTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctg                                              <  2:3618522/91‑1 (MQ=255)
                                               ttATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtc                                            <  2:520368/91‑1 (MQ=255)
                                                tATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtct                                           <  2:2110324/91‑1 (MQ=255)
                                                tATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtct                                           >  2:853785/1‑91 (MQ=255)
                                                 aTCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAAtgtctctgtctc                                          <  2:18238/91‑1 (MQ=255)
                                                  tCAAATACATTTGAAAAAAGCTGACTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATgtctctgtctct                                         <  2:1515520/91‑1 (MQ=255)
                                                   cAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCtt                                        <  1:3293856/91‑1 (MQ=255)
                                                   cAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCtt                                        >  2:3628603/1‑91 (MQ=255)
                                                    aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                       <  1:573536/91‑1 (MQ=255)
                                                    aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                       >  1:1009539/1‑91 (MQ=255)
                                                    aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                       >  1:1541994/1‑91 (MQ=255)
                                                    aaaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCttt                                       >  1:3207414/1‑91 (MQ=255)
                                                     aaTACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTg                                      <  2:1195822/91‑1 (MQ=255)
                                                      agaCATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTgg                                     <  1:882012/89‑1 (MQ=255)
                                                         cATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAAc                                  <  2:342243/91‑1 (MQ=255)
                                                         cATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAAc                                  <  1:37484/91‑1 (MQ=255)
                                                          aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                                 >  1:2769213/1‑91 (MQ=255)
                                                          aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                                 <  2:165138/91‑1 (MQ=255)
                                                          aTTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACt                                 >  1:240474/1‑91 (MQ=255)
                                                           tttGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACtt                                <  1:2689938/91‑1 (MQ=255)
                                                              gAAAAAAGCTGATTTCGATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGCACTTAATGTCTCTGTCTCTTTGGAACTTtat                             <  2:2970464/91‑1 (MQ=255)
                                                               aaaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtata                            >  1:1283234/1‑91 (MQ=255)
                                                               aaaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtata                            <  1:3588707/91‑1 (MQ=255)
                                                               aaaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTtata                            <  1:254664/91‑1 (MQ=255)
                                                                aaaaaGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTatat                           >  1:24424/1‑91 (MQ=255)
                                                                    aGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTACATTTCa                       >  2:3876806/1‑91 (MQ=255)
                                                                        gATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAg                   >  2:1573353/1‑91 (MQ=255)
                                                                        gATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCAGAAg                   >  1:3009726/1‑91 (MQ=255)
                                                                         aTTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt                  <  1:36560/91‑1 (MQ=255)
                                                                         aTTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTTGAACTTTATATTTCATAAGt                  >  2:2770791/1‑91 (MQ=21)
                                                                         aTTTCTATAAATTATTTGATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGt                  >  2:1408046/1‑91 (MQ=25)
                                                                           ttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTAc                >  1:685357/1‑91 (MQ=255)
                                                                           ttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTAc                <  1:2573783/91‑1 (MQ=255)
                                                                           ttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTAc                >  2:3010623/1‑91 (MQ=255)
                                                                           ttCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTAc                >  2:1986155/1‑91 (MQ=255)
                                                                            tCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACt               <  1:959795/91‑1 (MQ=255)
                                                                             cTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTa              >  2:2874945/1‑91 (MQ=255)
                                                                              tATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAt             <  2:911465/91‑1 (MQ=255)
                                                                              tATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAt             <  2:3793273/91‑1 (MQ=255)
                                                                              tATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAt             <  2:2510077/91‑1 (MQ=255)
                                                                               aTCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTAtt            <  2:712461/91‑1 (MQ=255)
                                                                                 cAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATCAGTACTCtttt          >  1:1286645/1‑91 (MQ=255)
                                                                                     taatTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTg      <  2:186584/91‑1 (MQ=255)
                                                                                       atTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGtt    >  1:3520084/1‑91 (MQ=255)
                                                                                       atTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGtt    >  2:1099911/1‑91 (MQ=255)
                                                                                       atTTATAGAAATCAGCTTTTTTACATCGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGtt    <  2:3357201/91‑1 (MQ=255)
                                                                                         ttATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGTTTa  >  1:2722227/1‑91 (MQ=255)
                                                                                         ttATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGTTTa  <  2:1294751/91‑1 (MQ=255)
                                                                                           |                                                                                         
GAATGGGACATAAATCCCTAGAAAACAGCAGCTCAACGAGCTGTAAATTATCAAATACATTTGAAAAAAGCTGATTTCTATCAAATAATTTATAGAAATCAGCTTTTTTACATTGTCTAAGAACTTAATGTCTCTGTCTCTTTGGAACTTTATATTTCATAAGTACTATTTTATTGTTTA  >  USA300TCH1516_ALE/2133404‑2133583

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: