Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,566,055 A→T 7.8% intergenic (+104/+36) 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,0550AT 7.8% 843.8 / 42.7 254intergenic (+104/+36)nhaK_2/plaPSodium, potassium, lithium and rubidium/H(+) antiporter/Low‑affinity putrescine importer PlaP
Reads supporting (aligned to +/- strand):  ref base A (134/99);  new base T (9/11);  total (144/110)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 266 total aligned reads displayed.

ATACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGATAAATT  >  USA300TCH1516_ALE/2565965‑2566145
                                                                                          |                                                                                          
aTACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTa                                                                                            <  1:1073423/91‑1 (MQ=255)
  aCCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTaaa                                                                                          <  1:587546/91‑1 (MQ=255)
   cccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAAt                                                                                         >  2:2292992/1‑91 (MQ=255)
   cccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAAt                                                                                         <  1:3379879/91‑1 (MQ=255)
    ccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATa                                                                                        <  1:1138223/91‑1 (MQ=255)
    ccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATa                                                                                        <  1:832065/91‑1 (MQ=255)
    ccccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTAATAAAAAAAAAATCAGATTGATGAAAAGTAAATa                                                                                        <  2:3803327/91‑1 (MQ=255)
     cccGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATaa                                                                                       >  1:2110365/1‑91 (MQ=255)
        gTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCt                                                                                    <  2:3642380/91‑1 (MQ=255)
            atatGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACCTCAAAAAACCAGATTGATGAAAAGTAAATAACCTTtca                                                                                >  2:1343996/1‑91 (MQ=255)
             tatGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTtcat                                                                               >  1:2588005/1‑91 (MQ=255)
              atGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTtcatc                                                                              >  1:3531963/1‑91 (MQ=255)
              atGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTtcatc                                                                              <  2:217743/91‑1 (MQ=255)
               tGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTtcatca                                                                             >  2:3488426/1‑91 (MQ=255)
                 agaCAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAAt                                                                           <  1:86847/91‑1 (MQ=255)
                    cAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAGGTTATTTACTTTTCATCAATctg                                                                        <  1:1464697/91‑3 (MQ=21)
                     aaaGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGa                                                                       <  2:3464164/91‑1 (MQ=255)
                      aaGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCGGATTGATGAAAGGTTATTTACTTTTCATCAATCTGAt                                                                      <  2:1732870/91‑1 (MQ=18)
                       aGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCGGATTGATGAAAGGTTATTTACTTTTCATCAATCTGAtt                                                                     <  1:2708094/91‑1 (MQ=17)
                        gTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGAttt                                                                    >  2:2936436/1‑91 (MQ=255)
                          catgatcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCGGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGa                                                                  <  2:1065076/87‑1 (MQ=255)
                          cATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGa                                                                  >  2:2081779/1‑91 (MQ=255)
                           aTTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATACCCTTTCATCAATCCGATTTGAt                                                                 >  2:3131001/1‑91 (MQ=255)
                           aTTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGAt                                                                 >  1:1593245/1‑91 (MQ=255)
                           aTTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGAt                                                                 >  2:824679/1‑91 (MQ=255)
                            ttATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGAtt                                                                >  1:2711322/1‑91 (MQ=255)
                              tttgtctcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                              >  2:2596962/5‑91 (MQ=255)
                              tttgtctcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                              >  1:1773353/5‑91 (MQ=255)
                              atcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                              >  2:2755801/1‑91 (MQ=255)
                              atcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtat                                                              >  2:2819111/1‑91 (MQ=255)
                               tcatcTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATtata                                                             <  2:3303636/91‑1 (MQ=255)
                                   ctcATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAAt                                                         <  1:657526/91‑1 (MQ=255)
                                     cATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCt                                                       <  1:2917588/91‑1 (MQ=255)
                                      atatGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTa                                                      >  1:1537358/1‑91 (MQ=255)
                                       tatGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAt                                                     >  2:3664915/1‑91 (MQ=255)
                                        atGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAtt                                                    <  2:2130823/91‑1 (MQ=255)
                                        atGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAtt                                                    <  2:3555645/91‑1 (MQ=255)
                                        atGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAtt                                                    >  1:2695288/1‑91 (MQ=255)
                                           tACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAttttt                                                 >  1:302789/1‑91 (MQ=255)
                                           tACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTAttttt                                                 >  2:1359902/1‑91 (MQ=255)
                                               gggTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAg                                             >  2:3167152/1‑91 (MQ=255)
                                               gggTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAg                                             <  1:500544/91‑1 (MQ=255)
                                               gggTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAg                                             >  1:1644899/1‑91 (MQ=255)
                                                 gTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGtt                                           >  1:3794726/1‑91 (MQ=255)
                                                 gTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGtt                                           <  2:2637433/91‑1 (MQ=255)
                                                   ttttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTa                                         <  1:3716826/91‑1 (MQ=255)
                                                   ttttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTa                                         >  1:2658824/1‑91 (MQ=255)
                                                    tttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTaa                                        >  1:2349645/1‑91 (MQ=255)
                                                     ttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTaaa                                       >  2:1724421/1‑91 (MQ=255)
                                                     ttATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTaaa                                       >  2:3850965/1‑91 (MQ=255)
                                                      tataTTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACTTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAAt                                      <  2:1614008/91‑1 (MQ=255)
                                                      tataTTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAAt                                      >  2:1769100/1‑91 (MQ=255)
                                                       atatTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATg                                     <  1:1376793/91‑1 (MQ=255)
                                                       atatTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATg                                     <  1:1867166/91‑1 (MQ=255)
                                                        tatTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATgg                                    >  1:2040518/1‑91 (MQ=255)
                                                        tatTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATgg                                    <  2:1856376/91‑1 (MQ=255)
                                                          ttCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGaa                                  >  2:2433843/1‑91 (MQ=255)
                                                           tCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAAt                                 <  1:1917016/91‑1 (MQ=255)
                                                           tCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAAt                                 <  1:2086330/91‑1 (MQ=255)
                                                            cAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAAtt                                >  1:2600149/1‑91 (MQ=255)
                                                             aaCATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTg                               >  1:2618228/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)
                                                                  cAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAACGGAATTGTACAt                          >  2:3228962/1‑91 (MQ=255)
                                                                   atcaaaTCGGATTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATa                         >  2:3790133/4‑91 (MQ=17)
                                                                   aaaaaaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGCTTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATa                         >  1:2791423/1‑91 (MQ=255)
                                                                     aaaaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACg                       >  1:2140593/1‑91 (MQ=255)
                                                                      aaaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGt                      >  2:1205993/1‑91 (MQ=255)
                                                                      aaaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGt                      <  2:2515499/91‑1 (MQ=255)
                                                                       aaTCAGGTTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGtt                     >  2:2372744/1‑91 (MQ=255)
                                                                       aaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGtt                     <  1:2936436/91‑1 (MQ=255)
                                                                       aaTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGtt                     <  2:1414548/91‑1 (MQ=255)
                                                                        aTCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTa                    <  2:1698094/91‑1 (MQ=255)
                                                                         tCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTaa                   >  1:2276642/1‑91 (MQ=255)
                                                                         tCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAAGTGTACATACGTTaa                   <  2:3292890/91‑1 (MQ=255)
                                                                          cAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAAc                  <  2:2024180/91‑1 (MQ=255)
                                                                           ggATTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACa                 >  2:241101/2‑91 (MQ=18)
                                                                           ggATTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACa                 >  1:3779579/2‑91 (MQ=18)
                                                                           aGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTg                               <  1:3010962/77‑1 (MQ=255)
                                                                           aGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACa                 <  1:417257/91‑1 (MQ=255)
                                                                           aGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACa                 <  2:605165/91‑1 (MQ=255)
                                                                            gATTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt                >  1:757460/1‑91 (MQ=21)
                                                                            gATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt                >  2:2714511/1‑91 (MQ=255)
                                                                            gATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATGTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt                >  2:3060321/1‑91 (MQ=255)
                                                                            gATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATACAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAt                >  1:3076767/1‑91 (MQ=255)
                                                                             aTTGATGAAAGGTTATTTACTTTTCATCAATCTGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAtt               >  2:1281082/1‑91 (MQ=21)
                                                                             aTTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAtt               >  2:3741048/1‑91 (MQ=255)
                                                                             aTTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACAtt               >  2:972544/1‑91 (MQ=255)
                                                                              ttGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTc              <  2:2973885/91‑1 (MQ=255)
                                                                              ttGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTc              >  1:1006692/1‑91 (MQ=255)
                                                                               tgatgaAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCt             >  1:3214511/1‑91 (MQ=255)
                                                                                gatgaAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCtt            <  2:1504878/91‑1 (MQ=255)
                                                                                    aaaaGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGa        <  2:1163629/91‑1 (MQ=255)
                                                                                    aaaaGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGa        <  1:988306/91‑1 (MQ=255)
                                                                                      aaGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGAta      >  2:338007/1‑91 (MQ=255)
                                                                                      aaGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGAta      <  1:1424781/91‑1 (MQ=255)
                                                                                         tAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGAtaaat   >  2:1795383/1‑91 (MQ=255)
                                                                                         tAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGAtaaat   >  1:2259049/1‑91 (MQ=255)
                                                                                          aaaTAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGAtaaatt  >  2:194086/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATACCCCCGTACATATGAGACAAAGTCATTATCATCTCATATGTACGGGGTTTTTATATTCAACATCAAAAAATCAGATTGATGAAAAGTAAATAACCTTTCATCAATCCGATTTGATTATAGAATCTATTTTTTAAGTTTAAATGGAATTGTACATACGTTAACATTCTTTTGATAAATT  >  USA300TCH1516_ALE/2565965‑2566145

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: