Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,751 A→T 23.1% L342* (TTA→TAA)  pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,372,7510AT23.1% 154.6 / 46.1 65L342* (TTA→TAA) pbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base A (37/13);  new base T (13/2);  total (50/15)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAG  >  USA300TCH1516_ALE/2372613‑2372888
                                                                                                                                          |                                                                                                                                         
aaGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAAt                                                                                                                                         >  2:410738/1‑141 (MQ=255)
aaGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCGGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTACAGAACAGCGCTAAt                                                                                                                                         >  2:212487/1‑141 (MQ=255)
      ggCACTGCACCTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGaaa                                                                                                                                   >  1:405080/1‑141 (MQ=255)
      ggCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGaaa                                                                                                                                   >  1:255877/1‑141 (MQ=255)
       gCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAaac                                                                                                                                  <  1:61300/141‑1 (MQ=255)
       gCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAaac                                                                                                                                  <  1:421388/141‑1 (MQ=255)
       gCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAaac                                                                                                                                  >  2:238322/1‑141 (MQ=255)
        cACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAaaca                                                                                                                                 >  2:51924/1‑141 (MQ=255)
        cACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAaaca                                                                                                                                 >  2:64932/1‑141 (MQ=255)
        cACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAaaca                                                                                                                                 >  2:396646/1‑141 (MQ=255)
        cACTGCAACTTCAAAACTTTTCCAGTTTATCGCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAaaca                                                                                                                                 <  2:30484/141‑1 (MQ=255)
              aaCTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAAcc                                                                                                                           >  1:198322/1‑141 (MQ=255)
                     aaaCTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACg                                                                                                                    >  2:327234/1‑141 (MQ=255)
                     aaaCTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACg                                                                                                                    >  1:314859/1‑141 (MQ=255)
                     aaaCTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACg                                                                                                                    >  1:199375/1‑141 (MQ=255)
                     aaaCTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACg                                                                                                                    >  2:37115/1‑141 (MQ=255)
                               aGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCaa                                                                                                          >  2:246499/1‑141 (MQ=255)
                               aGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCaa                                                                                                          >  2:325850/1‑141 (MQ=255)
                                   tttttcAGAAAAGTTTGCAGCCATTAATACGCCTACCACCACTAAGGCTGGTGTCGTAACTGCGCTTGTTACAACAGCGATTAACGGACTAAATAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATcc                                                                                                      >  2:158275/5‑141 (MQ=255)
                                   tATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATcc                                                                                                      >  1:25424/1‑141 (MQ=255)
                                   tATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATcc                                                                                                      >  2:214320/1‑141 (MQ=255)
                                   tATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATcc                                                                                                      >  2:385930/1‑141 (MQ=255)
                                      ctcAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGt                                                                                                   >  2:339385/1‑141 (MQ=255)
                                      ctcAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGt                                                                                                   >  2:333808/1‑141 (MQ=255)
                                               gTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCAcc                                                                                          >  1:422738/1‑141 (MQ=255)
                                               gTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCAcc                                                                                          >  1:166677/1‑141 (MQ=255)
                                                   gCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACa                                                                                      >  1:56263/1‑141 (MQ=255)
                                                   gCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACa                                                                                      >  2:43217/1‑141 (MQ=255)
                                                         aTTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAaca                                                                                >  2:314118/1‑141 (MQ=255)
                                                                cacCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAg                                                                         <  1:317531/141‑1 (MQ=255)
                                                                 acCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGt                                                                        <  1:428359/141‑1 (MQ=255)
                                                                 acCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGt                                                                        <  2:232813/141‑1 (MQ=255)
                                                                 acCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGt                                                                        <  2:107847/141‑1 (MQ=255)
                                                                  ccTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTTCAGCAACACCTGAAGtt                                                                       >  1:371108/1‑141 (MQ=255)
                                                                     accaccACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACAAAAGAACAGCGCTAATAGGAAACAAAAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCAGAAGTTGAt                                                                    >  2:320719/1‑141 (MQ=255)
                                                                      ccaccaCTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGAtt                                                                   >  1:411102/1‑141 (MQ=255)
                                                                      ccaccaCTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGAtt                                                                   >  1:327852/1‑141 (MQ=255)
                                                                      ccaccaCTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGAtt                                                                   >  2:111094/1‑141 (MQ=255)
                                                                             tAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGAtata                                                            <  2:289511/141‑1 (MQ=255)
                                                                                 gCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATaagaa                                                        >  2:267134/1‑141 (MQ=255)
                                                                                 gCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATaagaa                                                        >  1:101315/1‑141 (MQ=255)
                                                                                     gtgtCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAgttg                                                    >  1:387063/1‑141 (MQ=255)
                                                                                     gtgtCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAgttg                                                    >  1:387052/1‑141 (MQ=255)
                                                                                     gtgtCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAgttg                                                    >  1:292709/1‑141 (MQ=255)
                                                                                         cGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTtgttgt                                                >  1:94057/1‑141 (MQ=255)
                                                                                         cGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTtgttgt                                                >  2:215435/1‑141 (MQ=255)
                                                                                         cGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTtgttgt                                                >  2:208739/1‑141 (MQ=255)
                                                                                         cGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTtgttgt                                                >  1:347085/1‑141 (MQ=255)
                                                                                         cGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTtgttgt                                                >  1:226942/1‑141 (MQ=255)
                                                                                            aaCTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgt                                             >  2:307194/1‑141 (MQ=255)
                                                                                            aaCTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgt                                             >  1:373460/1‑141 (MQ=255)
                                                                                            aaCTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgt                                             >  2:375957/1‑141 (MQ=255)
                                                                                            aaCTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgt                                             >  2:397559/1‑141 (MQ=255)
                                                                                               tGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTTATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgttgt                                          <  1:140171/141‑1 (MQ=255)
                                                                                               tGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgttgt                                          <  1:238322/141‑1 (MQ=255)
                                                                                                gCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgttgtt                                         <  2:422738/141‑1 (MQ=255)
                                                                                                gCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgatgtt                                         <  2:166677/141‑1 (MQ=255)
                                                                                                    ttGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCaa                                     <  1:234381/141‑1 (MQ=255)
                                                                                                      gTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTATCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCaaaa                                   >  1:135689/1‑141 (MQ=255)
                                                                                                                  aTTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACa                       <  1:111094/141‑1 (MQ=255)
                                                                                                                       cGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGt                  >  1:44545/1‑141 (MQ=255)
                                                                                                                           cTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCt              >  1:148882/1‑141 (MQ=255)
                                                                                                                             aaaGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTaa            <  2:29550/141‑1 (MQ=255)
                                                                                                                                       gcTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTgag  >  1:220592/1‑141 (MQ=255)
                                                                                                                                       gcTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTgag  >  1:440157/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGAG  >  USA300TCH1516_ALE/2372613‑2372888

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