Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 677,733 A→G 5.2% I204V (ATT→GTT)  USA300HOU_RS03265 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE677,7330AG 5.2% 255.0 / 3.1 96I204V (ATT→GTT) USA300HOU_RS03265hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/50);  new base G (3/2);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GCGTTATACCAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACTATTCACTTGAT  >  USA300TCH1516_ALE/677593‑677864
                                                                                                                                             |                                                                                                                                     
gCGTTATTCCAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGa                                                                                                                                       <  1:129318/141‑1 (MQ=255)
gCGTTATACCAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGa                                                                                                                                       <  2:433188/141‑1 (MQ=255)
  gTTATACCAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGAtt                                                                                                                                     >  1:748395/1‑141 (MQ=255)
   ttATACCAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTc                                                                                                                                    >  2:708502/1‑141 (MQ=255)
    tataCCAATTACAAATTTTGGACATCAATTTTTCTTTTGCATTGGTATTATTTTAGTAATTGCCGACATATTTTATAAGCGACTTGAGAT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCa                                                                                                                                   >  2:396773/1‑141 (MQ=255)
     ttaCCAATTAAAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTACTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGCTTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCAt                                                                                                                                  <  1:526879/140‑1 (MQ=255)
     ataCCAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCAt                                                                                                                                  <  1:157149/141‑1 (MQ=255)
        ccAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCCTATGTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTGTTTTGTCAAAACACCGTCGATGCCATGGTTCATAAt                                                                                                                               >  1:63359/1‑141 (MQ=255)
        ccAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGGTTCATAAt                                                                                                                               >  1:850078/1‑141 (MQ=255)
        ccAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGGTTCATAAt                                                                                                                               >  1:8839/1‑141 (MQ=255)
          aaTTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAACAACTGCCGATGCAATGATTCATAATcc                                                                                                                             <  1:670973/141‑1 (MQ=255)
          aaTTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATTCAATGATTCATAATcc                                                                                                                             <  1:671547/141‑1 (MQ=255)
              aCAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACaa                                                                                                                         >  1:376599/1‑141 (MQ=255)
                aaaTTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATc                                                                                                                       >  1:723981/1‑141 (MQ=255)
                aaaTTTTGGAAACCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATc                                                                                                                       >  2:886055/1‑141 (MQ=255)
                 aaTTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCa                                                                                                                      >  1:527881/1‑141 (MQ=255)
                    tttGGACATCAATTTTTCTTTTTCATTGGTCTTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATACTCCACAATCAGcc                                                                                                                   >  1:78722/1‑141 (MQ=255)
                    tttGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGcc                                                                                                                   >  1:775157/1‑141 (MQ=255)
                     ttGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGccc                                                                                                                  <  2:748395/141‑1 (MQ=255)
                     ttGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGccc                                                                                                                  <  1:325632/141‑1 (MQ=255)
                      tGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCa                                                                                                                 <  1:843901/141‑1 (MQ=255)
                        gACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGa                                                                                                               <  1:241373/141‑1 (MQ=255)
                        gACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGa                                                                                                               >  2:705211/1‑141 (MQ=255)
                        gACATCAATTCTTCTTTGTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGa                                                                                                               <  1:24129/141‑1 (MQ=255)
                         aCATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTTTTTTGTCAAAAAACAGTCGATGCAATGATTCATAATCCACAATCAGCCCAGAt                                                                                                              >  2:278918/1‑141 (MQ=255)
                         aCATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGAt                                                                                                              >  1:722400/1‑141 (MQ=255)
                         aCATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGAt                                                                                                              <  2:314018/141‑1 (MQ=255)
                         aCATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGAt                                                                                                              >  1:536214/1‑141 (MQ=255)
                         aCATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATAAGCCCAGAt                                                                                                              >  2:190293/1‑141 (MQ=255)
                         aCATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCATCGATGCAATGATTCATAATCCACAATCAGCCCAGAt                                                                                                              >  1:802555/1‑141 (MQ=255)
                             cAAGTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGAttttt                                                                                                          <  1:379392/141‑1 (MQ=255)
                                tttttctttttcATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGCtttttttt                                                                                                       >  2:19657/1‑139 (MQ=255)
                                tttttctttttcATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCa                                                                                                       >  1:788382/1‑141 (MQ=255)
                                tttttctttttcATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCa                                                                                                       <  1:42702/141‑1 (MQ=255)
                                tttttctttttcATTGGTATTATTTTAGTAAGTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGTGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCa                                                                                                       <  2:781268/141‑1 (MQ=255)
                                tttttctttgtcATTGGTATTATTTTAGTAATCGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCa                                                                                                       <  1:195459/141‑1 (MQ=255)
                                tttttctatatcAATGGTATTAATTCACTAATTGCCGACATATTTTACAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATAAAAAAAAAAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCa                                                                                                       <  2:709551/141‑1 (MQ=255)
                                    tttttttcATTACTGTTGTTTTCTTAATTGCAGGCATTCTTAATAAGCCA‑TTCACTTCTTCGGGAACCGTGTTCTTATTTCGTCAAAAAACCTTCGATGCAATGATTCATAATCCACAATCCTCCCAGATTTTTTCattaa                                                                                                   <  2:441488/139‑1 (MQ=255)
                                    tctttttcATTGGTATTATTTTAGTAATTGCCGGTATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCattaa                                                                                                   <  1:537031/141‑1 (MQ=255)
                                     ctttttcATTGGTATTATTTTAGTAATTGCCGGCATATTTTATACGCGTCTGGTGTT‑GTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCattaat                                                                                                  <  1:252710/141‑1 (MQ=255)
                                      tttctcATGGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCattaatt                                                                                                 <  1:32965/141‑1 (MQ=255)
                                       ttttcATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCattaatta                                                                                                <  2:723981/141‑1 (MQ=255)
                                             ttGGCATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAATCCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGa                                                                                          <  1:64352/141‑1 (MQ=255)
                                               ggTATTATTTTAGAAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGata                                                                                        <  2:86355/141‑1 (MQ=255)
                                               ggTATAATTTTCGTAATGGCCGGCATATTTTATAAGCGACTTGAGGT‑TTCGGGATTCGGGTTAATAATTTGAAAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGata                                                                                        <  2:12400/141‑1 (MQ=255)
                                                     tttttACTAATTGACGCCATATTTTATAAGCGACTTGAGTT‑TTCGGGACTCGTGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAt‑tatta                                                                                 <  2:835700/140‑1 (MQ=255)
                                                     attTTAGGAATGGCCGTCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAt‑tatta                                                                                 <  2:239852/141‑1 (MQ=255)
                                                        ttAGTAATTGCAGGCATATTTTTTAAGCGTATTGAGTT‑TTCGGCAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGAAATATTCAGGAATTATATGGATAT‑TATtagta                                                                              <  2:417482/141‑1 (MQ=255)
                                                         tAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATCT‑TATtagtag                                                                             >  2:797332/1‑141 (MQ=255)
                                                         tAGGAATTGGCGCAATATGTTATAGGCGACTTGAGTT‑TTCGGGAATCGGGTCATTATTTTGTCGAAAAACCGTCGATGCAACGATTCATAGTCGACAACCGGCCCAGATTTTTTCCTTTATTATATGGATAT‑TATtagtgg                                                                             >  2:309316/1‑139 (MQ=255)
                                                          aGTATTTCCCGGCATATTTTATAAGCGACCTCAGTT‑TTCGGGAATCGAGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCATATTTTTTCATTACTTATATGGATAT‑TATTagtag‑t                                                                           <  1:248180/141‑1 (MQ=255)
                                                          aGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTagtag‑t                                                                           >  2:121167/1‑141 (MQ=255)
                                                          aGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGACTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTagtag‑t                                                                           <  1:312143/141‑1 (MQ=255)
                                                                 gCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGtt                                                                    >  1:14612/1‑141 (MQ=255)
                                                                  ccGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTta                                                                   <  1:885813/141‑1 (MQ=255)
                                                                  ccGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTta                                                                   <  1:221713/141‑1 (MQ=255)
                                                                  ccGCCATATTTTATAAGGGACTTAAGTT‑TTCTGGAATCGGGTTATAAAAAAAAAAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTta                                                                   <  2:732648/141‑1 (MQ=255)
                                                                    ggCATATTTTATACGCGACTTGAGTT‑TTCGGGAACCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTtata                                                                 <  1:3898/141‑1 (MQ=255)
                                                                    ggCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGCTTTTTTCATTAATTCTATGGCTCT‑TCCTCGTAG‑TTCTAGttctt                                                                 >  2:606238/1‑138 (MQ=255)
                                                                    ggCATATTTTATAAGCGACTCGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTtata                                                                 >  1:640582/1‑141 (MQ=255)
                                                                     gCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTTATTATATGGATAT‑TATTAGTAG‑TTCTAGTtatat                                                                <  1:309316/141‑1 (MQ=255)
                                                                         atTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATAttt                                                            <  2:479092/141‑1 (MQ=255)
                                                                         atTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATAttt                                                            <  1:64233/141‑1 (MQ=255)
                                                                          ttttATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTc                                                           >  1:28412/1‑141 (MQ=255)
                                                                            ttATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCCTCGATGCAATGATTCATAATCCACAATCAGCACAGATTTTTTCATTAATTATATGGATGT‑TATTAGGAG‑TTATGGTTAAATATTTcac                                                         >  1:410468/1‑141 (MQ=255)
                                                                            ttATAAGCGACTTGAGTT‑TTCGGGAAACGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTcac                                                         <  2:753698/141‑1 (MQ=255)
                                                                                   ggACTGGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGa                                                  <  1:729645/140‑1 (MQ=255)
                                                                                   cGACTTGAGTT‑TTCGGGAATGGGGTTATTATTTTGTCAAAAAACAGTCGATGCAATGACTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑GATTAGTAG‑TTCTAGTTATGTATTTCACAATTAGa                                                  >  1:551764/1‑141 (MQ=255)
                                                                                   cGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACGGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGa                                                  >  1:497181/1‑141 (MQ=255)
                                                                                   cGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGa                                                  <  1:326409/141‑1 (MQ=255)
                                                                                    gACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGAt                                                 <  1:356943/141‑1 (MQ=255)
                                                                                    gACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGAt                                                 >  1:814014/1‑141 (MQ=255)
                                                                                        tGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAACACCGTCGATGCAATGATTCATAAGCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTGG‑TTCTAGTTATATATTTCACAATTACATTAtt                                             >  1:396326/1‑140 (MQ=255)
                                                                                            tt‑ttCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCa                                         >  2:802437/1‑141 (MQ=255)
                                                                                             t‑ttCGGGAATCGGGTTATTATTTTTTCAAAAAACCGTCGATGCAATGATTCTTAACCCAACATCATCCCAGTTTTTTTCATTAATTATATGGATTT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCAc                                        <  2:214089/141‑1 (MQ=255)
                                                                                             t‑ttCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCAc                                        <  2:28412/141‑1 (MQ=255)
                                                                                             t‑ttCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCAc                                        <  2:634571/141‑1 (MQ=255)
                                                                                               ttCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGGTTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACg                                       <  2:850078/141‑1 (MQ=255)
                                                                                               ttCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGGTTCATAATACACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCttgtac                                       <  2:8839/141‑5 (MQ=255)
                                                                                               ttCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATCTGGATCT‑TCTTAGTAG‑ATCGATTTATATATTTCACAATTAGATTATCTTCACg                                       >  2:411322/1‑141 (MQ=255)
                                                                                                  gggAAACGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAGTTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTAc                                    <  2:518763/141‑1 (MQ=255)
                                                                                                   ggAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAACCCAGATTTTTTCATTAATTAAATGGATTTCTGTTA‑TAG‑ATCTCGTTATATATTTCACAATTAGATTATCTTCACTTACa                                   >  1:206864/1‑141 (MQ=255)
                                                                                                              attattTTTTTAAAAAACCGCCGATGCAATTTATCATTATGTACAAGCATTGCAGATTGTTTCATTAAGTATATGGATAT‑GATTAGTAG‑TTTTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAat                        <  2:26087/141‑1 (MQ=255)
                                                                                                              attattTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAat                        <  1:183317/141‑1 (MQ=255)
                                                                                                              attattTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCAAGCACCTTGCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAat                        <  1:662221/141‑1 (MQ=255)
                                                                                                              attattTTGTCAAAAAACCGTCGATGAAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAat                        <  2:410548/141‑1 (MQ=255)
                                                                                                                     tGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATaa                 >  1:53479/1‑141 (MQ=255)
                                                                                                                      gTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATaaa                <  1:768574/141‑1 (MQ=255)
                                                                                                                      aaaaaaaaaCCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATaaa                <  2:688545/138‑1 (MQ=255)
                                                                                                                         aaaaaaCCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTCGCTACTTG‑AATAGATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACt             >  1:209951/1‑141 (MQ=255)
                                                                                                                         aaaaaaCCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACt             >  1:8195/1‑141 (MQ=255)
                                                                                                                           aaaaCCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACTAt           <  2:896137/141‑1 (MQ=255)
                                                                                                                             aaCCGTCGATGCAATGATTCATAATCACCAATCCGCCCAGATTTTTTCATTAATTATATGGATCT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATATATCTGATAAAACTATTc         >  2:736241/1‑141 (MQ=255)
                                                                                                                              aCCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACTATTCa        >  1:564929/1‑141 (MQ=255)
                                                                                                                               ccGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACTATTCAc       >  2:694021/1‑141 (MQ=255)
                                                                                                                                 gTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACTATTCACtt     >  1:718828/1‑141 (MQ=255)
                                                                                                                                 ctCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACTATTCACtt     <  1:239840/140‑1 (MQ=255)
                                                                                                                                    gATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGCTAAAACTATTCACTTGAt  >  1:183624/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                     
GCGTTATACCAATTACAAATTTTGGACATCAATTTTTCTTTTTCATTGGTATTATTTTAGTAATTGCCGGCATATTTTATAAGCGACTTGAGTT‑TTCGGGAATCGGGTTATTATTTTGTCAAAAAACCGTCGATGCAATGATTCATAATCCACAATCAGCCCAGATTTTTTCATTAATTATATGGATAT‑TATTAGTAG‑TTCTAGTTATATATTTCACAATTAGATTATCTTCACGTACAAGATTATAAATATGATAAAACTATTCACTTGAT  >  USA300TCH1516_ALE/677593‑677864

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: ‑