Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,563,759 A→G 9.0% intergenic (+24/‑114) 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,7590AG 9.0% 730.3 / 44.5 244intergenic (+24/‑114)cycA_2/nhaK_2D‑serine/D‑alanine/glycine transporter/Sodium, potassium, lithium and rubidium/H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base A (117/105);  new base G (4/18);  total (121/123)
Fisher's exact test for biased strand distribution p-value = 2.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 249 total aligned reads displayed.

TGCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACT  >  USA300TCH1516_ALE/2563669‑2563849
                                                                                          |                                                                                          
 gCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAg                                                                                           <  1:332439/91‑1 (MQ=255)
 gCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAg                                                                                           >  1:518860/1‑91 (MQ=255)
 gCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAg                                                                                           >  2:2015081/1‑91 (MQ=255)
 gCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAg                                                                                           >  1:2235795/1‑91 (MQ=255)
 gCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAg                                                                                           <  1:1721231/91‑1 (MQ=255)
 gCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAg                                                                                           >  2:2147845/1‑91 (MQ=255)
  cTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGt                                                                                          <  2:3953654/91‑1 (MQ=255)
   tGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGtt                                                                                         <  1:3631296/91‑1 (MQ=255)
   tGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGtt                                                                                         <  2:1883374/91‑1 (MQ=255)
   tGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGtt                                                                                         <  2:1989939/91‑1 (MQ=255)
     aaaaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCa                                                                                       <  1:643013/91‑1 (MQ=255)
     aaaaGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCa                                                                                       <  1:2646237/91‑1 (MQ=255)
        aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                    <  1:1092653/91‑2 (MQ=255)
        aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                    <  1:3659694/91‑2 (MQ=255)
        aGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAAct                                                                                    <  1:1625704/91‑2 (MQ=255)
         gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAActg                                                                                   <  1:144372/91‑3 (MQ=255)
         gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCg                                                                                   >  1:3518830/1‑91 (MQ=255)
         gCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCg                                                                                   >  1:448090/1‑91 (MQ=255)
          cTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCGGTTGAActga                                                                                  <  1:3557210/91‑4 (MQ=255)
          cTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGa                                                                                  >  1:2964576/1‑91 (MQ=255)
           tCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGAt                                                                                 >  1:3269465/1‑91 (MQ=255)
           tCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGAt                                                                                 <  2:1481165/91‑1 (MQ=255)
            cgcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATa                                                                                >  1:1375669/1‑91 (MQ=255)
             gcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATAc                                                                               >  2:1722134/1‑91 (MQ=255)
             gcgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATAc                                                                               <  1:291590/91‑1 (MQ=255)
              cgTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACt                                                                              <  2:2972378/91‑1 (MQ=255)
               gTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTa                                                                             <  1:3600999/91‑1 (MQ=255)
                tATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTag                                                                            <  2:1368627/91‑1 (MQ=255)
                tATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTag                                                                            <  1:2853629/91‑1 (MQ=255)
                 aagAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTaga                                                                           <  2:2345687/89‑1 (MQ=255)
                  tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                          <  1:400379/91‑1 (MQ=255)
                  tGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTagag                                                                          <  1:1293341/91‑1 (MQ=255)
                   gAACAGCGACAAAATGACAGCGGTCATCTCGGATATGACAATCCATAATGTAGACGTCTCGCGCTAGTATCAGATCAACCGATACTAGAGt                                                                         <  2:377592/91‑1 (MQ=255)
                   gAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGt                                                                         >  2:215544/1‑91 (MQ=255)
                    aacaacGACAAAATGACTGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtt                                                                        <  1:3387618/91‑1 (MQ=255)
                    aacaacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtt                                                                        <  1:3909611/91‑1 (MQ=255)
                      caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                      <  1:1909334/91‑1 (MQ=255)
                      caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                      >  2:2766209/1‑91 (MQ=255)
                      caacGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGtttt                                                                      >  1:3734995/1‑91 (MQ=255)
                         cGACAAAATGACAGCGGTCATTTCGGATATGACAATCAATAATTTAGACATCTCGCTCTAGTATCAGTTCAACCGAGACTAGAGttttttt                                                                   <  2:878594/91‑1 (MQ=255)
                         cGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGttttttt                                                                   <  2:3187368/91‑1 (MQ=255)
                          gACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTg                                                                  <  2:3624270/91‑1 (MQ=255)
                           aCAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGa                                                                 >  2:2933092/1‑91 (MQ=255)
                              aaaTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaaaa                                                              >  1:1119760/1‑91 (MQ=255)
                               aaTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAAt                                                             >  1:1355151/1‑91 (MQ=255)
                               aaTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAAt                                                             <  1:1816310/91‑1 (MQ=255)
                                aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGaaattt                                                            <  2:51763/91‑4 (MQ=255)
                                aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATa                                                            >  1:3550506/1‑91 (MQ=255)
                                aTGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATa                                                            <  2:2700665/91‑1 (MQ=255)
                                  gACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAAc                                                          >  1:3180641/1‑91 (MQ=255)
                                   aCAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACt                                                         <  1:1135235/91‑1 (MQ=255)
                                   aCAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACt                                                         >  1:835982/1‑91 (MQ=255)
                                     cgctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTc                                                       >  2:2538616/5‑91 (MQ=255)
                                     aGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTATAGTCTTTTTGAAAATAACtca                                                       >  2:2625580/1‑89 (MQ=255)
                                      gctgTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCa                                                      >  1:2372351/4‑91 (MQ=255)
                                       cGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCAt                                                     >  1:205743/1‑91 (MQ=255)
                                            aTTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTAt                                                <  1:748362/91‑1 (MQ=255)
                                                 aGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAata                                           <  1:2437440/91‑1 (MQ=255)
                                                 aGATATGACAACCAATAATTTAAACACCTCGCTCTAGTATCAGTTCAACCGCTAATAGAGTTTTTTTGTAAATAACTTCATTTTCTAAata                                           >  2:3516844/1‑91 (MQ=255)
                                                  gATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAATCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAatat                                          <  2:3048889/91‑1 (MQ=255)
                                                   atatGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTCAAAATTATCTCTTTTTATTAttata                                         >  2:2423957/1‑91 (MQ=255)
                                                    tatGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAatatat                                        <  2:22280/91‑1 (MQ=255)
                                                     atGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATc                                       >  1:59512/1‑91 (MQ=255)
                                                          aatcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGt                                  <  1:2267186/91‑1 (MQ=255)
                                                           atcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTc                                 >  2:1303988/1‑91 (MQ=255)
                                                            tcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCa                                <  1:2479449/91‑1 (MQ=255)
                                                            tcaatAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCa                                >  1:991586/1‑91 (MQ=255)
                                                              aataatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCAtt                              <  1:2702119/91‑1 (MQ=255)
                                                               ataatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTTTAATTATCTTCATTTTATCAATATATCTATGTCATTa                             >  2:1362362/1‑91 (MQ=255)
                                                                taatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTAt                            >  1:1748871/1‑91 (MQ=255)
                                                                 aatTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATg                           <  1:850262/91‑1 (MQ=255)
                                                                  atTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGc                          <  1:3025159/91‑1 (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)
                                                                       aaCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaa                     <  2:2800623/91‑1 (MQ=255)
                                                                        aCATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTaaa                    >  1:2858438/1‑91 (MQ=255)
                                                                          aTCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAg                  <  1:1242294/91‑1 (MQ=255)
                                                                           tctcGCTCTCGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                 <  1:3779313/91‑1 (MQ=255)
                                                                           tctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGt                 <  2:64052/91‑1 (MQ=255)
                                                                            ctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTg                <  1:1526810/91‑1 (MQ=255)
                                                                            ctcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTg                >  2:601112/1‑91 (MQ=255)
                                                                             tcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGa               <  1:567786/91‑1 (MQ=255)
                                                                             tcGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGa               <  2:1646571/91‑1 (MQ=255)
                                                                              cGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGAt              <  1:1737857/91‑1 (MQ=255)
                                                                              atctctAGTATCGGTTGAACTGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGAt              >  2:2719207/3‑91 (MQ=38)
                                                                               gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa             <  1:1934018/91‑1 (MQ=255)
                                                                               gCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATa             >  1:3290063/1‑91 (MQ=255)
                                                                                ctctAGTATCAGTTCAACCGATACTAGAGTTTTTTTGTAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATaa            >  2:1814722/1‑91 (MQ=255)
                                                                                ctctAGTATCAGTTCAACCGATACTAGAGTTTTGTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTGAAGAGTGATaa            <  2:547575/91‑1 (MQ=255)
                                                                                   tAGTATCGGTTGAACTGATACTAGAGTCTTTTTGCAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAat         >  2:1535876/1‑91 (MQ=38)
                                                                                      tATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGGTAAAATAtt      >  1:902769/1‑91 (MQ=255)
                                                                                       aTCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTa     <  1:1433057/91‑1 (MQ=255)
                                                                                       aTCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAATATTaa    >  2:1914185/1‑91 (MQ=255)
                                                                                        tCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTaa    >  1:2568852/1‑91 (MQ=255)
                                                                                        tCAGTTCAACCGATACTAGAGTTTTTCTGAAAATAACTTCATTTTATCAATATACCTATGTCATTATGTCTTACAAGTGATAAAACATTaa    >  2:3196052/1‑91 (MQ=255)
                                                                                        tCAGTTCAACCGATACCAGAGTTTTCTTGACAATAACTTCACTTCATCAATCTCTCTATCTCATTATGCCTTAAAAGCGATAAAatatcac    >  2:3373442/1‑88 (MQ=255)
                                                                                         cAGTTCAACCGATACTAGAGTTTTTTTGGAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAAc   <  2:3035469/91‑1 (MQ=255)
                                                                                         cAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAAc   <  1:661769/91‑1 (MQ=255)
                                                                                          |                                                                                          
TGCTGAAAAGCTCGCGTATGAACAACGACAAAATGACAGCGGTCATTTCAGATATGACAATCAATAATTTAAACATCTCGCTCTAGTATCAGTTCAACCGATACTAGAGTTTTTTTGAAAATAACTTCATTTTATCAATATATCTATGTCATTATGCCTTAAAAGTGATAAAATATTAACT  >  USA300TCH1516_ALE/2563669‑2563849

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: