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

AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512503‑512783
                                                                                                                                            |                                                                                                                                            
aGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                                                                              >  1:226105/1‑141 (MQ=255)
    gtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAATTATTGACTTAAATGTtaataa                                                                                                                                          >  2:200141/1‑141 (MQ=255)
    gaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTtgataa                                                                                                                                          >  1:23551‑M1/114‑141 (MQ=255)
       aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  1:315719/1‑141 (MQ=255)
         taagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAGTTATTGACTTAAATGTTGATAAAATGta                                                                                                                                    >  2:279346‑M1/109‑141 (MQ=255)
          aagctacatagcgtagatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGta                                                                                                                                    >  2:317149‑M1/108‑141 (MQ=255)
               tgttgtTAAGCAATTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                                                                               <  2:34261/141‑1 (MQ=255)
               cgtgagcgtcgagggtagtcggacttagggttcgctagagtagaacgttgccaggaagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatt                                                                                                                               >  2:282103‑M1/103‑141 (MQ=255)
                gttgtaAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                              <  1:56020/141‑1 (MQ=255)
                  tagcgtcgatggtagtcgaacttacggtccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATGGACTTAAATGGTGATAAAATGTAtaattaat                                                                                                                            <  1:197670‑M1/42‑1 (MQ=255)
                   gttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                                                                           <  1:224136/141‑1 (MQ=255)
                    gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc                                                                                                                          <  1:31327‑M1/44‑1 (MQ=255)
                    gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGATATTGACTTAAATGTTGATAAACTGTGTAATTAATTc                                                                                                                          <  1:54437‑M1/44‑1 (MQ=255)
                    acgttgatggtagtcgaacttacgttccactagagtagaaagttgcaaggcagttttttaatcaaattttgtttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc                                                                                                                          >  2:174714‑M1/98‑141 (MQ=255)
                      gtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCtt                                                                                                                        <  2:321347‑M1/46‑1 (MQ=255)
                            ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  <  2:234847‑M1/52‑1 (MQ=255)
                            ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  <  1:334006‑M1/52‑1 (MQ=255)
                            ggtagacgaacttacgttcctcgagaggagaacattaccaggcagtttattaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  >  2:66500‑M1/90‑141 (MQ=255)
                                     cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                         >  1:382223/1‑141 (MQ=255)
                                     cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTAATTGACTAAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                         >  2:185860/1‑141 (MQ=255)
                                         acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                     <  1:42432‑M1/65‑1 (MQ=255)
                                          cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGCAAGAAAAATAAACa                                                                                                    <  2:141341‑M1/66‑1 (MQ=255)
                                          atgtacgatatagtagaaggatgcaaggcagttttttaaacaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                    >  2:304235‑M1/76‑141 (MQ=255)
                                          agttccgaaagagtagaacgtagaaaggcagttttttaatcaaattgtggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTCATTCATGTCGGTAAGAAAAATGAACa                                                                                                    >  2:402952‑M1/76‑141 (MQ=255)
                                           gttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                                                                   >  1:144664‑M1/75‑141 (MQ=255)
                                            ttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                  <  1:113782‑M1/68‑1 (MQ=255)
                                                gctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGCAAAATGAACATTGaaa                                                                                              <  1:45777‑M1/72‑1 (MQ=255)
                                                  tagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                            >  2:281819‑M1/68‑141 (MQ=255)
                                                             gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAat                                                                                 <  1:127581‑M1/85‑1 (MQ=255)
                                                                acaaaacaaCTCAACTTAAGAAATAACTTGACTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAACATATg                                                                              >  2:274312/1‑141 (MQ=255)
                                                                  caggcagttttttaatcaaattttggttaaaagataaaatggacaagataaAAAAAGTTATTGACTTAAATGTCGATAAAATGTATAATTAATTCTTGTCGTTAAGAAAAATGAACATTGAAAACTGAATGACACTATGTc                                                                            <  2:362686‑M1/90‑1 (MQ=255)
                                                                  caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                            <  1:11472‑M1/90‑1 (MQ=255)
                                                                   aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                           >  1:101827‑M1/51‑141 (MQ=255)
                                                                     acaacTCAACTTAAGAAATGACGTGAATTACTAACGAAAATTAATGTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTGGTCGGTAAGAAAAATGAACATTGAAAACCGAATGACAATAGATCaaa                                                                         >  2:325679/1‑140 (MQ=255)
                                                                        gttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                                                                      <  2:194592‑M1/96‑1 (MQ=255)
                                                                         ttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGAATTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTa                                                                     <  1:401822‑M1/97‑1 (MQ=255)
                                                                           taactTAAAAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                   <  2:95759/140‑1 (MQ=255)
                                                                            aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                  <  1:53816/141‑1 (MQ=255)
                                                                             ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                                 <  1:390921‑M1/101‑1 (MQ=255)
                                                                               aatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTtaattaca                                                               <  1:288154‑M1/103‑4 (MQ=255)
                                                                               aatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                               <  1:142255‑M1/103‑1 (MQ=255)
                                                                                tAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGTATGTCAAGTTTTCAACGTTAATTCCaa                                                              <  1:392028/141‑1 (MQ=255)
                                                                                tAAGAAATAACATGAATTAATAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                              <  2:113583/141‑1 (MQ=255)
                                                                                 acaatatttggttaaaaaataaaatcgaaaagataaAAATAGTTATTGTCTTAAATGTTGATAAAATGTATAATTAATTCTTTTCGGTAAGACAAATGAACCTTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                             >  2:111225‑M1/37‑141 (MQ=255)
                                                                                   aaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                           <  1:372238‑M1/107‑1 (MQ=255)
                                                                                                ttACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAg                                              <  1:44214/141‑1 (MQ=255)
                                                                                                  aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCTGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACttatagtt                                            <  2:169380/141‑8 (MQ=255)
                                                                                                    tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAAATTGAAAACTGAATGACAATATGTCAAAGTTTATTCTAAAAAACGTAAATATAAGTTaa                                          >  1:111173/1‑140 (MQ=38)
                                                                                                     aatgAAAATTAATTTTAAAAAGTTATTTACTTAAATGTAAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         <  2:20495/138‑1 (MQ=255)
                                                                                                     aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         >  2:103555‑M1/17‑141 (MQ=255)
                                                                                                     aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTCATTCCGTAAATCGTAACTATAAGTTACa                                         <  1:324124/141‑1 (MQ=255)
                                                                                                     aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         <  2:220541/141‑1 (MQ=255)
                                                                                                       cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  1:8725/1‑141 (MQ=255)
                                                                                                       aatgaacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  2:67059‑M1/15‑141 (MQ=255)
                                                                                                        atggacaagataaAAAAAGTTATTTACTTAAATGTTGATAAAATGTATAATTAATTTTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                      >  2:137244‑M1/14‑141 (MQ=255)
                                                                                                         aaaaTTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGACCATTGAAATCTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                     <  1:28698/141‑1 (MQ=255)
                                                                                                               tattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               >  1:316570/2‑141 (MQ=255)
                                                                                                               aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACGTTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               >  1:219673/1‑141 (MQ=255)
                                                                                                               aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                               <  1:102018/141‑1 (MQ=255)
                                                                                                                   ttAAAAATTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  2:226105/141‑1 (MQ=255)
                                                                                                                   ttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  1:162437/141‑1 (MQ=255)
                                                                                                                    aAAAAAGTTATTGAGTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAATGAATGAAAATATGTCAACATTAATTCCAAAAAACGGAACTATAAACTACAAAAATTATTTAGTAt                          <  2:369536‑M1/140‑1 (MQ=255)
                                                                                                                      aaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACGCAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt                        >  1:237702/1‑141 (MQ=255)
                                                                                                                         aGTTATTGACTTAAATGGTGATAAAATGTATAATTAATTTTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTATTATTTATg                     >  1:60677/1‑141 (MQ=255)
                                                                                                                          gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                    >  2:121349/1‑141 (MQ=255)
                                                                                                                            tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                  <  2:399920/141‑1 (MQ=255)
                                                                                                                            tATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGGAAATTTAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                  >  1:368844/1‑141 (MQ=255)
                                                                                                                                       aTGTTGATAAAATGTATAATTAATTCTTGTCGGTCAGAAAAATAAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAAAAACATTATTTAGTACTAATTATCTAATAAAAcat       >  2:246691/1‑141 (MQ=18)
                                                                                                                                        tGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc      <  2:11472/141‑1 (MQ=255)
                                                                                                                                            aataaAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGGTACAAAACTTATTTAGTATTTATGAGCTAATCAAACATCATaa  >  1:151014/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512503‑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: