Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 49.2% 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,6430AG49.2% 3.1 / 124.1 93intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (28/19);  new base G (27/18);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512503‑512782
                                                                                                                                            |                                                                                                                                           
aGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                                                                             <  2:410981/141‑1 (MQ=255)
 ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                            >  2:276605/1‑141 (MQ=255)
  gtttAATTTTTACTGTTGTAAAGCAGGTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAATTTATTGACTTAAATGTtaat                                                                                                                                           <  1:13320/138‑1 (MQ=255)
  gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                           >  2:446329/1‑141 (MQ=255)
   agaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTtgata                                                                                                                                          <  2:554099‑M1/27‑1 (MQ=255)
       gttaagctccttagcgtcgattgtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaagcaaactttggttgaaaaagaaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAt                                                                                                                                      >  1:532298‑M1/111‑141 (MQ=255)
       aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATACAAt                                                                                                                                      >  2:182228/1‑141 (MQ=255)
       aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                      >  2:464940/1‑141 (MQ=255)
        ttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg                                                                                                                                     >  2:32354‑M1/110‑141 (MQ=255)
            tACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtata                                                                                                                                 >  2:205616/1‑141 (MQ=255)
                gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                             <  1:275158/141‑1 (MQ=255)
                   gttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                                                                          <  2:358153/141‑1 (MQ=255)
                   agcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaaaataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATaattaatt                                                                                                                          >  2:498692‑M1/99‑141 (MQ=255)
                     cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt                                                                                                                        <  2:449269‑M1/45‑1 (MQ=255)
                        cgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAAttctggt                                                                                                                     <  1:499029‑M1/48‑4 (MQ=255)
                          aGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCg                                                                                                                   >  1:96732/1‑141 (MQ=255)
                            tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAATAAAAAGCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTTGGt                                                                                                                 >  2:118238/1‑141 (MQ=255)
                            ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                 <  1:357068‑M1/52‑1 (MQ=255)
                             ttGAAAGCCTGTAAAGTATTTATTTGTTGAGGCAAACCAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                <  1:153569/141‑1 (MQ=255)
                             gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                <  2:49262‑M1/53‑1 (MQ=255)
                             gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                <  2:283872‑M1/53‑1 (MQ=255)
                               gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTaaa                                                                                                              >  2:517971/1‑140 (MQ=255)
                               gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                              <  1:213245/141‑1 (MQ=255)
                                  cgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggtaaaaaaaaaaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaa                                                                                                           <  2:41342‑M1/58‑1 (MQ=255)
                                   gCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACTAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                                                                          <  2:60102/141‑1 (MQ=255)
                                     cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAAGTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                        >  2:269467/1‑141 (MQ=255)
                                        tacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa                                                                                                     >  2:46539‑M1/78‑141 (MQ=255)
                                          cgttccgctagaggagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                   >  2:185443‑M1/76‑141 (MQ=255)
                                            ttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATATGTTACATTATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                 >  1:555354‑M1/74‑141 (MQ=255)
                                            ttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATCTTGATCCACTGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                 >  2:422623‑M1/74‑141 (MQ=255)
                                             tatttattTGTTGAAGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                <  1:496330/141‑1 (MQ=255)
                                                ttatttGTTGAGGCAAACAAAACAACTCGACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaa                                                                                             >  1:363087/1‑141 (MQ=255)
                                                 tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            >  2:130315/1‑141 (MQ=255)
                                                 ctggagtagaaagttgccaggtagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            >  1:63792‑M1/69‑141 (MQ=255)
                                                 caagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            >  1:422712‑M1/69‑141 (MQ=255)
                                                  tagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaatagattgggcaagaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGATCATTGAAAAc                                                                                           >  1:408181‑M1/68‑141 (MQ=255)
                                                      gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                                       <  1:464940/141‑1 (MQ=255)
                                                        agaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                     >  2:418719‑M1/62‑141 (MQ=255)
                                                          aacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAc                                                                                   >  1:230814‑M1/60‑141 (MQ=255)
                                                          aGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAc                                                                                   <  1:446329/141‑1 (MQ=255)
                                                           ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                  >  2:262854/1‑141 (MQ=255)
                                                            gCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                 <  1:244392/141‑1 (MQ=255)
                                                            gCAAACAAAACAACTCAAATTCTCCAATGATTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAACATGAACATTGAAAACTGAATGACaa                                                                                 >  2:120372/1‑141 (MQ=255)
                                                            cgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                 >  1:124217‑M1/58‑141 (MQ=255)
                                                               aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                              >  2:410911/1‑141 (MQ=255)
                                                                cccaggaagttttttaatcaaatttaggataaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGTAAGAACAGTGAACATTGAAAACTGAATGACAATATg                                                                             <  2:393400‑M1/88‑1 (MQ=255)
                                                                 caaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                            <  1:426751/141‑1 (MQ=255)
                                                                   aggcagttttttaatcaaatgttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGGTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                          >  2:8004‑M1/51‑141 (MQ=255)
                                                                   aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                          <  1:560127/141‑1 (MQ=255)
                                                                     gcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATTACAATATGTCAAc                                                                        <  2:44302‑M1/93‑1 (MQ=255)
                                                                       aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                      >  2:535115/1‑141 (MQ=255)
                                                                         ttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAAGATGTGAACGTTa                                                                    <  1:184392‑M1/97‑1 (MQ=255)
                                                                           ttttaagaaaattatggttaaaaaataaaatagacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                  >  2:116799‑M1/43‑141 (MQ=255)
                                                                           caactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                  <  1:147854/141‑1 (MQ=255)
                                                                            tttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTcaacgtaaaaa                                                                 >  2:357068‑M1/42‑136 (MQ=255)
                                                                            aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                 >  2:26505/1‑141 (MQ=255)
                                                                            aactTAAGAAATAACTTGAATTACTAACGAAAATTAATCTTAAGAAGATATTGACTTAAATGTTAATAAAATGTATAATTAATTCTAGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGTTAAtt                                                                 >  2:513102/1‑141 (MQ=39)
                                                                               ttAACAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTTTCGCGGGAAATAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                              >  1:382702/1‑141 (MQ=38)
                                                                                  aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGATTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                           <  1:117067/141‑1 (MQ=255)
                                                                                   gAAATAACTCGAATTACTAACGAAAATTAATTTTAATAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                          >  1:411540/1‑141 (MQ=255)
                                                                                     aaTAACTTGAATTACTAACTAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc                                                        >  1:418959/1‑141 (MQ=255)
                                                                                      ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                       >  1:374102‑M1/32‑141 (MQ=255)
                                                                                      aTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                       <  1:511861/141‑1 (MQ=255)
                                                                                       tAACTTGAATTACTAACGAAAATTAATTTTAAAACGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGCCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                      <  2:43019/141‑1 (MQ=255)
                                                                                        tttgttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa                                                     >  2:522360‑M1/30‑141 (MQ=255)
                                                                                          gtgtaaaaactaatatgaagaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAATGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                   <  2:122804‑M1/114‑1 (MQ=255)
                                                                                          ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAAGTAAc                                                   <  2:466112‑M1/114‑1 (MQ=255)
                                                                                            ttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTAAACGTTAATTCAAAAAAACGTAACta                                                 <  2:86433‑M1/116‑1 (MQ=255)
                                                                                            ttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTGCAAAAAAAGTAACta                                                 <  1:209445‑M1/116‑1 (MQ=255)
                                                                                             gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTacatat                                                >  2:223023/1‑136 (MQ=255)
                                                                                             gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtat                                                >  2:58848/1‑141 (MQ=255)
                                                                                              aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAACTGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                               >  1:337979‑M1/24‑141 (MQ=255)
                                                                                               aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATCTATCAAAGTTAATTCCAAAAAACGTAACTATaa                                              >  2:402183/1‑141 (MQ=255)
                                                                                                 tACTAACGAAAATTAATTTTGAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                            >  1:342630/1‑141 (MQ=255)
                                                                                                    tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCACAGTTAATACCAAAAAACCTACCCAAAAGGTAc                                         >  1:549167/1‑141 (MQ=38)
                                                                                                    tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc                                         >  1:523610/1‑141 (MQ=255)
                                                                                                     aaaatggacaagataaAAAAATTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAATCTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                        >  1:254839‑M1/17‑141 (MQ=255)
                                                                                                     aaaaatgacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                        >  1:254842‑M1/17‑141 (MQ=255)
                                                                                                     aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                        <  1:276605/141‑1 (MQ=255)
                                                                                                        atggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACGAAc                                     <  1:306594‑M1/128‑1 (MQ=255)
                                                                                                        atggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAAATAAAAAAAAAAAAAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                     >  2:506507‑M1/14‑141 (MQ=255)
                                                                                                           aattaattTTAAAAAGTTATTGCATTAGTTGGTAAGAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt                                  >  1:557169/1‑141 (MQ=21)
                                                                                                             ttaattTTAAATAGTTTTTTACTTAAATGTTAATAAAATGTATAATTAATTCTCGGTGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGACAACGTTAATTCCAAAAAACGTAGCTATAAGTTACAAACattat                                <  2:307686/141‑1 (MQ=21)
                                                                                                               agataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                              >  1:469294‑M1/7‑141 (MQ=255)
                                                                                                               aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAAGACAATATGTCAACGTAAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                              <  1:269467/141‑1 (MQ=255)
                                                                                                                attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAAAAACAGTATTTa                             >  2:507132/1‑141 (MQ=255)
                                                                                                                 ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                            >  2:199800‑M1/5‑141 (MQ=255)
                                                                                                                  taaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGATACAAACATTATTTAGt                           <  1:486480‑M1/138‑1 (MQ=255)
                                                                                                                     aacaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATACGTCCACGTTGATGCCAAAAAACGTAACTATAAGTTACAAACATTCTTTAGTAtt                        >  2:277591/4‑141 (MQ=21)
                                                                                                                         aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACATAACTATAAGTTACAAACATTATTTATAATTTAta                    >  2:164150/1‑140 (MQ=255)
                                                                                                                           ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                  <  1:41342/141‑1 (MQ=255)
                                                                                                                              ttGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa               <  1:535279/141‑1 (MQ=255)
                                                                                                                               tGACTTAAATGTTTATAAAATGTATAATTAATTCTTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCACAAAACGCAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa              <  2:365338/141‑1 (MQ=255)
                                                                                                                               tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa              <  1:111940/141‑1 (MQ=255)
                                                                                                                               tGACTTAAATGTTAATAAAATGTATAATTAATTCTGGTCGGTAAAAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa              <  2:496886/141‑1 (MQ=255)
                                                                                                                                       aTGTTAATAAAATGTATAATTAATTCTAGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat      >  1:558646/1‑141 (MQ=255)
                                                                                                                                        tGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcata     >  2:92184/1‑140 (MQ=255)
                                                                                                                                         gTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca    >  1:59914/1‑141 (MQ=255)
                                                                                                                                           taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa  <  1:535115/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512503‑512782

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