Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 9,710 T→C 14.9% intergenic (+15/+72) gyrA → / ← nnrD DNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE9,7100TC14.9% 533.6 / 76.3 216intergenic (+15/+72)gyrA/nnrDDNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (96/87);  new base C (12/20);  total (109/107)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-01

Only 100 of 221 total aligned reads displayed.

TGGACGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTC  >  USA300TCH1516_ALE/9620‑9800
                                                                                          |                                                                                          
 ggACGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGCAGAATAATAAAAAATAAGACTTc                                                                                           >  2:776396/1‑91 (MQ=255)
   aCGTATTGAAGTAAGACACGATTTCATGGATCGTGTTGAAGAAGATCTAAAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTcac                                                                                         >  1:3084452/1‑89 (MQ=255)
   aCGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTccc                                                                                         >  2:476541/1‑91 (MQ=255)
    cGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCGGATGAAGAATAATAAAAAATAAGACTCCCCt                                                                                        <  1:3045155/91‑1 (MQ=255)
    cGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCt                                                                                        >  2:2812418/1‑91 (MQ=255)
      tATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCtat                                                                                      >  2:1170350/1‑91 (MQ=255)
            aGTAAGTCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAg                                                                                <  2:2743096/91‑1 (MQ=38)
                aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                            >  1:265505/1‑91 (MQ=255)
                aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                            <  1:2980280/91‑1 (MQ=255)
                aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                            >  2:3090906/1‑91 (MQ=255)
                 gACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAg                                                                           >  2:121608/1‑91 (MQ=255)
                 gACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgggaag                                                                           <  2:1025914/91‑4 (MQ=37)
                   cAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTc                                                                         <  1:2637415/91‑1 (MQ=255)
                   cAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTc                                                                         <  1:3813777/91‑1 (MQ=37)
                    aaGATTTCATGGATCGTGTTGAAGACTATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                        <  2:3359971/91‑1 (MQ=255)
                    aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                        >  2:2063255/1‑91 (MQ=255)
                    aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                        <  2:2104534/91‑1 (MQ=255)
                    aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                        <  2:3172923/91‑1 (MQ=37)
                    aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                        <  2:3765986/91‑1 (MQ=37)
                     aGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCtt                                                                       <  1:1747524/91‑1 (MQ=255)
                     aGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCtt                                                                       <  1:1085083/91‑1 (MQ=37)
                      gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTa                                                                      <  1:2743624/91‑1 (MQ=255)
                        tttCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtt                                                                    >  1:473160/1‑91 (MQ=255)
                         ttCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAttt                                                                   <  1:874948/91‑1 (MQ=255)
                          tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCTTAtttt                                                                  <  1:415307/91‑1 (MQ=25)
                           cATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAttttt                                                                 >  2:1883424/1‑91 (MQ=255)
                            aTGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTa                                                                >  1:2327656/1‑91 (MQ=255)
                             tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                               >  1:3233299/1‑91 (MQ=255)
                             tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                               >  2:2468132/1‑91 (MQ=255)
                             tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                               >  2:886419/1‑91 (MQ=255)
                               gATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAGAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGc                                                             >  1:3866287/1‑91 (MQ=255)
                                aTCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCt                                                            >  1:336553/1‑91 (MQ=255)
                                 tCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTa                                                           <  1:429565/91‑1 (MQ=255)
                                   gtgtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGa                                                         <  1:204773/91‑1 (MQ=255)
                                   gtgtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGa                                                         >  2:1774661/1‑91 (MQ=255)
                                   gtgtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGa                                                         <  2:1426554/91‑1 (MQ=255)
                                      ttGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGACTTATTTTTATGCTAGAAAg                                                      >  2:3313990/1‑91 (MQ=255)
                                       tGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGt                                                     <  1:1006314/91‑1 (MQ=255)
                                       tGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGt                                                     <  1:2276083/91‑1 (MQ=255)
                                       tGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGt                                                     <  2:1398590/91‑1 (MQ=255)
                                        gaagaaGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTa                                                    <  2:2600629/91‑1 (MQ=255)
                                        gaagaaGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTa                                                    <  2:2553211/91‑1 (MQ=255)
                                         aagaagATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTaa                                                   <  2:3883957/91‑1 (MQ=255)
                                           gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATg                                                 <  2:109857/91‑1 (MQ=255)
                                              gaTATACAACAATCATGAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCtt                                              <  2:2978508/91‑1 (MQ=255)
                                                 ataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTAc                                           >  2:2051497/1‑91 (MQ=255)
                                                 ataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTAc                                           >  1:254462/1‑91 (MQ=255)
                                                  taCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACt                                          >  1:1446843/1‑91 (MQ=255)
                                                  taCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACt                                          >  1:858944/1‑91 (MQ=255)
                                                  taCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACt                                          >  1:3576586/1‑91 (MQ=255)
                                                     aacaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACtata                                       <  1:1336344/91‑1 (MQ=255)
                                                      acaaTCATCATATGAAGAATAATAAATAATAAGACTTCCCTATGTTGAGGGGGGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTatat                                      <  2:44048/91‑1 (MQ=255)
                                                       caaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATAtt                                     >  2:1349691/1‑91 (MQ=255)
                                                        aaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTc                                    <  2:1700273/91‑1 (MQ=255)
                                                        aaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTc                                    >  1:2753485/1‑91 (MQ=255)
                                                         atcatcAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCa                                   >  1:1456121/1‑91 (MQ=255)
                                                           catcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAAt                                 >  2:571054/1‑91 (MQ=255)
                                                            atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                <  1:1127577/91‑1 (MQ=255)
                                                            atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                >  2:2409963/1‑91 (MQ=255)
                                                            atcaGATGAAGAATAAAAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                >  2:1583886/1‑91 (MQ=255)
                                                             tcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGa                               >  2:752515/1‑91 (MQ=255)
                                                             tcaGATGAAGAATAATAAAAAATAAGACTTCCCTAGATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGa                               >  1:3243079/1‑91 (MQ=255)
                                                              caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                              <  1:571054/91‑1 (MQ=255)
                                                              caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                              >  2:2783743/1‑91 (MQ=255)
                                                                gATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTa                            >  1:3494485/1‑91 (MQ=255)
                                                                aaTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTa                            <  2:1617634/90‑1 (MQ=255)
                                                                   gaagaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTa                         <  2:482854/91‑1 (MQ=255)
                                                                   gaagaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTGTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTa                         <  2:1653154/91‑1 (MQ=255)
                                                                     agaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAAt                       <  1:2061705/91‑1 (MQ=255)
                                                                     agaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAAt                       >  2:1512540/1‑91 (MQ=255)
                                                                     agaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAAt                       >  1:1187103/1‑91 (MQ=255)
                                                                      gaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATg                      >  1:3287858/1‑91 (MQ=255)
                                                                       aataatAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGa                     <  1:1661965/91‑1 (MQ=255)
                                                                       aataatAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGa                     <  2:2659026/91‑1 (MQ=255)
                                                                        ataataAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGAc                    <  1:1263104/91‑1 (MQ=255)
                                                                        ataataAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGAc                    <  1:3552024/91‑1 (MQ=255)
                                                                         taataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACt                   >  2:798831/1‑91 (MQ=255)
                                                                          aataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTa                  <  2:3266897/91‑1 (MQ=255)
                                                                          aataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTa                  >  1:634478/1‑91 (MQ=255)
                                                                           ataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTaa                 >  1:366135/1‑91 (MQ=255)
                                                                            taaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAAc                <  2:1047868/91‑1 (MQ=255)
                                                                               aaaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACttt             <  2:1881626/91‑1 (MQ=255)
                                                                               aaaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACttt             >  1:1845795/1‑91 (MQ=255)
                                                                               aaaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACttt             >  1:1690198/1‑91 (MQ=255)
                                                                                  aTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTa          >  1:1190257/1‑91 (MQ=25)
                                                                                   tAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa         <  2:1981519/91‑1 (MQ=255)
                                                                                    aaGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt        <  1:2002934/91‑1 (MQ=255)
                                                                                    aaGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt        <  1:516995/91‑1 (MQ=255)
                                                                                    aaGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt        <  1:248101/91‑1 (MQ=255)
                                                                                    aaGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt        >  1:3839438/1‑91 (MQ=37)
                                                                                     aGACTTCCCTATATGTAGGGGCGACTTATTTTTATGCTAGAAAGTCATGCCTTACTATATTCAATGATTAGTAATGACTAAATTTCTAAtt       >  2:1809056/1‑91 (MQ=255)
                                                                                     aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  2:3471470/1‑91 (MQ=255)
                                                                                     aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  1:782159/1‑91 (MQ=255)
                                                                                     aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  2:530721/1‑91 (MQ=255)
                                                                                     aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  1:3410263/1‑91 (MQ=255)
                                                                                     aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  1:1519004/1‑91 (MQ=37)
                                                                                     aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  2:2085558/1‑91 (MQ=37)
                                                                                      gACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg      >  1:195310/1‑91 (MQ=37)
                                                                                        cTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGtt    >  2:3704827/1‑91 (MQ=255)
                                                                                          tCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTc  >  2:3110502/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGGACGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTC  >  USA300TCH1516_ALE/9620‑9800

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: