Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,563,763 C→G 10.1% intergenic (+28/‑110) cycA_2 → / → nhaK_2 D‑serine/D‑alanine/glycine transporter/Sodium, potassium, lithium and rubidium/H(+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,563,7630CG10.1% 789.1 / 50.7 242intergenic (+28/‑110)cycA_2/nhaK_2D‑serine/D‑alanine/glycine transporter/Sodium, potassium, lithium and rubidium/H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (110/107);  new base G (7/18);  total (117/125)
Fisher's exact test for biased strand distribution p-value = 3.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

GAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTA  >  USA300TCH1516_ALE/2563673‑2563853
                                                                                          |                                                                                          
 aaaaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCa                                                                                           <  1:625988/91‑1 (MQ=255)
  aaaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGaa                                                                                          <  1:451013/91‑1 (MQ=255)
  aaaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCaa                                                                                          >  1:2654451/1‑91 (MQ=255)
   aaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAAc                                                                                         >  1:3197970/1‑91 (MQ=255)
    aGCTGGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATAAGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                        <  1:3292792/91‑2 (MQ=255)
    aGCTCGCGTCTGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAAcc                                                                                        <  2:2720218/91‑1 (MQ=255)
    aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                        <  1:1092653/91‑2 (MQ=255)
    aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                        <  1:426087/91‑2 (MQ=255)
    aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                        <  1:1625704/91‑2 (MQ=255)
    aGCTCGCGTATAAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCACTCTAGTATCGGTTGAAct                                                                                        <  2:1023100/91‑2 (MQ=255)
     gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACCATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCg                                                                                       >  1:3518842/1‑91 (MQ=255)
     gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAActg                                                                                       <  1:144372/91‑3 (MQ=255)
     gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCg                                                                                       >  2:1304298/1‑91 (MQ=255)
     gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCg                                                                                       >  2:2658767/1‑91 (MQ=255)
      cTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAActga                                                                                      <  1:3557114/91‑4 (MQ=255)
      cTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGa                                                                                      >  1:2964576/1‑91 (MQ=255)
       tCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGAt                                                                                     <  2:1481165/91‑1 (MQ=255)
        cgcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATa                                                                                    >  1:1375669/1‑91 (MQ=255)
        cgcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATa                                                                                    >  1:2880383/1‑91 (MQ=255)
         gcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATAc                                                                                   >  2:1722134/1‑91 (MQ=255)
            tATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTag                                                                                <  1:3717489/91‑1 (MQ=255)
              tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                              <  1:1293341/91‑1 (MQ=255)
              tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                              <  1:2400033/91‑1 (MQ=255)
              tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                              >  1:3809591/1‑91 (MQ=255)
              tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                              <  1:537381/91‑1 (MQ=255)
               gAACAGCGACAAAATGACAGCGGTCATCTCGGATATGACAATCCATAATGTAGACGTCTCGCGCTAGTATCAGATCAACCGATACTAGAGt                                                                             <  2:377592/91‑1 (MQ=255)
                aacaacGACAAAATGACTGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtt                                                                            <  1:3387618/91‑1 (MQ=255)
                 acaacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGttt                                                                           >  1:2460425/1‑91 (MQ=255)
                  caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                          >  2:2065001/1‑91 (MQ=255)
                  caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                          >  2:298669/1‑91 (MQ=255)
                  caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                          >  2:2766209/1‑91 (MQ=255)
                     cGACAAAATGACAGCGGTCATTTCGGATATGACAATCAATAATTTAGACATCTCGCTCTAGTATCAGTTCAACCGAGACTAGAGttttttt                                                                       <  2:878594/91‑1 (MQ=255)
                     cGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATTATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGttttttt                                                                       <  2:3187311/91‑1 (MQ=255)
                     cGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGttttttt                                                                       <  2:2523106/91‑1 (MQ=255)
                      gACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTg                                                                      <  2:371059/91‑1 (MQ=255)
                       aCAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGa                                                                     >  2:2933092/1‑91 (MQ=255)
                        cAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaa                                                                    >  2:1364205/1‑91 (MQ=255)
                          aaaTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaaaa                                                                  >  1:3592562/1‑91 (MQ=255)
                          aaaTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaaaa                                                                  >  1:1119760/1‑91 (MQ=255)
                          aaaTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaaaa                                                                  >  1:1920166/1‑91 (MQ=255)
                           aaTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAAt                                                                 >  1:1355151/1‑91 (MQ=255)
                            aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaaattt                                                                <  2:51763/91‑4 (MQ=255)
                              gACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAAc                                                              <  1:3602159/91‑1 (MQ=255)
                              gACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAAc                                                              >  1:3180641/1‑91 (MQ=255)
                               aCAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACt                                                             <  1:1135235/91‑1 (MQ=255)
                               aCAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACt                                                             >  1:835996/1‑91 (MQ=255)
                                 cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                           >  1:2595195/5‑91 (MQ=255)
                                 aGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTATAGTCTTTTTGAAAATAACtca                                                           >  2:2625580/1‑89 (MQ=255)
                                  gctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCa                                                          >  1:946063/4‑91 (MQ=255)
                                   cGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAt                                                         <  1:3004769/91‑1 (MQ=255)
                                    ggTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAtt                                                        >  1:584297/1‑91 (MQ=255)
                                      tCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAATATATCTTCAtttt                                                      >  2:2198866/1‑91 (MQ=255)
                                          ttCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCa                                                  >  2:2803389/1‑91 (MQ=255)
                                             aGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAata                                               <  1:1722842/91‑1 (MQ=255)
                                               atatGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAatata                                             >  2:2423820/1‑91 (MQ=255)
                                                tatGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAatatat                                            <  2:22280/91‑1 (MQ=255)
                                                 atGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATc                                           >  1:750121/1‑91 (MQ=255)
                                                      aatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGt                                      <  1:743659/91‑1 (MQ=255)
                                                      aatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGt                                      <  1:65776/91‑1 (MQ=255)
                                                      aatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGt                                      <  1:2267186/91‑1 (MQ=255)
                                                           atgatTTGAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCACTTTATCAATATATCTATGTCATTa                                 <  2:3568270/88‑1 (MQ=255)
                                                            taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                <  1:2013848/91‑1 (MQ=255)
                                                            taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                >  1:1748871/1‑91 (MQ=255)
                                                             aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                               <  1:850273/91‑1 (MQ=255)
                                                             aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                               >  1:275665/1‑91 (MQ=255)
                                                              atTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGc                              >  1:1547333/1‑91 (MQ=255)
                                                                ttAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGATTATATCTTCATTTTATCAATATATCTATGTCATTATGctt                            >  2:694612/1‑89 (MQ=255)
                                                                   aaCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTTTAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaa                         >  2:772068/1‑91 (MQ=255)
                                                                   aaCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaa                         >  2:2486561/1‑91 (MQ=255)
                                                                    aCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaa                        >  1:2857564/1‑91 (MQ=255)
                                                                     cATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaaa                       <  1:3916740/91‑1 (MQ=255)
                                                                      aTCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAg                      <  1:1242294/91‑1 (MQ=255)
                                                                       tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                     <  2:64052/91‑1 (MQ=255)
                                                                        ctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTg                    >  2:601112/1‑91 (MQ=255)
                                                                        ctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTg                    <  1:1526810/91‑1 (MQ=255)
                                                                         tcGCTCTCGTGTCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGa                   <  2:3513712/91‑1 (MQ=255)
                                                                         tcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGa                   <  2:1646571/91‑1 (MQ=255)
                                                                          cGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGAt                  <  1:1737857/91‑1 (MQ=255)
                                                                            ctctAGTATCAGTTCAACCGATACTAGAGTTTTTTTGTAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATaa                >  2:1814722/1‑91 (MQ=255)
                                                                            ctctAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATaa                <  1:3281260/91‑1 (MQ=255)
                                                                            ctctAGTATCAGTTCAACCGATACTAGAGTTTTGTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTGAAGAGTGATaa                <  2:547575/91‑1 (MQ=255)
                                                                               tAGTATCGGTTGAACTGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat             >  1:254690/1‑91 (MQ=39)
                                                                               tAGTATCGGTTGAACTGATACTAGAGTCTTTTTGCAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat             >  2:1535876/1‑91 (MQ=38)
                                                                               tAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat             >  2:2383632/1‑91 (MQ=255)
                                                                               tAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat             >  1:3772824/1‑91 (MQ=255)
                                                                                   aTCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTa         <  1:1433057/91‑1 (MQ=255)
                                                                                   aTCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAATATTaa        >  2:1914185/1‑91 (MQ=255)
                                                                                    tCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTaa        >  1:222230/1‑91 (MQ=255)
                                                                                    tCAGTTCAACCGATACTAGAGTTTTTCTGAAAATAACTTCATTTTATCAATATACCTATGTCATTATGTCTTACAAGTGATAAAACATTaa        >  2:3196052/1‑91 (MQ=255)
                                                                                     cAGTTCAACCGATACTAGAGTTTTTTTGGAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAAc       <  2:3035469/91‑1 (MQ=255)
                                                                                      ggTTGAACTGTTACTAGAGTTTTTTTGAACATATCTTTATTTGATCCATATATCCATGTCATTCGGCCTTAAACTTGCTAAAATATTAACt      >  2:2704634/2‑91 (MQ=255)
                                                                                      aGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACt      <  1:1861986/91‑1 (MQ=255)
                                                                                       gTTGAACTGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTa     >  1:2346989/1‑91 (MQ=255)
                                                                                        ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt    <  2:3101138/91‑1 (MQ=255)
                                                                                        ctCAACCGATACTAGAGTTTTTTTGAAAATAACTTCTTTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt    <  2:3919400/90‑1 (MQ=255)
                                                                                         tCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAtt   <  2:1379464/91‑1 (MQ=255)
                                                                                          cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAACGTGATAAAATATTAACTATTa  >  2:3305117/1‑91 (MQ=255)
                                                                                          cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa  >  1:181042/1‑91 (MQ=255)
                                                                                          cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa  >  2:1867168/1‑91 (MQ=255)
                                                                                          cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa  >  2:2138949/1‑91 (MQ=255)
                                                                                          |                                                                                          
GAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTA  >  USA300TCH1516_ALE/2563673‑2563853

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: