Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 47.0% intergenic (+711/‑6052) recR → / → speA_2 Recombination protein RecR/Arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE512,6430AG47.0% 177.0 / inf 740intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (222/168);  new base G (185/161);  total (408/332)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 773 total aligned reads displayed.

AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512503‑512783
                                                                                                                                            |                                                                                                                                            
aGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCATGTATAGTATTTATTTGTTGAGGCAAACAAAACAACACAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                                                                              <  2:3166648/141‑1 (MQ=255)
 ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                             >  1:554014/1‑141 (MQ=255)
  gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            <  2:1387909/141‑1 (MQ=255)
  gcgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            <  2:682237/139‑1 (MQ=255)
    gtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTCTAAAAAGCTATTGACTTAAATGTtaataa                                                                                                                                          >  2:3706572/1‑141 (MQ=255)
    atAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGATGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                                                                          <  1:1412967/140‑1 (MQ=255)
      aaTTTTTACTGTTGTTAAGCAGTTTGAAAGCATGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaaa                                                                                                                                        >  2:502063/1‑141 (MQ=255)
          aagctccttagcgtcgatggtagtcgaacttacgttccgccagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaantaagatggacaatataaAAAAAGTTATTGACTTAAGTGTTGATAAAATGta                                                                                                                                    <  2:421220‑M1/34‑1 (MQ=255)
            gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                  <  1:500874‑M1/36‑1 (MQ=255)
             ctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatctggttcttcactaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaa                                                                                                                                 <  1:1026452‑M1/37‑1 (MQ=255)
               ggttgttAATCAGTTTTAAAGCATTTATAGTATTTAGTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAATTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                                                                               <  2:90957/140‑1 (MQ=255)
                 ttgttAAGGAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGTCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                                                                             <  2:39148/141‑1 (MQ=255)
                  tagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaat                                                                                                                            <  1:3763218‑M1/42‑1 (MQ=255)
                  tagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaat                                                                                                                            <  1:306375‑M1/42‑1 (MQ=255)
                             ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCACACCAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATCAATTTTAAAAAGTTATTGACTTACATGTTACTCAAATCTATAATTAATTCTTGTCGGTa                                                                                                                 >  2:617584/1‑141 (MQ=255)
                              tagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTaa                                                                                                                >  1:2399849‑M1/88‑141 (MQ=255)
                               gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTATCGGTaaa                                                                                                               >  2:2350390/1‑140 (MQ=255)
                               agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                               <  1:2000157‑M1/55‑1 (MQ=255)
                                ggcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                              >  1:445416‑M1/86‑141 (MQ=255)
                                aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                              >  1:3811036/1‑141 (MQ=255)
                                 aaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                             >  1:2840280/1‑141 (MQ=255)
                                   gaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                                                                           >  1:2975880‑M1/83‑141 (MQ=255)
                                      cttacgttccgctagagtagaacgttgccaggcagttaattaaataatttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg                                                                                                        >  1:824009‑M1/80‑141 (MQ=255)
                                         acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGGAAAATGAAc                                                                                                     <  1:3577764‑M1/65‑1 (MQ=255)
                                         acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                     <  2:3213926‑M1/65‑1 (MQ=255)
                                            gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTCAGAAAAATGAACAtt                                                                                                  >  2:1761666/1‑141 (MQ=255)
                                             tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 >  1:1745301‑M1/73‑141 (MQ=255)
                                             tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 <  2:3795751/141‑1 (MQ=255)
                                             gccactagattagaacgttgccaggcagtgttttaataaaattttagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCCGTAAGAAATATGAACATTg                                                                                                 >  2:555225‑M1/73‑141 (MQ=255)
                                             aacgatagagtagaaagttgacaggaagtatattaatcaaattatagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 >  2:2559653‑M1/73‑141 (MQ=255)
                                              ccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGa                                                                                                <  1:1863477‑M1/70‑1 (MQ=255)
                                              atttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATTAACATTgc                                                                                                >  1:2496861/1‑140 (MQ=255)
                                               tttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaa                                                                                               >  1:651059/1‑141 (MQ=255)
                                                 atagagtagaacgataaaaggcagttttataatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                             >  2:19210‑M1/69‑141 (MQ=255)
                                                   tttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAACATGAACATTGAAAACt                                                                                           >  2:4150691/1‑141 (MQ=255)
                                                    ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                          >  2:2702753/1‑141 (MQ=255)
                                                      tgagaacgtggacagggagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACAGaa                                                                                        >  2:1840019‑M1/64‑141 (MQ=255)
                                                      gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                                        <  1:2137123/141‑1 (MQ=255)
                                                        tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                      >  1:3151649/1‑141 (MQ=255)
                                                         gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                     <  2:1099245‑M1/81‑1 (MQ=255)
                                                         gAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                     >  2:3651570/1‑141 (MQ=255)
                                                           ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                   >  2:2783096/1‑141 (MQ=255)
                                                            cgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGCATGACaa                                                                                  <  2:460102‑M1/84‑1 (MQ=255)
                                                             gttgccaggcagttttttaatcacattttggataaaaaataaaatgcacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGGCGTCACGAAAAATGAACATGGCACAGTGCATGacagt                                                                                 <  2:3554271‑M1/85‑3 (MQ=255)
                                                             gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAat                                                                                 <  2:2110031‑M1/85‑1 (MQ=255)
                                                              ttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                                >  1:2067835‑M1/56‑141 (MQ=255)
                                                              ttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                                <  1:3653733‑M1/86‑1 (MQ=255)
                                                                   aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                           >  1:2300534‑M1/51‑141 (MQ=255)
                                                                     acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAATCAGAAAAACACTTTGTCaaa                                                                         >  2:1964934/1‑140 (MQ=255)
                                                                          tcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTAAATCAAATGTATAAAAAATTCTTGTCGGTAAGAAAAATGAACATCGAAAACTGAATGACAATATGCCAACGTTaa                                                                    >  2:1578886/1‑141 (MQ=39)
                                                                           caactTAAGAAATAACTTGAATTACTAACGAAAAGTAATTTTAAAAAGTTCTTGACGAAAATGTGAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                   <  1:3221555/141‑1 (MQ=38)
                                                                            aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                  >  2:3394024/1‑141 (MQ=255)
                                                                             actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                                 >  2:3284306/1‑141 (MQ=255)
                                                                              taatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                                <  1:4292067‑M1/102‑1 (MQ=255)
                                                                               ttAAGAAAGAACGTGAATTACAAAAGAAAATTAAGTTTAAAAGGTTATGGACTGAAATGTGAATAAAATGTATAAATAATTATTGTCGGTAAAAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                               <  2:1657564/141‑1 (MQ=25)
                                                                               aatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                               >  2:4298952‑M1/39‑141 (MQ=255)
                                                                                 aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                             <  1:1197425/141‑1 (MQ=255)
                                                                                 aaGAAATAACTTGAATTACGAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                             <  2:407179/141‑1 (MQ=255)
                                                                                 aaGAAATAACATGAAATACTATAGAAAAATAAAATAAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAGTATGTCAACGTTAATTCCaaa                                                             <  2:3978022/141‑1 (MQ=37)
                                                                                    aattgtggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                          >  1:2583278‑M1/34‑141 (MQ=255)
                                                                                         tggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                     <  1:1291746‑M1/113‑1 (MQ=255)
                                                                                          ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    >  2:1036469‑M1/28‑141 (MQ=255)
                                                                                          ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    <  1:1455743‑M1/114‑1 (MQ=255)
                                                                                          ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    <  1:3469350‑M1/114‑1 (MQ=255)
                                                                                           ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   <  1:1544231/141‑1 (MQ=255)
                                                                                           ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATAAACATTGAAAACTGAATGACACTATGTCAACATTAATTCCAAAAAAAGTAACt                                                   >  2:120467/1‑141 (MQ=255)
                                                                                           ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTCAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   >  1:956346/1‑141 (MQ=255)
                                                                                            ttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACta                                                  <  1:3344910‑M1/116‑1 (MQ=255)
                                                                                              aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                                >  1:3289162‑M1/24‑141 (MQ=255)
                                                                                               aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                               >  1:712820/1‑141 (MQ=255)
                                                                                                ttACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtataca                                              >  2:3464745/1‑139 (MQ=255)
                                                                                                     aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         >  1:1580520‑M1/17‑141 (MQ=255)
                                                                                                     aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         >  1:394432‑M1/17‑141 (MQ=255)
                                                                                                       cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  2:2647445/1‑141 (MQ=255)
                                                                                                        aaaaaTTAATTTTAAAAAGTTAATGACGTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                      <  1:2761957/140‑1 (MQ=255)
                                                                                                             ttaagtTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTAAAAACTGAATGACAATATGTCAACGTTAATTCAACAAAGTGTAGAGATAAGTTACAAAcatcat                                 >  2:4026403/1‑141 (MQ=25)
                                                                                                               agataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               <  1:903856‑M1/135‑1 (MQ=255)
                                                                                                               aattTTAAACACATATTATCTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               >  2:3867495/1‑141 (MQ=21)
                                                                                                               aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               <  1:2204874/141‑1 (MQ=255)
                                                                                                                attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                              >  1:3628379/1‑141 (MQ=255)
                                                                                                                 ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAAATTATTTAg                             >  2:1500916/1‑141 (MQ=255)
                                                                                                                  gaaAAAAAGTTATTGAGTGAAATGGTGATAAAATGTATAATTAATTCCTGTCGGTAAGAAAAATGAACATTCAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                            >  1:4034016‑M1/4‑141 (MQ=255)
                                                                                                                   ttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  1:1135807/141‑1 (MQ=255)
                                                                                                                    aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTTTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                          >  2:1329573‑M1/2‑141 (MQ=255)
                                                                                                                        aaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAATTAt                      >  2:2216676/1‑141 (MQ=255)
                                                                                                                         aGTTATTGACTTAAATGTAGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATg                     >  1:2909649/1‑141 (MQ=255)
                                                                                                                           ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                   <  2:2037718/141‑1 (MQ=255)
                                                                                                                             aTTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCTAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt                 >  2:836863/1‑141 (MQ=255)
                                                                                                                             aTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATACGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt                 <  1:4114106/141‑1 (MQ=255)
                                                                                                                              ttGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa                <  2:2955730/141‑1 (MQ=255)
                                                                                                                              ttGACTTAAATGGTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa                <  1:1946434/141‑1 (MQ=255)
                                                                                                                               tGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGACCTaa               >  1:2688960/1‑141 (MQ=255)
                                                                                                                               tGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGATAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTCATGAACTaa               >  2:2830893/1‑141 (MQ=255)
                                                                                                                                   ttAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaa           >  1:3024583/1‑141 (MQ=255)
                                                                                                                                      aaTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAca        >  2:3050604/1‑141 (MQ=255)
                                                                                                                                       aTGTTGATAAAATGTATAATTAATACTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat       >  1:1373697/1‑141 (MQ=255)
                                                                                                                                         gTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca     <  1:261909/141‑1 (MQ=255)
                                                                                                                                          ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat    <  1:4205618/141‑1 (MQ=255)
                                                                                                                                          ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat    <  1:775414/141‑1 (MQ=255)
                                                                                                                                           taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATAAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATTAACTAATCAAAAATCATa   >  2:2447669/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512503‑512783

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