Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,563,763 C→G 8.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,7630CG 8.1% 495.1 / 28.0 148intergenic (+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 (85/51);  new base G (5/7);  total (90/58)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 153 total aligned reads displayed.

AAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTA  >  USA300TCH1516_ALE/2563675‑2563853
                                                                                        |                                                                                          
aaaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCaa                                                                                          >  1:1616358/1‑91 (MQ=255)
aaaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCaa                                                                                          >  1:371331/1‑91 (MQ=255)
  aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                        <  1:2292212/91‑2 (MQ=255)
  aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                        <  1:2309919/91‑2 (MQ=255)
   gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAActg                                                                                       <  1:2834841/91‑3 (MQ=255)
    cTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAActga                                                                                      <  1:1053412/91‑4 (MQ=255)
     tCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGAt                                                                                     >  1:735710/1‑91 (MQ=255)
      cgcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATa                                                                                    >  1:2238348/1‑91 (MQ=255)
       gcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATAc                                                                                   >  1:3392044/1‑91 (MQ=255)
         gTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTa                                                                                 <  1:1550807/91‑1 (MQ=255)
          tATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTag                                                                                <  1:1518766/91‑1 (MQ=255)
          tATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTag                                                                                <  1:1518746/91‑1 (MQ=255)
              aacaacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATAGTAGAGtt                                                                            >  1:2082535/1‑91 (MQ=255)
               acaacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGttt                                                                           <  1:1088670/91‑1 (MQ=255)
                caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                          >  1:3842843/1‑91 (MQ=255)
                caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTACCAGTTCAACCGATACTAGAGtttt                                                                          >  1:1947709/1‑91 (MQ=255)
                caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATCTAAACATCTCGCCCTAGTATCAGTTCAACCGATACTAGAGattt                                                                          >  1:2745043/1‑91 (MQ=255)
                 aacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGttttt                                                                         >  1:1360836/1‑91 (MQ=255)
                    gACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTg                                                                      <  1:1994080/91‑1 (MQ=255)
                    gACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTg                                                                      >  1:956047/1‑91 (MQ=255)
                      cAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaa                                                                    <  1:976067/91‑1 (MQ=255)
                      cAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaa                                                                    >  1:1923625/1‑91 (MQ=255)
                      cAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaa                                                                    >  1:3440414/1‑91 (MQ=255)
                      cAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaa                                                                    >  1:102273/1‑91 (MQ=255)
                            gACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAAc                                                              <  1:1042930/91‑1 (MQ=255)
                             aCAGCGGTCATTTCAGATATGACAATCAATAATTTAAGCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACt                                                             >  1:3502743/1‑91 (MQ=255)
                             aCAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACt                                                             <  1:2516323/91‑1 (MQ=255)
                              cAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACtt                                                            <  1:3457920/91‑1 (MQ=255)
                              cAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACtt                                                            >  1:2545917/1‑91 (MQ=255)
                              cAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTATTTTCAAAAAAACtc                                                            <  1:630480/91‑2 (MQ=255)
                               cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                           >  1:1871853/5‑91 (MQ=255)
                               cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                           >  1:1871854/5‑91 (MQ=255)
                                gCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCa                                                          <  1:4105484/91‑1 (MQ=255)
                                 cGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAt                                                         >  1:2820632/1‑91 (MQ=255)
                                   gTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAttt                                                       >  1:2539068/1‑91 (MQ=255)
                                     cATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTa                                                     >  1:1148647/1‑91 (MQ=255)
                                     cATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTa                                                     <  1:3770832/91‑1 (MQ=255)
                                     cATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTa                                                     <  1:1368074/91‑1 (MQ=255)
                                     cATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTa                                                     <  1:1371671/91‑1 (MQ=255)
                                      aTTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTAt                                                    <  1:2395583/91‑1 (MQ=255)
                                         tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                 >  1:354397/1‑91 (MQ=255)
                                         tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                 >  1:1717/1‑91 (MQ=255)
                                         tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                 >  1:434842/1‑91 (MQ=255)
                                          cAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAat                                                >  1:2852173/1‑91 (MQ=255)
                                           aGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAata                                               <  1:2288264/91‑1 (MQ=255)
                                            gATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAatat                                              >  1:870968/1‑91 (MQ=255)
                                            gATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAatat                                              >  1:1323760/1‑91 (MQ=255)
                                               atGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATc                                           <  1:501966/91‑1 (MQ=255)
                                               atGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATc                                           <  1:2042423/91‑1 (MQ=255)
                                                tGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCt                                          <  1:1238712/91‑1 (MQ=255)
                                                tGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCt                                          >  1:216644/1‑91 (MQ=255)
                                                   caatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATg                                       <  1:2579458/91‑1 (MQ=255)
                                                   caatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATg                                       <  1:1247666/91‑1 (MQ=255)
                                                   caatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATg                                       >  1:2997561/1‑91 (MQ=255)
                                                    aatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGt                                      >  1:88088/1‑91 (MQ=255)
                                                     atcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTc                                     >  1:480189/1‑91 (MQ=255)
                                                      tcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCa                                    >  1:1165317/1‑91 (MQ=255)
                                                        aataatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCAtt                                  >  1:1478202/1‑91 (MQ=255)
                                                          taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                <  1:1723241/91‑1 (MQ=255)
                                                          taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                >  1:2078970/1‑91 (MQ=255)
                                                          taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                >  1:1805379/1‑91 (MQ=255)
                                                          taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                >  1:3089606/1‑91 (MQ=255)
                                                          taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                <  1:2934395/91‑1 (MQ=255)
                                                           aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                               <  1:1807719/91‑1 (MQ=255)
                                                           aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                               <  1:1808524/91‑1 (MQ=255)
                                                                 aaCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaa                         >  1:3069212/1‑91 (MQ=255)
                                                                  aCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaa                        >  1:3283302/1‑91 (MQ=255)
                                                                   cATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaaa                       >  1:3535499/1‑91 (MQ=255)
                                                                   cATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaaa                       >  1:3254736/1‑91 (MQ=255)
                                                                    aTCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAg                      <  1:3995388/91‑1 (MQ=255)
                                                                     tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                     >  1:1939828/1‑91 (MQ=255)
                                                                     tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                     >  1:607799/1‑91 (MQ=255)
                                                                     tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                     <  1:3997263/91‑1 (MQ=255)
                                                                     tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                     <  1:2579045/91‑1 (MQ=255)
                                                                      aaaactctAGTATCGGTTGAACTGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTg                    >  1:3766097/5‑91 (MQ=37)
                                                                       tcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGa                   <  1:2824873/91‑1 (MQ=255)
                                                                        cGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGAt                  >  1:2306429/1‑91 (MQ=255)
                                                                         gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa                 <  1:1825238/91‑1 (MQ=255)
                                                                         gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa                 >  1:363936/1‑91 (MQ=255)
                                                                          ctctAGTATCGGTTGAACTGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATaa                >  1:4079511/1‑91 (MQ=38)
                                                                             tAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat             >  1:4069680/1‑91 (MQ=255)
                                                                             tAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat             >  1:2473241/1‑91 (MQ=255)
                                                                              aGTATCGGTTGAGCTGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAata            >  1:3182006/1‑91 (MQ=255)
                                                                               gTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAatat           >  1:2412320/1‑91 (MQ=255)
                                                                               gTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAatat           >  1:575513/1‑91 (MQ=255)
                                                                               gTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGCGATAAAatat           >  1:575453/1‑91 (MQ=255)
                                                                                  tCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTaa        >  1:1616274/1‑91 (MQ=255)
                                                                                   cAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAAc       >  1:1464111/1‑91 (MQ=255)
                                                                                    aGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACt      >  1:2031382/1‑91 (MQ=255)
                                                                                    aGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACt      >  1:2034992/1‑91 (MQ=255)
                                                                                      ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt    >  1:3216647/1‑91 (MQ=255)
                                                                                      ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt    >  1:3216772/1‑91 (MQ=255)
                                                                                      ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt    <  1:3571853/91‑1 (MQ=255)
                                                                                      ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt    <  1:38983/91‑1 (MQ=255)
                                                                                      ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt    >  1:3948344/1‑91 (MQ=255)
                                                                                       tCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAtt   >  1:644092/1‑91 (MQ=255)
                                                                                       tCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAtt   <  1:3969324/91‑1 (MQ=255)
                                                                                        cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa  <  1:3182518/91‑1 (MQ=255)
                                                                                        cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa  >  1:628701/1‑91 (MQ=255)
                                                                                        cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa  <  1:804351/91‑1 (MQ=255)
                                                                                        |                                                                                          
AAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTA  >  USA300TCH1516_ALE/2563675‑2563853

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: