Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,530,786 A→G 10.6% intergenic (‑37/+236) yefM ← / ← bdbD Antitoxin YefM/Disulfide bond formation protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,530,7860AG10.6% 644.8 / 49.8 232intergenic (‑37/+236)yefM/bdbDAntitoxin YefM/Disulfide bond formation protein D
Reads supporting (aligned to +/- strand):  ref base A (102/103);  new base G (3/22);  total (107/125)
Fisher's exact test for biased strand distribution p-value = 2.35e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 246 total aligned reads displayed.

TGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTC  >  USA300TCH1516_ALE/2530696‑2530876
                                                                                          |                                                                                          
 gTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACaaa                                                                                           >  2:645756/1‑91 (MQ=255)
  tCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAat                                                                                          >  1:1063881/1‑91 (MQ=255)
   cATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAata                                                                                         <  1:2256391/91‑1 (MQ=255)
    aTAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatat                                                                                        <  2:2667672/91‑1 (MQ=255)
    aTAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatat                                                                                        <  1:3718666/91‑1 (MQ=255)
      aaGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatattt                                                                                      <  1:826977/91‑3 (MQ=255)
      aaGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatattt                                                                                      <  1:486486/91‑3 (MQ=255)
       aGTTCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAATTATCTg                                                                                     <  1:3070397/91‑1 (MQ=255)
       aGTGCTTTTAAATTCTGTCGCTCGTTTGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTCCAAAGATCTg                                                                                     >  2:3643415/1‑91 (MQ=255)
       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatatttg                                                                                     <  1:1366346/91‑4 (MQ=255)
       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTg                                                                                     >  2:3840901/1‑91 (MQ=255)
       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTg                                                                                     >  2:2411990/1‑91 (MQ=255)
       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTg                                                                                     >  2:461800/1‑91 (MQ=255)
       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAAGATCTg                                                                                     >  2:3644266/1‑91 (MQ=255)
        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                    <  2:893024/91‑1 (MQ=255)
        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                    <  2:1087995/91‑1 (MQ=255)
        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                    <  1:3877546/91‑1 (MQ=255)
        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGt                                                                                    <  1:1395844/91‑1 (MQ=255)
         tGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTa                                                                                   <  1:3293825/91‑1 (MQ=255)
         tGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTa                                                                                   <  2:2716601/91‑1 (MQ=255)
          gCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTAc                                                                                  >  1:2797406/1‑91 (MQ=255)
          gCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTAc                                                                                  >  2:3334505/1‑91 (MQ=255)
           ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACa                                                                                 >  1:3768824/1‑91 (MQ=255)
             tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                               >  2:1850286/1‑91 (MQ=255)
              ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGTACAAtt                                                                              <  1:127452/91‑2 (MQ=255)
              ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATa                                                                              >  1:1450556/1‑91 (MQ=255)
               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             >  2:2322468/1‑91 (MQ=255)
               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             <  2:2583349/91‑1 (MQ=255)
                aaaTTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaaa                                                                            >  2:1794155/1‑91 (MQ=255)
                 aaTTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaaaa                                                                           <  1:2665400/91‑1 (MQ=255)
                     cTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAg                                                                       >  1:1185340/1‑91 (MQ=255)
                     cTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAg                                                                       >  1:2314998/1‑91 (MQ=255)
                      tGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGc                                                                      >  2:3919397/1‑91 (MQ=255)
                       gTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCg                                                                     >  1:3106565/1‑91 (MQ=255)
                        tCGAGCGTAAGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                    >  2:2561251/1‑91 (MQ=255)
                           aGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGacac                                                                 <  1:3529262/91‑1 (MQ=255)
                            gCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGAcaca                                                                <  1:432983/91‑1 (MQ=255)
                              gTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGa                                                              <  1:2368334/91‑1 (MQ=255)
                              gTATGAATAATTTTCAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAAAATCTGCACAATAAAACAAGTGACACAGa                                                              >  2:3355653/1‑91 (MQ=255)
                                 tGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaa                                                           >  1:544037/1‑91 (MQ=255)
                                  gAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaa                                                          <  1:1017877/91‑1 (MQ=255)
                                  gAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaa                                                          <  1:2070746/91‑1 (MQ=255)
                                    ataatTTTTAATAATCATATTTATCGCCTCCTTTTAACACAGCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaaaa                                                        >  1:3430132/1‑91 (MQ=255)
                                    ataatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaaaa                                                        >  1:225334/1‑91 (MQ=255)
                                     taatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAac                                                       >  1:899745/1‑91 (MQ=255)
                                     taatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAac                                                       >  1:3059534/1‑91 (MQ=255)
                                      aatTTTTAATAGTCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAact                                                      <  2:2669555/91‑1 (MQ=255)
                                      aatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAact                                                      <  2:1897999/91‑1 (MQ=255)
                                      aatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCGGTACAATAAAACAAGCGACACAGAAAAAAact                                                      <  2:1195429/91‑1 (MQ=255)
                                        tttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAactac                                                    <  2:1764470/91‑1 (MQ=255)
                                        tttatAATAAACATATTTATCGCGTCCTTTAATCACAACTCAAGTGTACAGATATCTGGACAATAAAACAAGCGACACAGAAAAAAactac                                                    <  2:2239666/91‑1 (MQ=255)
                                            taataaTCATATTTATCGCCGCCTTTTAACACACCTCAATTGTACAAAGATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAAt                                                <  2:1565714/91‑1 (MQ=255)
                                             aataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATc                                               <  1:3543244/91‑1 (MQ=255)
                                             aataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATc                                               >  2:3487181/1‑91 (MQ=255)
                                             aataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATc                                               <  1:2966401/91‑1 (MQ=255)
                                                 atCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGa                                           >  1:2393167/1‑91 (MQ=255)
                                                 atCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATAAAGa                                           >  1:934595/1‑91 (MQ=255)
                                                  tCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAt                                          <  1:2518837/91‑1 (MQ=255)
                                                   cATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAtt                                         <  2:1671809/91‑1 (MQ=255)
                                                      atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                      <  1:995847/91‑1 (MQ=255)
                                                      atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                      <  2:442437/91‑1 (MQ=255)
                                                      atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                      <  2:242866/91‑1 (MQ=255)
                                                          aTCGCCTCCTTTTAACACACCGCAATTGTACAAATATCTGTACAAGAAAACAAGCGACACAGAAAAAAACTGCTAATCAAGATTTCGCACg                                  <  2:3460444/91‑1 (MQ=255)
                                                           tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgg                                 <  2:3311073/91‑1 (MQ=255)
                                                           tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgg                                 <  2:658133/91‑1 (MQ=255)
                                                           tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgg                                 >  2:3860348/1‑91 (MQ=255)
                                                           tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgg                                 >  1:2068999/1‑91 (MQ=255)
                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAATTACTAATCAAGATTTCGCACgggg                               >  2:3688706/1‑91 (MQ=255)
                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               <  1:2459118/91‑1 (MQ=255)
                                                                tcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                            <  2:3065731/91‑1 (MQ=255)
                                                                tcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                            <  1:3597757/91‑1 (MQ=255)
                                                                 cctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaaga                           <  1:2696022/91‑1 (MQ=255)
                                                                   ttttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAAt                         >  1:598849/1‑91 (MQ=255)
                                                                     ttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCa                       <  1:2885693/91‑1 (MQ=255)
                                                                     ttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCa                       <  2:2040600/91‑1 (MQ=255)
                                                                      tAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCaa                      >  1:2987352/1‑91 (MQ=255)
                                                                         cacacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAc                   >  2:1359431/1‑91 (MQ=255)
                                                                         cacacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAc                   <  2:2118291/91‑1 (MQ=255)
                                                                          acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                  >  1:2587491/1‑91 (MQ=255)
                                                                          acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                  <  1:933321/91‑1 (MQ=255)
                                                                          acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                  <  1:303915/91‑1 (MQ=255)
                                                                           cacCTCAATTGTACAAAGATCTGTACGATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAGTTCGCACGGGGAAGAATCAATACCt                 <  1:1679622/91‑1 (MQ=255)
                                                                             ccTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTAt               >  2:2942133/1‑91 (MQ=255)
                                                                             ccTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTAt               >  1:1805618/1‑91 (MQ=255)
                                                                              cTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATc              <  1:3586486/91‑1 (MQ=255)
                                                                              cTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATc              >  1:2853541/1‑91 (MQ=255)
                                                                                 aaTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTg           <  1:2937238/91‑1 (MQ=255)
                                                                                  aTTGTACAAATCTCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATCTCGCACGGGGAAGAATCAATACCTATCCTGt          <  2:3024394/91‑1 (MQ=255)
                                                                                   ttGTACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa         >  1:388447/1‑91 (MQ=255)
                                                                                    tGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc        >  2:2903447/1‑91 (MQ=255)
                                                                                    tGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc        <  2:2445279/91‑1 (MQ=255)
                                                                                     gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTGCTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt       <  1:3141634/91‑1 (MQ=255)
                                                                                     gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt       <  2:1050812/91‑1 (MQ=255)
                                                                                     gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt       <  2:411354/91‑1 (MQ=255)
                                                                                      tACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTg      <  2:3717026/91‑1 (MQ=255)
                                                                                       aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa     >  1:2623189/1‑91 (MQ=255)
                                                                                       aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa     >  2:3235651/1‑91 (MQ=255)
                                                                                       aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa     >  1:3625915/1‑91 (MQ=255)
                                                                                        cAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGAt    >  2:3647465/1‑91 (MQ=255)
                                                                                          aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc  >  2:854296/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTC  >  USA300TCH1516_ALE/2530696‑2530876

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: ‑