Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 40.4% 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,6430AG40.4% 61.1 / 79.9 86intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base G (18/16);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512504‑512781
                                                                                                                                           |                                                                                                                                          
ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                           <  1:45883/141‑1 (MQ=255)
   gtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTCAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTAATTAACTTAAATGTtaataa                                                                                                                                        >  2:257430/1‑141 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                     >  2:119723/1‑141 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                     >  1:233237/1‑141 (MQ=255)
          agctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagtgttttaatcaaattttggttaaaaaataaaatggacaagataaAACAAGTTATTGACTTAAATGTTGATAAAATGtat                                                                                                                                 <  2:26439‑M1/35‑1 (MQ=255)
           gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                <  2:188837‑M1/36‑1 (MQ=255)
             cTGTTGTTAAGCAGTTTGACAGCCTGTATAGTATTTCTTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaat                                                                                                                              >  2:387586/1‑141 (MQ=255)
                   gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc                                                                                                                        <  1:86978‑M1/44‑1 (MQ=255)
                    tAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCt                                                                                                                       <  1:84910/141‑1 (MQ=255)
                    cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt                                                                                                                       >  1:78599‑M1/97‑141 (MQ=255)
                    cgtcaatagtagtcgaaattacgtaacgctagaaaagaacgttgacaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt                                                                                                                       >  2:323199‑M1/97‑141 (MQ=255)
                      tcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTCGATAAAATgtataattaattattt                                                                                                                     <  2:202307‑M1/47‑5 (MQ=255)
                        gatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTc                                                                                                                   <  2:348154‑M1/49‑1 (MQ=255)
                        cAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTc                                                                                                                   >  2:72560/1‑141 (MQ=255)
                         aGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCg                                                                                                                  >  1:364429/1‑141 (MQ=255)
                            ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                               >  1:263159/1‑141 (MQ=255)
                            ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATCCTTGTCGGTa                                                                                                               >  1:251297/1‑141 (MQ=255)
                               gtcgaacttacgtgccgctagagtagaacgttgccaagcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAACAATTTATTGACTTAAATTTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                            >  2:162466‑M1/86‑141 (MQ=255)
                               gtagaaattccgatcagctagagtagaacgttgtcaggcagttttttaatcaaattttggttaaaaaataaaataaaaaaaataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                            >  2:130230‑M1/86‑141 (MQ=255)
                                  gCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                                                                         >  1:111750/1‑141 (MQ=255)
                                    cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                       >  1:31684/1‑141 (MQ=255)
                                    acttacgtttcgctagagtagaacgttgctaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTCATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                       >  1:16713‑M1/81‑141 (MQ=255)
                                            tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                               <  1:43576/141‑1 (MQ=255)
                                            tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                               <  1:147040/141‑1 (MQ=255)
                                            tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                               <  1:320492/141‑1 (MQ=255)
                                            tatttatgTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                               <  1:33937/141‑1 (MQ=255)
                                               cctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaa                                                                                            >  1:395909‑M1/70‑141 (MQ=255)
                                                tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                           <  1:296720/141‑1 (MQ=255)
                                                atagagtagaacgatgacaagcagtttattaataaaatattggataaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                           >  2:358246‑M1/69‑141 (MQ=255)
                                                  tttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                         >  1:35615/1‑141 (MQ=255)
                                                   ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                        <  1:119723/141‑1 (MQ=255)
                                                    agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                       <  1:360929‑M1/77‑1 (MQ=255)
                                                        gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagctaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATTAACATTGGAAACTGAAtgt                                                                                   <  1:133740‑M1/81‑2 (MQ=255)
                                                        gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                   >  1:102274‑M1/61‑141 (MQ=255)
                                                                 caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTCAAAACTGAATGACAATATGTc                                                                          <  2:4329‑M1/90‑1 (MQ=255)
                                                                  aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAGAATGTATAATTAATTCTTGTCGGTAAGAAAAATGGACATTGAATACTTACTGACAATATGTCa                                                                         <  1:13917‑M1/91‑1 (MQ=255)
                                                                   ggcagttttttaatcaaattttggttaaaaaataaaatagacaaaataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                        >  1:284728‑M1/50‑141 (MQ=255)
                                                                    acaacTCAACTTAAGAAATAACTTGAATTACTACAGAAACTTCACCACAAAAAGCTCTTGTTTTTTGAAGGAAGAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                       >  1:93687/1‑141 (MQ=16)
                                                                      aactcaactTAAGAAATAACTTGAATTTTTAACGAAAATTAATTTTAAAAAGTTAATGACATAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                     <  1:96410/141‑1 (MQ=255)
                                                                        ttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTa                                                                   <  2:241657‑M1/97‑1 (MQ=255)
                                                                        ctcaactTAAGAAATAACTTGAATTACTAACGAACATTAATTTTAAAAAGGTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTa                                                                   >  1:16168/1‑141 (MQ=255)
                                                                          caactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                 >  1:80501/1‑141 (MQ=255)
                                                                           aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                <  1:72560/141‑1 (MQ=255)
                                                                            ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                               <  1:322745‑M1/101‑1 (MQ=255)
                                                                             ctcaataaaTAACTTGAATTACTAACGAAAATTAATTTTAAAAATTTATTGACTTAAATGTTAATAAAATGTATAATTAAGTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                              <  1:25525/135‑1 (MQ=255)
                                                                              aatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                             >  1:188837‑M1/39‑141 (MQ=255)
                                                                               tAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTGGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                            >  2:318152/1‑141 (MQ=255)
                                                                                tcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                           >  1:386650‑M1/37‑141 (MQ=255)
                                                                                taaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                           >  1:368385‑M1/37‑141 (MQ=255)
                                                                                  gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                         >  1:58250/1‑141 (MQ=255)
                                                                                   aaatatgtttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                        >  2:192500‑M1/34‑141 (MQ=255)
                                                                                     ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                      >  2:86978‑M1/32‑141 (MQ=255)
                                                                                       ttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa                                                    <  2:370531‑M1/112‑1 (MQ=255)
                                                                                             aaTTACTAACGAAAATTAATTCTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                              <  1:53158/141‑1 (MQ=255)
                                                                                              aaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGAACACGTTAGTCCTCTCTGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                             >  1:290116‑M1/23‑141 (MQ=255)
                                                                                              aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                             <  1:88924/141‑1 (MQ=255)
                                                                                               ttACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTGAACGTTAATTCCAAAAAACGTAACTATAAg                                            >  1:279933/1‑141 (MQ=255)
                                                                                                aaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATAGTCGGTAAGAAAAATGCAAATCAAAAAGTTAATGACAATATGTCAACGTTAATTAAAAAAAACATAACTATAaga                                           <  2:64858‑M1/121‑2 (MQ=255)
                                                                                                 aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                          >  1:120374/1‑141 (MQ=255)
                                                                                                 aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                          <  2:76373/141‑1 (MQ=255)
                                                                                                  cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                         >  2:24235/1‑141 (MQ=255)
                                                                                                    aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                       <  1:119901‑M1/125‑1 (MQ=255)
                                                                                                    aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                       <  1:396655‑M1/125‑1 (MQ=255)
                                                                                                    aaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                       >  2:214328‑M1/17‑141 (MQ=255)
                                                                                                    aaCGAAAATTAATTTTAAAAAGTTATTGAATTAAATGTTAATAAAATGTATAATTCATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTGTTTTTTTTTTCCAAAAAACGTAACTATAAGTTACa                                       <  2:123807/141‑1 (MQ=37)
                                                                                                      cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                     <  1:264485/141‑1 (MQ=255)
                                                                                                      cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                     <  2:35615/141‑1 (MQ=255)
                                                                                                      aaaggacaagaaaaAGAAAAGTATTGACGTAAAGGTAGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                     >  1:271159‑M1/15‑141 (MQ=255)
                                                                                                          gacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt                                 >  1:274940‑M1/11‑141 (MQ=255)
                                                                                                          aattaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTAAATCTCACAAACagt                                 >  2:137714/1‑139 (MQ=37)
                                                                                                              aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTTCAAAAAACGTAACTATAAGTTACAAACATTAttt                             <  2:67701/141‑1 (MQ=255)
                                                                                                              aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTCACTATAGGTGACAAACATTAttt                             >  1:247745/1‑141 (MQ=255)
                                                                                                              aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                             <  1:47129/141‑1 (MQ=255)
                                                                                                                ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAAACACAAAAAATATTTaa                           >  2:42339/1‑140 (MQ=25)
                                                                                                                 tttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATACCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          >  2:177890/1‑141 (MQ=255)
                                                                                                                 ttaataaaGTGATTGACGAAAAAGTAAATAGAAAGTATAATTAAATCGTGACNGNAAGAAAAATGAACATAGAAATCTGAATGACTATATGTCAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          <  2:183549/136‑1 (MQ=12)
                                                                                                                  atAAAAAGATATAGAATTAAATGTTAATAAAATGTATAATTAATTCTTGTAAGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                         <  2:111750/140‑1 (MQ=21)
                                                                                                                      aaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa                     <  2:373713/141‑1 (MQ=255)
                                                                                                                       aaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt                    >  2:322858/1‑141 (MQ=255)
                                                                                                                         attattaACTAAAATGTGAATAAAAAGTATAATAAAAACTTAGAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAAATAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                  <  2:82603/140‑1 (MQ=17)
                                                                                                                          ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                 <  1:295108/141‑1 (MQ=255)
                                                                                                                          ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                 >  2:152167/1‑141 (MQ=255)
                                                                                                                              tGAATTAAATGATAATAAAATGTAAAATTAATTCTTATCGGTAAGAAAAATGAAAATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa             <  2:246587/141‑1 (MQ=18)
                                                                                                                               gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt            >  1:113583/1‑141 (MQ=255)
                                                                                                                                  gtAAATGTGTATAGAATGTATAAGTAATTGTGGTCAGTAAGAAAAATGAACATTGAAAACTGAATGATAATATGTCAACGTTAATTCCAAAAAACGTAAAGATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaa         <  2:13917/140‑1 (MQ=17)
                                                                                                                                      aTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat     >  1:283612/1‑141 (MQ=255)
                                                                                                                                      aTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat     >  1:79378/1‑141 (MQ=255)
                                                                                                                                         ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat  <  2:93687/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512504‑512781

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: