Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,566,052 A→G 7.6% intergenic (+101/+39) nhaK_2 → / ← plaP Sodium, potassium, lithium and rubidium/H(+) antiporter/Low‑affinity putrescine importer PlaP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,566,0520AG 7.6% 788.1 / 38.3 253intergenic (+101/+39)nhaK_2/plaPSodium, potassium, lithium and rubidium/H(+) antiporter/Low‑affinity putrescine importer PlaP
Reads supporting (aligned to +/- strand):  ref base A (131/100);  new base G (8/11);  total (141/112)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 265 total aligned reads displayed.

ACAATACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGATA  >  USA300TCH1516_ALE/2565962‑2566141
                                                                                          |                                                                                         
aCAATACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGaaaa                                                                                           >  2:1339266/1‑91 (MQ=255)
aCAATACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGaaaa                                                                                           >  2:11928/1‑91 (MQ=255)
 cAATACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAg                                                                                          >  1:1620911/1‑91 (MQ=255)
   aTACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTa                                                                                        <  1:1073423/91‑1 (MQ=255)
    tACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTaa                                                                                       <  2:657521/91‑1 (MQ=255)
    tACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTaa                                                                                       <  1:209679/91‑1 (MQ=255)
      cccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAAt                                                                                     >  2:2292992/1‑91 (MQ=255)
      cccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAAt                                                                                     <  2:2485515/91‑1 (MQ=255)
       ccccGTACATATGATACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATAATCAGCAACAAAAAATCAGATTGATGAAAAGTAAATa                                                                                    <  2:3803850/91‑1 (MQ=255)
       ccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATa                                                                                    <  1:1138223/91‑1 (MQ=255)
       ccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATa                                                                                    <  1:832065/91‑1 (MQ=255)
        cccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATaa                                                                                   >  1:183837/1‑91 (MQ=255)
           gTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCt                                                                                <  2:3642380/91‑1 (MQ=255)
                tatGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTtcat                                                                           >  1:2588005/1‑91 (MQ=255)
                 atGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTtcatc                                                                          <  1:3792032/91‑1 (MQ=255)
                  tGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATACAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTtcatca                                                                         <  2:1546796/91‑1 (MQ=255)
                   gagaCAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCaa                                                                        >  2:2215128/1‑91 (MQ=255)
                    agaCAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAAt                                                                       <  1:2077304/91‑1 (MQ=255)
                    agaCAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAAt                                                                       <  1:86814/91‑1 (MQ=255)
                       cAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAGGTTATTTACTTTTCATCAATctg                                                                    <  1:1464697/91‑3 (MQ=21)
                        aaaGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCGGATTGATGAAAGGTTATTTACTTTTCATCAATctga                                                                   <  2:1716861/91‑4 (MQ=18)
                          aGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCGGATTGATGAAAGGTTATTTACTTTTCATCAATCTGAtt                                                                 <  1:2708078/91‑1 (MQ=17)
                          aGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGAtt                                                                 >  1:3594731/1‑91 (MQ=255)
                             catgatcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCGGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGa                                                              <  2:1065076/87‑1 (MQ=255)
                             cATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGa                                                              >  2:2316204/1‑91 (MQ=255)
                             cATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGa                                                              >  2:2082184/1‑91 (MQ=255)
                              aTTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGAt                                                             >  1:1593245/1‑91 (MQ=255)
                              aTTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGAt                                                             >  2:3292026/1‑91 (MQ=255)
                               ttATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGAtt                                                            >  1:2711322/1‑91 (MQ=255)
                                 tttgtctcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                          >  1:3019547/5‑91 (MQ=255)
                                 tttgtctcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                          >  2:2596962/5‑91 (MQ=255)
                                 atcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                          >  2:2819111/1‑91 (MQ=255)
                                 atcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                          >  2:2755801/1‑91 (MQ=255)
                                  tcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtata                                                         <  2:3303636/91‑1 (MQ=255)
                                  tcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtata                                                         <  2:1679935/91‑1 (MQ=255)
                                   catcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAg                                                        <  2:3123153/91‑1 (MQ=255)
                                    atcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGa                                                       >  1:614641/1‑91 (MQ=255)
                                      ctcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAAt                                                     <  2:1409453/91‑1 (MQ=255)
                                      ctcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAAt                                                     >  2:3494955/1‑91 (MQ=255)
                                       tcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCATTCCGATTTGATTATAGAATc                                                    <  1:1518361/91‑1 (MQ=255)
                                       tcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATc                                                    <  2:3346369/91‑1 (MQ=255)
                                        cATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCt                                                   <  1:2916765/91‑1 (MQ=255)
                                           atGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAtt                                                <  2:2130823/91‑1 (MQ=255)
                                           atGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAtt                                                <  2:3555645/91‑1 (MQ=255)
                                           atGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAtt                                                >  1:2695288/1‑91 (MQ=255)
                                              tACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAttttt                                             >  2:1359881/1‑91 (MQ=255)
                                               aCGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAtttttt                                            <  2:1420836/91‑1 (MQ=255)
                                                  gggTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAg                                         >  1:2171769/1‑91 (MQ=255)
                                                  gggTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAg                                         >  2:3167152/1‑91 (MQ=255)
                                                   ggTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTCTTTAAGt                                        >  1:1948520/1‑91 (MQ=255)
                                                    gTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGtt                                       <  2:2637433/91‑1 (MQ=255)
                                                    gTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGtt                                       >  1:3794740/1‑91 (MQ=255)
                                                      ttttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTa                                     >  1:2658824/1‑91 (MQ=255)
                                                      ttttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTa                                     <  1:3716826/91‑1 (MQ=255)
                                                       tttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTaa                                    >  1:2349645/1‑91 (MQ=255)
                                                        ttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTaaa                                   >  2:3850965/1‑91 (MQ=255)
                                                        ttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTaaa                                   >  1:1699073/1‑91 (MQ=255)
                                                         tataTTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAAt                                  >  2:1768591/1‑91 (MQ=255)
                                                         tataTTCAACATCAAAAAATCAGAGTGATGAAAAGTAAATAACCTTTCATCAATCCGACTTGATTACAGAATCTATTTTTTAAGTTTAAAt                                  >  2:1769151/1‑91 (MQ=255)
                                                          atatTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATg                                 <  1:1376793/91‑1 (MQ=255)
                                                           tatTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATgg                                <  2:1856376/91‑1 (MQ=255)
                                                             ttCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAATTGGaa                              >  2:1075931/1‑91 (MQ=255)
                                                             ttCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGaa                              >  2:2433843/1‑91 (MQ=255)
                                                              tCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAAt                             <  2:1909663/91‑1 (MQ=255)
                                                               cAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAAtt                            >  1:2600149/1‑91 (MQ=255)
                                                                 aCATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGt                          <  1:3420835/91‑1 (MQ=255)
                                                                  cATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTa                         >  1:2481538/1‑91 (MQ=255)
                                                                    tCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACa                       >  1:1173179/1‑91 (MQ=255)
                                                                    tCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACa                       >  1:1911245/1‑91 (MQ=255)
                                                                      atcaaaTCGGATTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATa                     >  2:3790133/4‑91 (MQ=17)
                                                                      aaaaaaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGCTTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATa                     >  1:2791423/1‑91 (MQ=255)
                                                                      aaaaaaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATa                     >  1:2137365/1‑91 (MQ=255)
                                                                         aaaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGt                  >  2:1728474/1‑91 (MQ=255)
                                                                          aaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGtt                 <  2:2695288/91‑1 (MQ=255)
                                                                          aaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGtt                 <  1:2936436/91‑1 (MQ=255)
                                                                           aTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTa                <  1:3293735/91‑1 (MQ=255)
                                                                            tCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTaa               >  1:2276642/1‑91 (MQ=255)
                                                                            tCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTaa               >  2:874909/1‑91 (MQ=255)
                                                                            tCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTaa               <  1:503077/91‑1 (MQ=255)
                                                                             cAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAAc              <  2:502732/91‑1 (MQ=255)
                                                                              ggATTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACa             >  2:241101/2‑91 (MQ=18)
                                                                              aGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTg                           <  1:3010962/77‑1 (MQ=255)
                                                                              aGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACa             <  1:417257/91‑1 (MQ=255)
                                                                               gATTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt            >  1:757460/1‑91 (MQ=21)
                                                                               gATTGATGAAAGGTTATTTACTTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTTAGTTTATTTGGAATTGTACATACGTTAACAt            >  2:322679/1‑91 (MQ=18)
                                                                               gATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt            <  2:2831828/91‑1 (MQ=255)
                                                                               gATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATGTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt            >  2:3060321/1‑91 (MQ=255)
                                                                               gATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATACAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt            >  1:3076767/1‑91 (MQ=255)
                                                                                aTTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAtt           >  2:2026564/1‑91 (MQ=25)
                                                                                aTTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAtt           >  2:1281082/1‑91 (MQ=21)
                                                                                aTTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAtt           >  2:3741048/1‑91 (MQ=255)
                                                                                aTTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAtt           >  2:972544/1‑91 (MQ=255)
                                                                                 ttGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTc          >  1:1006692/1‑91 (MQ=255)
                                                                                 ttGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTc          <  2:2973885/91‑1 (MQ=255)
                                                                                  tgatgaAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCt         >  1:3214511/1‑91 (MQ=255)
                                                                                   gatgaAACGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCtt        <  2:3451189/91‑1 (MQ=255)
                                                                                   gatgaAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCtt        <  2:366732/91‑1 (MQ=255)
                                                                                     tgaAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCtttt      <  1:261072/91‑1 (MQ=255)
                                                                                       aaaaGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGa    <  1:988306/91‑1 (MQ=255)
                                                                                         aaGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGAta  <  1:1424781/91‑1 (MQ=255)
                                                                                          |                                                                                         
ACAATACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGATA  >  USA300TCH1516_ALE/2565962‑2566141

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: