Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 52.9% 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,6430AG52.9% ‑1.7 / 83.2 68intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (20/12);  new base G (22/14);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512508‑512782
                                                                                                                                       |                                                                                                                                            
aaattaagctccttagcgacaatggtagtcaaacttacgttccgctagagcagaacgttgccaggcagttttataaacaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa                                                                                                                                         >  2:393462‑M1/113‑141 (MQ=255)
  aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  1:101211/1‑141 (MQ=255)
  aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  2:55384/1‑141 (MQ=255)
  aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  2:339281/1‑141 (MQ=255)
           cttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta                                                                                                                              <  1:146724‑M1/40‑1 (MQ=255)
            ttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaatataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaa                                                                                                                             >  1:55724‑M1/101‑141 (MQ=255)
               gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc                                                                                                                          <  1:178679‑M1/44‑1 (MQ=255)
                    cAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTc                                                                                                                     >  1:170968/1‑141 (MQ=255)
                      gTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCgg                                                                                                                  <  2:371319/141‑1 (MQ=255)
                       tttGAAAGCCTGTATAGTATTTATTTGTTCAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                 <  1:295305/141‑1 (MQ=255)
                        gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAACATGTATAATTAATTCTTGTCGGTa                                                                                                                <  1:226840‑M1/53‑1 (MQ=255)
                        gtagtcgaacttacgttccgatagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATCAAATGTATAATTAATTCTTGTCGGTa                                                                                                                <  1:241224‑M1/53‑1 (MQ=255)
                        aagAAACCCTGTATATTATTTATTTGTGGAGGCAAAGAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTCCATACAACCGTAATAAAATGTATAATTAATCCTTGTTGGTa                                                                                                                >  2:180452/3‑141 (MQ=255)
                          attcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACGTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                              >  2:87972‑M1/87‑141 (MQ=255)
                           ggcgaacttacgttccgctagagtagaacgttgccatgcagttttttaatcaaatttttgttaaaaaataaaatggacaagataaAAAAAGTTATTGACGTAAATGTTGATAAAATGTATGATTATTTCTTGTCGGTAAGa                                                                                                             >  2:14844‑M1/86‑141 (MQ=255)
                           aaaGCATGTACAGTATTCATTTGCTGAGGCAAACAAAACAACTCAACTTAAGAAATAAATTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                             <  1:27522/141‑1 (MQ=255)
                            aaGCCTGGATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAAAGTTACTAGAAGTTTAAATTAATTCTTGTCGGAAAGaa                                                                                                            >  2:180308/1‑141 (MQ=255)
                             cgaactttcgttcctctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaa                                                                                                           >  2:177295‑M1/84‑141 (MQ=255)
                                 tGTATAGTATTTATTTGTTGAGGCCCACACAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAATTATTGACTTAAATATCATTAAAATGTGTAATTAATACTTGTAGATAAAAAAAAta                                                                                                       >  2:357316/1‑140 (MQ=255)
                                   tataGTATTTATTTGACGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa                                                                                                     >  1:199503/1‑141 (MQ=255)
                                    acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaACAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                    <  1:268819‑M1/65‑1 (MQ=255)
                                            tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            <  1:335990/141‑1 (MQ=255)
                                            ctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            >  1:22170‑M1/69‑141 (MQ=255)
                                             tagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                           >  2:173111‑M1/68‑141 (MQ=255)
                                             atttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                           <  1:396517/141‑1 (MQ=255)
                                              agagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                          <  2:39134‑M1/75‑1 (MQ=255)
                                              agagtagaacgttgccaagcagttttttaaacaaattttgtttaaaaaataaaacggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                          >  2:59247‑M1/67‑141 (MQ=255)
                                              agagtagaaccttcccaggaagtattttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTAGTCGATAAGAAAAATGAAAATTGAAAACt                                                                                          >  2:36574‑M1/67‑141 (MQ=255)
                                               ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTACCGAAAATTAATTTTAAAAAGTTATTGACTTAACTGTTAATAAAATGTATAATTCATTCTTGTCGGTAAGAAAAATGAACATTGacaccta                                                                                         >  1:171080/1‑135 (MQ=255)
                                                tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTGATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                        >  1:240100/1‑141 (MQ=255)
                                                   agaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                     <  1:259509‑M1/80‑1 (MQ=255)
                                                     aGGCAAACAAAACAACTCAAATTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAc                                                                                   >  1:242803/1‑141 (MQ=255)
                                                      ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                  <  1:55660/141‑1 (MQ=255)
                                                      agCAAAAAAAACACAGCAACCTAAGAAATAACTTGAATTAATAACGAAAATTAATTTTAAAAAGTTATTGACTTAACTGTTAATAAAATGTATAAATAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                  <  1:18365/140‑1 (MQ=255)
                                                      acgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                  <  1:362693‑M1/83‑1 (MQ=255)
                                                          aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATAAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAAATAATACTGGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                              >  2:236260/1‑141 (MQ=255)
                                                          aacaaaacaaCTAAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTAAAAAAGTTATTGACATAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                              <  1:379398/141‑1 (MQ=255)
                                                             caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                           <  1:336196‑M1/90‑1 (MQ=255)
                                                             aaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTTTCGGTAAGAAAAATGAACATTGAAAACTGAATGACACTATGTc                                                                           >  2:304754/1‑141 (MQ=255)
                                                              aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                          >  1:336075‑M1/51‑141 (MQ=255)
                                                              aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                          >  2:104387/1‑141 (MQ=255)
                                                                acaacTCTACGTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                        <  2:275467/141‑1 (MQ=255)
                                                                      ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                  >  1:6710‑M1/43‑141 (MQ=255)
                                                                       tttaatcaaattttgggtaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAACTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                 >  2:208051‑M1/42‑141 (MQ=255)
                                                                        actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                                >  2:336159/1‑141 (MQ=255)
                                                                         ctTAAGCAATAACTTGAATTACTAACGAAGATCAATTTTGAAAAGTTATTGACGTAAATGTTAATAAAATGTATAATAAATTCTTGTCGAGGAACAAGAAGAACATCGAAAACTGAATGGCAATATGAACAAGTTAATTcc                                                               >  2:187427/1‑141 (MQ=255)
                                                                           tAAGAAATAACTAGAATCACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                             <  1:42772/141‑1 (MQ=255)
                                                                             aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                           >  2:39061/1‑141 (MQ=255)
                                                                              aaattttggtaaaaaaaaaaaagggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                          >  2:313426‑M1/35‑141 (MQ=255)
                                                                               aattttggttaaaaaataaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGGTAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATACCaaaaaa                                                         >  2:125444‑M1/34‑141 (MQ=255)
                                                                                aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTACATATTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGGAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc                                                        >  1:132897/1‑141 (MQ=255)
                                                                                      ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                  >  1:339612/1‑141 (MQ=255)
                                                                                       ttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAGGTAACta                                                 <  2:61712‑M1/116‑1 (MQ=255)
                                                                                        taaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtat                                                >  1:365954‑M1/25‑141 (MQ=255)
                                                                                         aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                               >  1:33317‑M1/24‑141 (MQ=255)
                                                                                         aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                               >  1:366736‑M1/24‑141 (MQ=255)
                                                                                          aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATGAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGATAAAAAAAAAAAACGTAACTATaa                                              <  2:170968/141‑1 (MQ=38)
                                                                                               taaaatgtacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAGTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACTGAACTATAAGTTAc                                         >  2:86418‑M1/18‑141 (MQ=255)
                                                                                                aaCGAAAATTAATTGTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                        <  1:32619/141‑1 (MQ=255)
                                                                                                  aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGCCAACTGAGTGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                      <  1:91714‑M1/127‑1 (MQ=255)
                                                                                                   atggacaagataaAAAAAGTTATTGACTTCAATGTTGATAAAATGTATAATTAATTCTGGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAAGATGTCAAGGTTAATTCCAAAAAACGTAAATAAAAGTTACAAAc                                     <  2:9460‑M1/128‑1 (MQ=255)
                                                                                                       acaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                 >  2:24370‑M1/10‑141 (MQ=255)
                                                                                                           gatcaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                             <  2:131630‑M1/136‑1 (MQ=255)
                                                                                                             taaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                           >  2:178679‑M1/4‑141 (MQ=255)
                                                                                                                 aaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt                       >  2:287812/1‑141 (MQ=255)
                                                                                                                          tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa              <  2:259509/141‑1 (MQ=255)
                                                                                                                               tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAAATTAATTAGTATTTATGAGCAAATCgaa         >  2:186896/1‑141 (MQ=18)
                                                                                                                                   tGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAACCATTATTTAGTAGTTTTGAGCTAATCCACcatc     >  1:146092/1‑141 (MQ=18)
                                                                                                                                      taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa  <  2:397906/141‑1 (MQ=255)
                                                                                                                                      taataaAATGAATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGATAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa  <  2:74232/141‑1 (MQ=18)
                                                                                                                                       |                                                                                                                                            
TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512508‑512782

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: