Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,435,075 T→A 5.8% G144G (GGT→GGA) ‡ dapH → 2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,435,0750TA 5.8% 920.8 / 2.1 277G144G (GGT→GGA) ‡dapH2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (80/175);  new base A (20/2);  total (100/177)
Fisher's exact test for biased strand distribution p-value = 4.25e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

GATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTA  >  USA300TCH1516_ALE/1434985‑1435165
                                                                                          |                                                                                          
gatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                            <  1:4094081/91‑1 (MQ=255)
gatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                            >  1:4356012/1‑91 (MQ=255)
 atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                           <  1:5399206/91‑1 (MQ=255)
  tgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaa                                                                                          >  2:3557344/1‑91 (MQ=255)
   gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                                         <  2:4992812/91‑1 (MQ=255)
   gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTAGTaaa                                                                                         <  1:6643/91‑1 (MQ=255)
     gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                                       <  2:352561/91‑1 (MQ=255)
       gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACAACTGGTAAAAATg                                                                                     >  2:2334034/1‑91 (MQ=255)
       gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                     <  2:3513246/91‑1 (MQ=255)
       gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                     <  2:3751188/91‑1 (MQ=255)
        caacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGt                                                                                    <  2:5255673/91‑1 (MQ=255)
         aacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACAACTGGTAAAAATGTa                                                                                   >  2:6076422/1‑91 (MQ=255)
          acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                                  >  1:2770359/1‑91 (MQ=255)
          acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                                  <  2:2874795/91‑1 (MQ=255)
              ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                              <  1:951155/91‑1 (MQ=255)
                aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAg                                                                            >  1:6683438/1‑91 (MQ=255)
                aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAg                                                                            >  1:1054134/1‑91 (MQ=255)
                  tatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                          <  1:4391024/91‑1 (MQ=255)
                   atTGGGGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgggg                                                                         <  2:6356706/91‑1 (MQ=255)
                       gcgcAGTCGTTGGCGAAGGTACGATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                     <  2:1738697/91‑1 (MQ=255)
                       gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                     <  2:2731824/91‑1 (MQ=255)
                       gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                     <  1:6169678/91‑1 (MQ=255)
                       gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                     >  1:3361744/1‑91 (MQ=255)
                        cgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCGACAACTGGTAAAAATGTACATGTAGGGGCTGgc                                                                    <  1:6946412/91‑1 (MQ=255)
                        cgcAGGCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgc                                                                    <  2:5613095/91‑1 (MQ=255)
                          cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                                  <  1:1954772/91‑1 (MQ=255)
                          cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                                  <  1:1269231/91‑1 (MQ=255)
                           ggacGTTGGCGAAGGTACAATGCTTGGTATGAATGGTACTCTCGGTGGTCGTGATAAAAATGATAAAAATGTACATGTAGGGGCTGGCGCa                                                                 <  2:6880218/88‑1 (MQ=255)
                             tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGt                                                               <  1:3467761/91‑1 (MQ=255)
                              cGTTGGCGAAGGTACAATGATTGATATGAATGCTGCTCTCGGGGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                              <  2:4348401/91‑1 (MQ=255)
                                ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                            >  1:6874712/1‑91 (MQ=255)
                                ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                            <  1:795635/91‑1 (MQ=255)
                                ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTATCATTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTAtt                                                            >  2:3064955/1‑91 (MQ=255)
                                ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCGGGCGCAGTAtt                                                            >  2:6550136/1‑91 (MQ=255)
                                  ggCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                          >  1:3038798/1‑91 (MQ=255)
                                  ggCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTAGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                          >  2:1269231/1‑91 (MQ=255)
                                    ggAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                        <  2:1892199/90‑1 (MQ=255)
                                    cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                        >  2:4181551/1‑91 (MQ=255)
                                    cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                        <  1:3867312/91‑1 (MQ=255)
                                     gAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAg                                                       <  1:3638796/91‑1 (MQ=255)
                                        ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                                    >  1:1628670/1‑91 (MQ=255)
                                         ttACAATAATTAATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGCGCTGGCGCAGTATTAGCAGgtgt                                                   <  2:3238530/90‑1 (MQ=255)
                                         gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                                   <  2:1445674/91‑1 (MQ=255)
                                         gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                                   <  1:1975505/91‑1 (MQ=255)
                                           aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                 <  1:4612904/91‑1 (MQ=255)
                                           aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                 <  1:6215147/91‑1 (MQ=255)
                                            cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                                <  2:5942369/91‑1 (MQ=255)
                                             aaTGATTGATATGAATGCTACTCTCGGTGGTCGCGGTATCATTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAtt                                               >  2:2511970/1‑91 (MQ=255)
                                              aTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                              <  1:4321311/91‑1 (MQ=255)
                                               tgattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                             >  1:5156832/1‑91 (MQ=255)
                                                 attgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAc                                           >  1:2439599/1‑91 (MQ=255)
                                                 attgatATGAATGCTACTCTCGGTGGTCGGGATACAACTAGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAc                                           <  2:2060470/91‑1 (MQ=255)
                                                  ttggtATGAATGCTACTCTCGGTGGTCGTGCGACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                          <  1:6504304/91‑1 (MQ=255)
                                                  ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                          <  1:1000535/91‑1 (MQ=255)
                                                  ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                          <  1:2663937/91‑1 (MQ=255)
                                                  ttgatATGAATGCTACTCTCGGTGGTCGTCGTACCATTGAAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                          >  2:77032/1‑91 (MQ=255)
                                                    gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                        <  2:1825719/91‑1 (MQ=255)
                                                      tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                      >  2:3891519/1‑91 (MQ=255)
                                                      tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                      >  1:5248279/1‑91 (MQ=255)
                                                      tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                      <  1:378501/91‑1 (MQ=255)
                                                      tatGAATGCTACTCTCGGTGGTCGCGGTACAATTGGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                      >  2:5321251/1‑91 (MQ=255)
                                                      tatGAATGCTACTCTCGGTGGTCGCCGTACCACTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                      >  2:2223208/1‑91 (MQ=255)
                                                        tGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACGTGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAg                                    <  1:5263017/91‑1 (MQ=255)
                                                        tGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAg                                    >  1:2914897/1‑91 (MQ=255)
                                                        tGAATGCTACTCTCGGTGGTCGCCGTACCATTGGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAg                                    >  2:265063/1‑91 (MQ=255)
                                                         gAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                                   <  2:2996322/91‑1 (MQ=255)
                                                          aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                  <  2:1149717/91‑1 (MQ=255)
                                                          aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                  <  2:2770359/91‑1 (MQ=255)
                                                            tGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCt                                <  1:37367/91‑1 (MQ=255)
                                                             gCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCtt                               <  1:2178832/91‑1 (MQ=255)
                                                               tACTCTCGGTGGTCGTCGTACCATTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCa                             >  2:3367989/1‑91 (MQ=255)
                                                                aCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAc                            <  1:5735736/91‑1 (MQ=255)
                                                                  tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                          >  2:1497362/1‑91 (MQ=255)
                                                                   ctaggtgggCGTGCTACAACAGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCgg                         <  2:679050/88‑1 (MQ=255)
                                                                     cGGTGGTCGCCCTATCACTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGtt                       >  2:4668186/1‑91 (MQ=255)
                                                                       gtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtat                     <  1:4938104/91‑1 (MQ=255)
                                                                        tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtata                    >  1:2812300/1‑91 (MQ=255)
                                                                        tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtata                    <  1:1451265/91‑1 (MQ=255)
                                                                          gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAAt                  <  1:617983/91‑1 (MQ=255)
                                                                           tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc                 <  2:2439599/91‑1 (MQ=255)
                                                                              tGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAg              <  1:5486193/91‑1 (MQ=255)
                                                                               gCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAgg             <  1:2334034/91‑1 (MQ=255)
                                                                               gCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGAAGGTGTGATTGAACTCCCTAGTGCTTCACCGGTTATAATCGAgg             <  2:5150037/91‑1 (MQ=255)
                                                                                 tACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgat           >  1:72229/1‑91 (MQ=255)
                                                                                   cAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatga         <  1:3557344/91‑1 (MQ=255)
                                                                                   cAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatga         <  2:3257841/91‑1 (MQ=255)
                                                                                   cAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatga         <  2:1628670/91‑1 (MQ=255)
                                                                                     agtGGTAAGAACGTACATGTAGGGGCGAGCGCAGTATCAGCAGGTGTGACTGAATCACGTAGTGCGTCACCGGTTATAATCGAGGatgatg       <  2:6959768/89‑1 (MQ=255)
                                                                                     aCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGatgatg       <  2:5459162/91‑1 (MQ=255)
                                                                                       tGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATATTCGAGGATGATGTa     >  1:5805165/1‑91 (MQ=255)
                                                                                       tGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGatgatg       <  2:4435673/89‑1 (MQ=255)
                                                                                       tGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTa     <  1:4477818/91‑1 (MQ=255)
                                                                                       tGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTa     <  2:1210308/91‑1 (MQ=255)
                                                                                        ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAt    <  2:4721152/91‑1 (MQ=255)
                                                                                        ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAt    <  1:3135487/91‑1 (MQ=255)
                                                                                         gTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAtt   >  2:5758384/1‑91 (MQ=255)
                                                                                         gTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAtt   <  1:185078/91‑1 (MQ=255)
                                                                                         gTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAtt   <  1:3410194/91‑1 (MQ=255)
                                                                                         gTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAtt   >  2:4038009/1‑91 (MQ=255)
                                                                                         gTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAtt   <  2:857598/91‑1 (MQ=255)
                                                                                          |                                                                                          
GATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTA  >  USA300TCH1516_ALE/1434985‑1435165

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