Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,544 T→C 44.4% Y43C (TAC→TGC)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,5440TC44.4% 24.9 / 91.5 78Y43C (TAC→TGC) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (27/16);  new base C (18/16);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

ACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTC  >  USA300TCH1516_ALE/1744405‑1744683
                                                                                                                                           |                                                                                                                                           
aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGATGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                                                                                                                                            >  1:589597/1‑141 (MQ=255)
aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGATGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAACTATGTCTTTTGAGTa                                                                                                                                            >  2:492252/1‑141 (MQ=255)
    aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGATCCAGGTCCAACAACGATATCAACATCTCTGTATTTTGCTTATTCt                                                                                                                                        >  1:638447/1‑141 (MQ=255)
    aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGCCAGCCAATGATAAAGCCACGACCGTCAGGTCCTACAACGATATCAACATCCCTGTCTTTTGCGCATTCt                                                                                                                                        >  2:108612/1‑141 (MQ=255)
      gaatgcACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGGCAATCCCCATGGAATAAGCTACAGGACAGCCAATGATAATGCCACGCGCTTCAGGTCGTACAACGAGATCAACATCTCTGTCTTTTGCGCATTCTAc                                                                                                                                      <  2:99820/136‑1 (MQ=255)
        cTTCACGAGGTAATTTCCCTTCTTTTCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGCCTTTTACGTATTCTACaa                                                                                                                                    >  2:635366/1‑141 (MQ=255)
           cacaaagTAATTTCTCTTCTTTTCTAACATGTGCAAAGCCAATCCTCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACTATATCAACATCTCTGTCTTTTGCGCATTCTACAAttt                                                                                                                                 <  2:330444/135‑1 (MQ=255)
           cACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAGTAAGCTACAGGACAGCCAATGTTAAAGCCACGCGCTTCAGGTCCTACAACGCTATCAAAATCTCTGTCTTTTGCGCATACAACAAtat                                                                                                                                 >  1:77887/1‑139 (MQ=255)
           cACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAttt                                                                                                                                 <  2:805107/141‑1 (MQ=255)
              gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTAt                                                                                                                              <  1:408263/141‑1 (MQ=255)
              gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTAt                                                                                                                              <  1:763882/141‑1 (MQ=255)
               aGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATc                                                                                                                             <  1:91442/141‑1 (MQ=255)
                          cttcttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGc                                                                                                                  <  2:449996/141‑1 (MQ=255)
                          cttcttTTCTAAACGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCAGAGc                                                                                                                  >  2:62659/1‑141 (MQ=255)
                            ttttttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCa                                                                                                                <  1:250602/139‑1 (MQ=255)
                            tcttTTCAAACAGGTGCAAAGCAAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCCTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCa                                                                                                                <  2:290809/141‑1 (MQ=255)
                              ttttCTAAAAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCata                                                                                                              <  2:760116/141‑1 (MQ=255)
                                ttCTAACACGTGCAAAGCCAACCCCCATCGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCTTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATg                                                                                                            <  1:655005/141‑1 (MQ=255)
                                  cTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATATCCATATGCt                                                                                                          >  1:427424/1‑141 (MQ=255)
                                     aCAGATGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGTGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCa                                                                                                       <  1:769344/141‑1 (MQ=255)
                                       aGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCAcc                                                                                                     >  2:796383/1‑141 (MQ=255)
                                        ggTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCa                                                                                                    >  2:254444/1‑141 (MQ=255)
                                        ggTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTACAGGTCCTACAACGACATAAACATAACTGTCTTTTGCGCATTCTACTATTTTATCTGTTGCATAGCCATATGCTTCACCa                                                                                                    <  2:651724/141‑1 (MQ=255)
                                            cAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCGATATGCTTCCCCATTAt                                                                                                >  1:867827/1‑141 (MQ=255)
                                            cAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGACTTCTCAAATATTATCTGTTGCATAGTCATATGCTTTACCAttct                                                                                                >  2:414451/1‑139 (MQ=255)
                                              aaGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGGTATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATcc                                                                                              <  2:804548/141‑1 (MQ=255)
                                              aaGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATcc                                                                                              <  1:614683/141‑1 (MQ=255)
                                              aaGCCAATCACCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGAGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGGCATATGCTTCACCATTAtta                                                                                              >  2:476910/1‑139 (MQ=255)
                                                gCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATccct                                                                                            >  1:604546/1‑139 (MQ=255)
                                                gCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCAt                                                                                            >  2:378848/1‑141 (MQ=255)
                                                        ccATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGTCATATGCTTCACCATTATCCATAATTgtag                                                                                    >  2:588492/1‑141 (MQ=255)
                                                         cATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTgtagt                                                                                   <  1:287908/141‑1 (MQ=255)
                                                             gAATAAGCTACAGGACAGCCATTGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAata                                                                               <  1:558676/141‑1 (MQ=255)
                                                               aTAAGCTACAGGACAGCCAATGATGATGGAACGCTCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATc                                                                             >  1:935232/1‑141 (MQ=255)
                                                                       cAGGACAACCAAAGATAAAGCCACGCGCTTAAGAGTCGACAAAGAAATAAAGATAACATACATATTCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAAc                                                                     >  2:665222/1‑141 (MQ=255)
                                                                        aGGACAGCCAATGATAAAGCCACGCGCTTCACGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTCCAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACt                                                                    >  2:308126/1‑141 (MQ=255)
                                                                         ggACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTa                                                                   >  1:429212/1‑141 (MQ=255)
                                                                          gACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTaa                                                                  >  2:428567/1‑141 (MQ=255)
                                                                           aCAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACTATTATCCATACTTGTAGTAATATCCTTTAAACTAac                                                                 >  1:499650/1‑141 (MQ=255)
                                                                             aGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAacac                                                               <  1:258946/141‑1 (MQ=255)
                                                                             aGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAaaac                                                               >  1:196992/1‑141 (MQ=255)
                                                                              gCCAATGATAAAGCCACGCGCTTCCGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACAcc                                                              <  2:762576/141‑1 (MQ=255)
                                                                                   tGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGtt                                                         >  2:485466/1‑141 (MQ=255)
                                                                                        aaGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCgtc                                                    >  2:474453/1‑139 (MQ=255)
                                                                                        aaGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATACTTGAAACTAACACCTGGTTTCGGc                                                    >  1:191771/1‑141 (MQ=255)
                                                                                         aGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGcc                                                   >  2:554256/1‑141 (MQ=255)
                                                                                         aGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACACTATTCTCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAAGAACTGGTTTCGGcc                                                   >  1:362267/1‑141 (MQ=255)
                                                                                         aGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGcc                                                   >  1:606315/1‑141 (MQ=255)
                                                                                          gCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCa                                                  >  1:789769/1‑141 (MQ=255)
                                                                                                cgcTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                            <  2:396263/141‑1 (MQ=255)
                                                                                                cgcTTCAGCTCCTACAACGATATCAACATCTCGGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTCTCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                            <  2:545522/141‑1 (MQ=255)
                                                                                                  cTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGaa                                          >  1:33966/1‑141 (MQ=255)
                                                                                                   ttCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAAc                                         <  1:622279/141‑1 (MQ=255)
                                                                                                    tCAGTTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACt                                        <  2:522617/141‑1 (MQ=255)
                                                                                                      aGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTc                                      >  2:933211/1‑141 (MQ=255)
                                                                                                       ggTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCt                                     <  2:126828/141‑1 (MQ=255)
                                                                                                         tCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGa                                   >  2:938369/1‑141 (MQ=255)
                                                                                                           cTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTTTGATa                                 >  2:210158/1‑141 (MQ=255)
                                                                                                              cAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGt                              >  1:536158/1‑141 (MQ=255)
                                                                                                                 cGATCTCAACATCTCTGTCTTCTGCCTTGTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTAtt                           <  2:946735/141‑1 (MQ=255)
                                                                                                                   atatCAACTTCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGc                         <  1:343371/141‑1 (MQ=255)
                                                                                                                   atatCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGc                         <  1:378848/141‑1 (MQ=255)
                                                                                                                   atatCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGc                         <  2:604546/141‑1 (MQ=255)
                                                                                                                        aaCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATTTTTTTTTTTTTTTTTTTTAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGGACTTCTGATCCGTATTTCTTTaa                    >  2:216835/1‑141 (MQ=255)
                                                                                                                        aaCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTaa                    >  1:160913/1‑141 (MQ=255)
                                                                                                                        aaCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTaa                    >  1:842642/1‑141 (MQ=255)
                                                                                                                        aaCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTaa                    >  1:804435/1‑141 (MQ=255)
                                                                                                                          cATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATTCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTTAACTTCTGAAACGTATTGCTTTAAAt                  >  1:877537/1‑141 (MQ=255)
                                                                                                                               ctGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCAtt             >  1:601352/1‑141 (MQ=255)
                                                                                                                               ctGTCTTTTGAGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCAtt             >  1:682715/1‑141 (MQ=255)
                                                                                                                                tGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTa            >  2:624932/1‑141 (MQ=255)
                                                                                                                                  tCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTCTTGCTTTAAATCCATTAat          <  1:613384/141‑1 (MQ=255)
                                                                                                                                  gcTTTTGCGTATTCTACAATGTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAat          <  2:650384/140‑1 (MQ=255)
                                                                                                                                   cTTCTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACGAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAata         <  1:130157/141‑1 (MQ=255)
                                                                                                                                    ttttGCGTTTTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAatat        <  2:397908/141‑1 (MQ=255)
                                                                                                                                    ttttGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAatat        >  2:78048/1‑141 (MQ=255)
                                                                                                                                        gcgcATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTtcc    <  1:607602/141‑1 (MQ=255)
                                                                                                                                          gTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTTTTACCCTTGAAACTGACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTtcctc  >  1:179556/1‑141 (MQ=255)
                                                                                                                                          gTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTtcctc  >  2:629212/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
ACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTC  >  USA300TCH1516_ALE/1744405‑1744683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: