Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 58.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,6430AG58.6% 40.3 / 61.9 65intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (8/19);  new base G (21/17);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 4.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512503‑512781
                                                                                                                                            |                                                                                                                                          
aGGTGTAATTTTTACCGGTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                                                                            <  2:222776/141‑1 (MQ=255)
 ggTGTAATTTTTACTGTTGTTAATCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAAATTTTCATGGTGATTGACTTAAATGTtaa                                                                                                                                           <  2:82641/141‑1 (MQ=255)
 ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                           <  2:5912/141‑1 (MQ=255)
 ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                           <  1:178054/141‑1 (MQ=255)
  gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                          <  2:226644/141‑1 (MQ=255)
     tattttttACTGTAGATAATAAGTTTGAAAGCCTGTATAGTATTTAATTGTTGAGTCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAAAAAAAAAAAAAAAGTTATTGACTTAAATGTTAATaaa                                                                                                                                       <  2:39350/138‑1 (MQ=255)
     aagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa                                                                                                                                       <  1:333362‑M1/29‑1 (MQ=255)
      agttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa                                                                                                                                      >  1:327477‑M1/112‑141 (MQ=255)
             aCTGTTGTAAAGCAGGTGGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAGACAACGCAACTTCAGAAATAACTTGAATTACTAACGAAAATGCACTTGAAAAAGTTATTGCCTTAAATGTTAATAAAATgcagaa                                                                                                                               >  2:131738/1‑136 (MQ=255)
              cTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATTTAtaat                                                                                                                              >  2:138612/1‑141 (MQ=255)
               ccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTCAATGTTAATAAAATGTAtaatt                                                                                                                             <  1:83554‑M1/39‑1 (MQ=255)
                 ttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgatcccgaggaaggttttaatcaaatttttgttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaa                                                                                                                           >  1:260729‑M1/101‑141 (MQ=255)
                     cgttgatggtagtcgaacttacgttccgctagagtagaacctggccaggtagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt                                                                                                                       >  2:29798‑M1/97‑141 (MQ=255)
                            ggtagtcgaacttacgttccggtagagtagaacgtggccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                <  1:296810‑M1/52‑1 (MQ=255)
                            ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                <  2:372662‑M1/52‑1 (MQ=255)
                             gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaagtaaaatggacaagataaAAAAAGTGATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                               <  2:254527‑M1/53‑1 (MQ=255)
                             gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                               <  1:133764‑M1/53‑1 (MQ=255)
                               gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                             <  1:84760/141‑1 (MQ=255)
                               agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                             <  2:219055‑M1/55‑1 (MQ=255)
                                         acgtgtcgctagagtagaacgttgcgaggcaggtatttaataaaattttggttaaaaaataaaagggacaagataaAAAAAGTTATTGACTTAAATGCTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                   >  2:228995‑M1/77‑141 (MQ=255)
                                         acggtcagatagagttggacgttgccaggcagttttttaatcaaattttggttaaaaaataaaataaaaaaaaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                   >  2:168827‑M1/77‑141 (MQ=255)
                                            ttccgctagagtagaacgttgccaggcagttttttttttatattttcgataaaaaataaaatggacaaaataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTGCGAAAAATaaacgtt                                                                                                <  2:13536‑M1/68‑4 (MQ=255)
                                            gTATATATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAATAATGAACAtt                                                                                                <  2:258849/141‑1 (MQ=255)
                                                 ctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                           >  1:192067‑M1/69‑141 (MQ=255)
                                                              ttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                              >  2:239219‑M1/56‑141 (MQ=255)
                                                              ttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                              >  2:239215‑M1/56‑141 (MQ=255)
                                                              aaacaaaacaACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                              >  1:218224/1‑141 (MQ=255)
                                                                  caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                          <  2:351796‑M1/90‑1 (MQ=255)
                                                                   aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                         >  1:1660‑M1/51‑141 (MQ=255)
                                                                    ggcagttttttaattaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATTATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAATCTGAATGACAATATGTCaa                                                                        >  2:355987‑M1/50‑141 (MQ=255)
                                                                     acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTACTTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                       >  2:318501/1‑141 (MQ=255)
                                                                          tcaactTAAAAAATAACGTTAATTACTAACGAAAAATAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                                                                  <  2:340827/141‑1 (MQ=255)
                                                                             acgtAAGAAATAACTTGAAGTACCAACGAAAATTAATTTTAAAAAGGTATTGCCTTAAATTGTAATAAAATGTATAATTAATTCTTGTCGTGAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTAAATTc                                                               <  2:163280/138‑1 (MQ=37)
                                                                              taatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCACCGTTAGTTcc                                                              <  1:185900‑M1/102‑1 (MQ=255)
                                                                               ttAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCACAGTTAATTCCa                                                             >  2:310689/1‑141 (MQ=255)
                                                                                tAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                            <  2:151067/141‑1 (MQ=255)
                                                                                    aattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                        >  1:219055‑M1/34‑141 (MQ=255)
                                                                                      aTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                      <  2:53949/141‑1 (MQ=255)
                                                                                       tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGTCAATATGTCAACGTTAATTCCAAAAAACGt                                                     <  2:280798‑M1/111‑1 (MQ=255)
                                                                                       tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                     <  2:160633‑M1/111‑1 (MQ=255)
                                                                                          ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                  <  2:244094‑M1/114‑1 (MQ=255)
                                                                                           ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCACAAAACGTAACt                                                 >  2:55083/1‑141 (MQ=255)
                                                                                           gttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                 <  2:57351‑M1/115‑1 (MQ=255)
                                                                                             taaaaaataaaacggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATCACAATATGTCAACGTTAATTCCAAAAAACGTAACtat                                               >  2:285130‑M1/25‑141 (MQ=255)
                                                                                             gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGAATTAAATGTTAATGAAATGTATAATTAATTCTTGTCGGTAAGAAAAGTGAACGTTGCAAACTGAATGACGATATGGCAACGTTAATTCTAAAAAACGTAACtat                                               >  1:122670/1‑141 (MQ=255)
                                                                                               aaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                             >  1:53650‑M1/23‑141 (MQ=255)
                                                                                               aaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                             >  1:24441‑M1/23‑141 (MQ=255)
                                                                                                 tACTAACTAAAATCAATTTTAATAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAATATTGAAACCTGAATGACAATCTGTCAACGTTAATAACAAAAAACGTAACTATAAGt                                           <  2:345255/141‑1 (MQ=37)
                                                                                                 tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                           <  2:55194/141‑1 (MQ=255)
                                                                                                 aaataaaataaaaaaaaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAACACTGAATGACAATATGTCAAAAATAAAAAAAAAAAACGTAACTATAAGt                                           >  2:223135‑M1/21‑141 (MQ=255)
                                                                                                  aataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                          >  1:342800‑M1/20‑141 (MQ=255)
                                                                                                   ataaaatggacaagataaAACAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAATCGTAACTATAAGTTa                                         >  1:282973‑M1/19‑141 (MQ=255)
                                                                                                   ataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGGAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                         <  1:289751‑M1/123‑1 (MQ=255)
                                                                                                   ataaaaaagacaagataaAAAAATTTATTGATTTAAATGTTGATAAAATGTATAATTAATTCCTGTCGGTAAGAAAAAAGAACATTGAAAACTGAATGACAATATGTCAACTTTAATTCCAAAAAACGTAACTACAAGTTa                                         >  2:267983‑M1/19‑141 (MQ=255)
                                                                                                      atatggacacgagaaAAAAAGTTATTTACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                      >  2:359785‑M1/16‑141 (MQ=255)
                                                                                                       aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                     >  1:302307‑M1/15‑141 (MQ=255)
                                                                                                        atggacaaggtaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                    <  2:56303‑M1/128‑1 (MQ=255)
                                                                                                        atggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                    <  1:20000‑M1/128‑1 (MQ=255)
                                                                                                            attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                <  2:42512/141‑1 (MQ=255)
                                                                                                                 ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGATAATATGTAAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                           <  2:49343/141‑1 (MQ=255)
                                                                                                                  tttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATCGAAAACTGCATGCCAATATTTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          <  1:166441/141‑1 (MQ=255)
                                                                                                                               tGAATTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa             <  2:120166/141‑1 (MQ=255)
                                                                                                                                gACTTAAATGTTGATAAAATGTATCATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt            <  1:137921/141‑1 (MQ=255)
                                                                                                                                gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt            <  1:313737/141‑1 (MQ=255)
                                                                                                                                  cTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCa          >  1:207496/1‑141 (MQ=255)
                                                                                                                                     aaaTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAc       <  1:342266/141‑1 (MQ=255)
                                                                                                                                       aTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat     <  2:207496/141‑1 (MQ=255)
                                                                                                                                        tGTTAATAAAATGTATAATTAATTCTTCTCGGTAAGAAAAATGAACATTAAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc    <  1:76441/141‑1 (MQ=255)
                                                                                                                                         gTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAAAAACATTATTTAGTATTTATGAGCTAATCAAAcagca   >  2:292909/1‑141 (MQ=255)
                                                                                                                                          ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat  <  1:279336/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512503‑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: