Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 63.6% 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,6430AG63.6% 49.7 / 63.6 69intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (10/15);  new base G (24/20);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

GTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512505‑512783
                                                                                                                                          |                                                                                                                                             
gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            >  1:67937/1‑141 (MQ=255)
gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            >  1:208478/1‑141 (MQ=255)
   aagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa                                                                                                                                         >  1:237141‑M1/113‑141 (MQ=255)
           agccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaa                                                                                                                                 >  2:14797‑M1/105‑141 (MQ=255)
            cTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaat                                                                                                                                >  2:113407/1‑141 (MQ=255)
                aagcgtcgatggtagtcgaaattacgttcagctagagtagaacgatgacaggcagttttttaataaaatttaggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaat                                                                                                                            >  2:67622‑M1/100‑141 (MQ=255)
                     tcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTg                                                                                                                       <  2:39680‑M1/47‑1 (MQ=255)
                      gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGGTTTGGTAACAAAAATTAATTTCGCCAATCTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGt                                                                                                                      <  1:92957/141‑1 (MQ=255)
                      cgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttactcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGGCTTAAATGTTGATAAAATATATAATTAATTCTTGt                                                                                                                      <  2:235182‑M1/48‑1 (MQ=255)
                           gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                 <  1:341214‑M1/53‑1 (MQ=255)
                             agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                               >  1:195459‑M1/87‑141 (MQ=255)
                             agtcgaacttacgttccgctagaggagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                               >  1:170660‑M1/87‑141 (MQ=255)
                              gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                              >  1:245760‑M1/86‑141 (MQ=255)
                                cgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTTTCGATAAGaaa                                                                                                            <  1:305615‑M1/58‑1 (MQ=255)
                                cgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaa                                                                                                            <  1:305397‑M1/58‑1 (MQ=255)
                                    cttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg                                                                                                        <  2:126853‑M1/62‑1 (MQ=255)
                                          ttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                  <  1:209064‑M1/68‑1 (MQ=255)
                                          ttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                  >  1:258527‑M1/74‑141 (MQ=255)
                                          ttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                  >  1:279625‑M1/74‑141 (MQ=255)
                                           tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 >  1:351557‑M1/73‑141 (MQ=255)
                                            atttatttGTTGAGGCAAACAAAACAACTCAACTAAATAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtgta                                                                                                >  2:27664/1‑138 (MQ=255)
                                               atagagtaaaacgttggcaggaagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                             >  2:273184‑M1/69‑141 (MQ=255)
                                                 tttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                           <  2:296024/141‑1 (MQ=255)
                                                  ttattgAGGCAAACAAAACAACACAACTTAAGAAATAACTTGAATTAATAACGAAAATTAATTTTAAAAAGTTATTGACGTATATGTTAATAAAATGTATAATTAAGTCTTGTCGGTAAGAAAAATGAACATTGAAACCTg                                                                                          <  2:300902/138‑1 (MQ=255)
                                                   agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                         >  1:334284‑M1/65‑141 (MQ=255)
                                                   agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                         <  2:68887‑M1/77‑1 (MQ=255)
                                                   agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTAACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                         <  1:378148‑M1/77‑1 (MQ=255)
                                                      tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATACACTGTATAAATAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                      >  1:267090/1‑141 (MQ=255)
                                                      agaacgttgccaggcagttttttaatcaaatgttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                      >  1:132719‑M1/62‑141 (MQ=255)
                                                       gAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                     <  1:145859/141‑1 (MQ=255)
                                                        aacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAc                                                                                    <  2:132211‑M1/82‑1 (MQ=255)
                                                         acgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTATCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                   <  2:90449‑M1/83‑1 (MQ=255)
                                                             aacaaaacaaCTTAACATAAGAAATAACTTGAATTAATAACGAAAAATAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                               <  2:355554/141‑1 (MQ=255)
                                                                  ggcagttttttaatcaaattttggttaaaaaatcaactggacgtgataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                          >  1:202164‑M1/50‑141 (MQ=255)
                                                                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                         <  2:62011/141‑1 (MQ=255)
                                                                     aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                       >  2:284483/1‑141 (MQ=255)
                                                                         taactTAAGAAATATCTTTAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                   <  1:64522/140‑1 (MQ=255)
                                                                           ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                                 <  1:198487‑M1/101‑1 (MQ=255)
                                                                                ttttttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGTTAAATTGTATAATTAATTCTTTTTTTTTATAAAAATGAACATTGAAAACTGAATGTCAATATATCAACGTTAATTCCaaaa                                                            <  2:212168‑M1/106‑1 (MQ=255)
                                                                                     tAACATGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGATAATAAAATGTAGAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACACTGTTGCAATGTTAATTCCAAAAAACGt                                                       <  2:281690/141‑1 (MQ=38)
                                                                                     aaaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                       <  1:153000/140‑1 (MQ=255)
                                                                                      taggtaaaaaaataaaatggaacagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa                                                      >  1:126216‑M1/30‑141 (MQ=255)
                                                                                         ttGAATTAAGAACGAAAAGTAATTTTAAAAAGTTAGTGACTTAAATGTTAATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   <  2:57149/141‑1 (MQ=37)
                                                                                         gttaaaaaataaaatggaaatgataaGAAAAATTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   >  1:265290‑M1/27‑141 (MQ=255)
                                                                                             attgcTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                               >  2:219135/5‑141 (MQ=255)
                                                                                             aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                               <  1:101353/141‑1 (MQ=255)
                                                                                              aaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAg                                              >  1:275205‑M1/22‑141 (MQ=255)
                                                                                                 cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                           >  1:135524/1‑141 (MQ=255)
                                                                                                 ataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                           >  2:37755‑M1/19‑141 (MQ=255)
                                                                                                  tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTaa                                          >  1:30684/1‑140 (MQ=255)
                                                                                                   aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         <  2:25340‑M1/125‑1 (MQ=255)
                                                                                                   aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         <  1:103295‑M1/125‑1 (MQ=255)
                                                                                                    aaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                        <  1:170109‑M1/126‑1 (MQ=255)
                                                                                                     cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGACCATTCCTCCCTGAAAGTTGTTCTACAAAAGGTCAACGTTAATTCCAACAAACGTAACTATAAGTTACaaa                                       <  2:279547/141‑1 (MQ=21)
                                                                                                     aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  1:39680‑M1/15‑141 (MQ=255)
                                                                                                      atgtataagataaAAAAAATTATTTAATTAAATAGTAATAAAATGTATAATTAATTCTTGTTGATAAGAAAAATGAACATTTAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                      >  2:66177‑M1/14‑141 (MQ=255)
                                                                                                       tggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                     >  1:266457‑M1/13‑141 (MQ=255)
                                                                                                       tggacaagataaAAAAAGTTATTGACTCAAACGTTGAGTAATTTGTTATTTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                     >  1:266455‑M1/13‑141 (MQ=255)
                                                                                                        aaaGTAATTTTAAAAACTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAATTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACat                                    <  1:243982/141‑1 (MQ=255)
                                                                                                          acaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                  <  1:76085‑M1/132‑1 (MQ=255)
                                                                                                             agataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               >  1:29097‑M1/7‑141 (MQ=255)
                                                                                                             aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAACTTACAAACATTAttt                               >  2:316629/1‑141 (MQ=255)
                                                                                                             aattTTAAAAAGTGATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTAGGGGTGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               <  1:286027/141‑1 (MQ=21)
                                                                                                               ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                             >  1:144107/1‑141 (MQ=255)
                                                                                                               ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAATTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                             >  1:143580‑M1/5‑141 (MQ=255)
                                                                                                                  aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                          <  1:320694‑M1/140‑1 (MQ=255)
                                                                                                                      aaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt                      <  2:84160/141‑1 (MQ=255)
                                                                                                                          tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                  <  2:371860/141‑1 (MQ=255)
                                                                                                                          tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAATTAATATTCATGaac                  >  2:89626/1‑139 (MQ=255)
                                                                                                                              gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTAGGTAAGAAAAATGAACATTGAAATCTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt              <  2:76227/141‑1 (MQ=255)
                                                                                                                               aCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATCACTTCCAAATATTTTTTAGTATTTATGAGCTAAta             >  1:200479/1‑140 (MQ=21)
                                                                                                                                         taataaAATGTATAATTAATTCTTGTCGGAAAGAAAAAAGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACTTTATTTAGTATTTATGAGCTAATCAAACATCATa   <  1:143023/141‑1 (MQ=255)
                                                                                                                                          aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATaa  >  2:154442/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                             
GTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512505‑512783

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: