Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,563,767 C→T 5.5% intergenic (+32/‑106) 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,7670CT 5.5% 727.8 / 16.9 219intergenic (+32/‑106)cycA_2/nhaK_2D‑serine/D‑alanine/glycine transporter/Sodium, potassium, lithium and rubidium/H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (115/91);  new base T (3/9);  total (119/100)
Fisher's exact test for biased strand distribution p-value = 6.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

AGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTAAAAT  >  USA300TCH1516_ALE/2563677‑2563857
                                                                                          |                                                                                          
   tCGCGTATGAACGACGACAAAATGACAGCGGTCATTTCATATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAACTGAt                                                                                         <  2:2543513/91‑1 (MQ=255)
   tCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAACTGAt                                                                                         <  1:996770/91‑1 (MQ=255)
    cgcgTATGAACAACGACAAAATGACAGCGGTCGTTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATa                                                                                        >  2:3801712/1‑91 (MQ=255)
    cgcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAACTGATa                                                                                        <  2:3123621/91‑1 (MQ=255)
    cgcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATa                                                                                        >  2:341708/1‑91 (MQ=255)
     gcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAACTGATAc                                                                                       <  2:2275380/91‑1 (MQ=39)
     gcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATAc                                                                                       <  1:347737/91‑1 (MQ=255)
      cgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACt                                                                                      <  2:1454554/91‑1 (MQ=255)
       gTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTa                                                                                     <  2:786115/91‑1 (MQ=255)
         aTGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTaga                                                                                   <  1:200926/91‑1 (MQ=255)
          tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                                  >  1:1476995/1‑91 (MQ=255)
          tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                                  <  1:3295056/91‑1 (MQ=255)
           gaagaaCGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGt                                                                                 <  1:3498401/91‑1 (MQ=255)
            aacaacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtt                                                                                >  2:2577886/1‑91 (MQ=255)
              caacGACAAAGTGACAGCGGTCATTTCAGATATGACGATCAATAATTTAAACAGCTCGCTCTAGTATCGGTTGAACTGATACTagagcgag                                                                              <  2:421145/91‑5 (MQ=21)
              caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                              <  2:2324695/91‑1 (MQ=255)
                acGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttttt                                                                            <  1:346553/91‑1 (MQ=255)
                   aCAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGa                                                                         <  1:2528933/91‑1 (MQ=255)
                   aCAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGa                                                                         >  2:1552214/1‑91 (MQ=255)
                   aCAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGa                                                                         <  1:196719/91‑1 (MQ=255)
                        aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATa                                                                    >  1:2413025/1‑91 (MQ=255)
                        aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATa                                                                    >  2:179351/1‑91 (MQ=255)
                        aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATa                                                                    >  1:1386537/1‑91 (MQ=255)
                        aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATa                                                                    >  1:111475/1‑91 (MQ=255)
                            cAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACtt                                                                <  2:1470751/91‑1 (MQ=255)
                            cAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACtt                                                                <  2:527573/91‑1 (MQ=255)
                             cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                               >  1:1351379/5‑91 (MQ=255)
                             cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                               >  2:3697781/5‑91 (MQ=255)
                             cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                               >  1:1351458/5‑91 (MQ=255)
                             cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                               >  1:357793/5‑91 (MQ=255)
                             aGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                               >  2:3668833/1‑91 (MQ=255)
                              gctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCa                                                              >  2:1712029/4‑91 (MQ=255)
                               cGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAt                                                             <  1:1076641/91‑1 (MQ=255)
                                ggTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAtt                                                            >  1:1274311/1‑91 (MQ=255)
                                  tCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAtttt                                                          >  2:76492/1‑91 (MQ=255)
                                    aTTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTAt                                                        >  2:881721/1‑91 (MQ=255)
                                     tttCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATc                                                       <  1:3381024/91‑1 (MQ=255)
                                     tttCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATc                                                       >  2:1824147/1‑91 (MQ=255)
                                     tttCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATc                                                       <  1:1059098/91‑1 (MQ=255)
                                     tttCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATc                                                       >  1:2138558/1‑91 (MQ=255)
                                       tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                     <  2:2395868/91‑1 (MQ=255)
                                       tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                     <  2:72149/91‑1 (MQ=255)
                                       tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                     <  1:2591392/91‑1 (MQ=255)
                                       tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                     <  1:876043/91‑1 (MQ=255)
                                       tCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCaa                                                     >  1:1273693/1‑91 (MQ=255)
                                        cAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAat                                                    >  1:587985/1‑91 (MQ=255)
                                         aGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAata                                                   <  2:609184/91‑1 (MQ=255)
                                         aGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAata                                                   >  2:1627511/1‑91 (MQ=255)
                                             atGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATc                                               <  1:1680506/91‑1 (MQ=255)
                                              tGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACtt                                                                >  1:542187/1‑73 (MQ=255)
                                              tGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCt                                              >  1:2452355/1‑91 (MQ=255)
                                                 aagccaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATTAATATATCTATg                                           <  1:2997747/87‑1 (MQ=255)
                                                  aatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGt                                          <  1:752582/91‑1 (MQ=255)
                                                        taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                                    >  2:1916414/1‑91 (MQ=255)
                                                         aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                                   >  2:3331035/1‑91 (MQ=255)
                                                         aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                                   <  1:2227695/91‑1 (MQ=255)
                                                         aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                                   <  1:1583364/91‑1 (MQ=255)
                                                          atTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGc                                  >  2:1132680/1‑91 (MQ=255)
                                                          atTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGc                                  <  2:1412187/91‑1 (MQ=255)
                                                          atTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGc                                  >  1:536335/1‑91 (MQ=255)
                                                           tttAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGcc                                 >  2:3235689/1‑91 (MQ=255)
                                                            ttAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCt                                >  1:3332718/1‑91 (MQ=255)
                                                            ttAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCt                                >  2:958414/1‑91 (MQ=255)
                                                             tAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCtt                               >  1:1868218/1‑91 (MQ=255)
                                                             tAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCtt                               <  1:692589/91‑1 (MQ=255)
                                                                aCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaa                            >  2:2446848/1‑91 (MQ=255)
                                                                 cATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaaa                           <  2:2242072/91‑1 (MQ=255)
                                                                 cATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaaa                           <  2:509900/91‑1 (MQ=255)
                                                                  aTCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAg                          <  2:492758/91‑1 (MQ=255)
                                                                   tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                         >  1:3871268/1‑91 (MQ=255)
                                                                   tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                         >  1:2642676/1‑91 (MQ=255)
                                                                     tcGCTCTAGTAGCAGTTCAACCGTTACTAGAGTTTTTTTGAAAATAACTTCATTTTATCGATATATCTATGTCATTATGCCTTAAAAGTGa                       <  2:2960256/91‑1 (MQ=255)
                                                                      cGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGAt                      <  2:3824084/91‑1 (MQ=255)
                                                                       gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa                     <  1:1585713/91‑1 (MQ=255)
                                                                       gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa                     >  2:2106807/1‑91 (MQ=255)
                                                                       gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa                     >  1:305769/1‑91 (MQ=255)
                                                                       gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa                     <  2:9694/91‑1 (MQ=255)
                                                                         tctAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATaaa                   <  2:1332774/91‑1 (MQ=255)
                                                                          ctAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATc                                                       >  1:3338932/1‑54 (MQ=255)
                                                                           tAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat                 <  1:635980/91‑1 (MQ=255)
                                                                             gTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAatat               >  1:3527906/1‑91 (MQ=255)
                                                                             gTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAatat               <  2:1244620/91‑1 (MQ=255)
                                                                              tATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATAtt              >  2:739962/1‑91 (MQ=255)
                                                                               aTCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTa             >  2:1770518/1‑91 (MQ=255)
                                                                                 cAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAAc           >  1:3073633/1‑91 (MQ=255)
                                                                                 cAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAAc           <  2:3549577/91‑1 (MQ=255)
                                                                                   gTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTa         >  1:582217/1‑91 (MQ=255)
                                                                                    ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt        <  1:1906249/91‑1 (MQ=255)
                                                                                    ttCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAt        <  2:2021649/91‑1 (MQ=255)
                                                                                     tCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTAtt       >  1:3031080/1‑91 (MQ=255)
                                                                                      cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa      <  2:1569586/91‑1 (MQ=255)
                                                                                      cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa      <  1:3849218/91‑1 (MQ=255)
                                                                                      cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa      >  1:423393/1‑91 (MQ=255)
                                                                                      cAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTa      >  1:2561070/1‑91 (MQ=255)
                                                                                       aaCCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTaa     <  2:2160140/91‑1 (MQ=255)
                                                                                       aaCCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTaa     <  1:1744686/91‑1 (MQ=255)
                                                                                       aaCCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTaa     <  2:2939946/91‑1 (MQ=255)
                                                                                        aCCGATACTAGAGATTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAATGATAACATATTAACTGTTaaa    >  2:3134836/1‑91 (MQ=255)
                                                                                         ccGATACTAGAGTTTTTTTGATATTTATTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTaaaa   >  2:1750375/1‑91 (MQ=255)
                                                                                         ccGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTaaaa   >  1:3669394/1‑91 (MQ=255)
                                                                                          |                                                                                          
AGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACTATTAAAAT  >  USA300TCH1516_ALE/2563677‑2563857

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: